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Englisch-französische Titelerläuterungen der in Band XIX, Heft 1 
der „LAMM“ (1939) veröffentlichten Hauptaufsätze. | 








A. Huber in Wien. 


The progress of the melting point in a linear conductor. The term „melting curve“ 
is used to designate the graphical representation of the dependence upon time of that part 


of a linear conductor, where the solid phase passes over into the liquid. The melting 
curve is approximated in numerous small intervals by straight lines. Its determination _ 


follows then from the solution of boundary problems for a series of trapezoidal areas for 
which purpose a special case of transformation, discovered by Appell, is used. 


Sur la propagation de la fusion dans un conducteur linsalre. On appelle „courbe de 
fusion“ la repr6sentation graphique de la loi, suivant laquelle le döplacement du point de 
-contiguit& des phases liquide et solide d&pend du temps. Pour aborder le problöme, l’auteur 
remplace la courbe,. dans des intervalles suffisamment petits, par un trait polygonal. De 
cette manidre la determination de la courbe se ramene & la r&solution des probl&mes de 
limite pour une suite de domaines de la forme d’un trapöze en employant un cas particulier 
d’une transformation indiquse par Appell. 


W. Bollay in Cambridge, Mass. U. $. A. 
A non-linear wing theory and its a lication to rectangular wings of small aspect 


ratio. The author develops a non-linear wing theory for calculating the forces acting. 


x the wing of small aspect ratio. This theory differs prineipally from the usual theory 
of lifting vortex line or lifting surface, in. that the trailing vortices leave the lifting 
surface at an angle, at the point of formation. The ensuing velocity tangential to the 
lifting surface, alters the resulting force on the surface. The theory agrees well with the 
experimental results obtained by Winter. f 


Une th6orie non-linsaire des alles et son application aux alles rectangulalires & 
grande profondeur. Pour döterminer les forces agissant sur une aile & grande profondeur, 
l'auteur d&veloppe une thöorie qui s’appuie sur la supposition que les tourbillons se d&tachent 
de la surface portante sous un angle & l’endroit m&me oü ils se forment. En ce cas nait 
une vitesse tangentielle & la surface portante modifiant la force rösultante qui agit sur 
l’aile. Les r&sultats de la thöorie s’accordent bien avec les exp6riences de Winter. 


R. Müller in Berlin. 


The numerical computation and application of certain functions related to Bessel’s 
together with comments on Bessel functions themselves. With the aid of the 
tables given for the expressions: 


> je x | a2 
Jrunas, [zmwanae, = fermumat Br 


a large number of functions in integral form become readily accessible; for example the 
Struve functions for half integers as parameters. Further these tables enable the solution 
of er gear Bessel differential equations with a SURSETORE function, in the form of 
a potential series with whole positive and negative exponents. renders the numerical 
treatment of wing vibrations and induced vibrations of cylindrical bodies simpler. Finally, 
a number of differential equations reducible to Bessel equations are given. 


Sur l’6valuation num6rique et l’application de certaines fonctions ressemblant 
aux fonctions besseliennes, avec quelques annotations sur le domeaine des 


fonctions bess6liennes. Ä l’aide des tables pour les fonctions 


x x 2 | 
[»unas, fızmwunas, et Kae et. 43 
0 0 


il est possible de calculer un assez grand nombre d’autres fonctions se pr&sentant sous les 
formes d’intögrales,’p. e. les fonctions de Struve aux paramötres de la forme + 3 ..En 


outre les tables permettent de r&soudre l’&quation non-homogöne de Bessel contenant une 
fonction perturbatrice par une s6rie de puissances ‘a &xposants positifs et nögatifs. Cette 
solution faeilite l’6tude des oscillations des ailes ainsi que des oseillations contraintes des 
corps eylindriques. L’auteur finit par indiquer quelques öquations difförentielles r&ductibles 
aux öquations diff6rentielles de Bessel. 
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HAUPTAUFSÄTZE 
Über das Fortschreiten der Schmelzgrenze in einem linearen 
Leiter. 


Von A. Huber in Freiburg (Schweiz), jetzt in Wien. 


W sich ein linearer homogener Leiter teils im festen, teils im flüssigen Aggregat- 
zustande befindet, dann werden die Trennungsstellen der beiden Zustände, an denen 
also gerade die Schmelztemperatur herrscht, im allgemeinen ihre Lage mit der Zeit ändern. 
In einem ., #-Diagramm werden daher jene Stellen eine aus einem oder mehreren getrennten 
Stücken bestehende Kurve beschreiben, die wir als Schmelzkurve') bezeiehnen wollen. Die 
besondere Schwierigkeit der zur Ermittlung dieser Kurve zu lösenden Randwertaufgabe besteht 
bekanntlich darın?), daß ein Teil des Randes selber zu den Unbekannten des Problemes gehört. 
Da mir abgesehen von den in der Enzyklopädie d. mathem. Wissensch., 5. Bd., 1. H., S. 204 f. 
angeführten älteren Lösungsversuchen eines besonderen Falles dieser Aufgabe aus neuerer 
Zeit nur eine Arbeit des Herrn K. Lachmann’) und ein Referat über einen Aufsatz des 
Herrn S. Kowner') bekannt geworden sind, so möchte ich im folgenden eine allgemein an- 
wendbare neue Methode zur näherungsweisen Bestimmung der Schmelzkurve mitteilen und 
durch ein vollständig durcehgerechnetes Beispiel illustrieren. Ihre praktische Anwendung er- 
fordert allerdings, wenn man auf eine bei Wärmeleitungsproblemen ungewöhnlich große 
(Genauigkeit Wert legt, einen erheblichen Aufwand an rechnerischer Arbeit, aber auch die 
von Herrn Lachmann angegebene Lösung des von ihm betrachteten teehnisch interessanten 
Sonderfalles, der sich dureh besondere Symmetrieverhältnisse auszeichnet, stellt in dieser 
Hinsicht ziemlich hohe Ansprüche. Weiters liegt dem von Herrn Lachmann mittels der 
Methode der Quellpunkte durchgeführten Ansatze, der sich überhaupt nur -in dem von ihm 
betrachteten besonderen Fall machen läßt, die wesentliche Annahme zugrunde, daß Temperatur- 
und Wärmeleitzahl vom Aggregatzustande unabhängig seien, so daß man für die ganze Aus- 
dehnung des Leiters nur eine einzige Differentialgleichung zu integrieren hat. Diese Annahme 
hat natürlich eine bedeutende Vereinfachung der Rechnung zur Folge, indessen dürften wohl 
nicht unbegründete Bedenken am Platze sein, ob sie in der Regel wenigstens einigermaßen 
den Beobachtungstatsachen entspricht. Nun zeigt aber schon eine flüchtige Betrachtung der 
wenigen Fälle, in denen überhaupt Beobachtungen der Wärmeleitzahl z aus einem Temperatur- 

1) Obwohl diese Bezeichnung in der Thermodynamik bekanntlich in einem ganz anderen Sinne bereits ver- 
wendet wird, wollen wir sie der Kürze wegen doch beibehalten, da ja hier wohl keine Verweehslung zu befürehten ist. 

®2\, Vel. z.B. Riemann-Weber: Die Diff.- u. Int.-Gl.d. Mech. u. Phys., 7. Aufl., hgg. v.P.Frank u.R.v. Mises 
(1927), 2. Bd., S. 20 If, 

3) ZAMM 15 (1935), S. 345 bis 358: Zum Problem des Erstarrens f. d. durch zwei parallele Ebenen begrenzten 
Körper. 

4) Journ. Geophys. Moskow 3 (1935), S. 32 bis 41 [Jahrb. d. Fortschr. d. Math. 59/11 (1933), S. 1477]. 
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intervall vorliegen, das auch den Schmelzpunkt der betreffenden Substanz enthält?’), daß siclı 
x» ım allgemeinen zwar wenig ändert, solange die Temperatur entweder ständig ober- oder 
ständig unterhalb des Schmelzpunktes bleibt, daß dagegen beim Überschreiten desselben eine 
oft bis über den halben Betrag anwachsende plötzliche Abnahme von x stattfindet. Nur 
Wismut zeiet eine starke Zunahme, indem # (200°) = 0,017 und x (300°) — 0,037 werden. 
Bei dieser Sachlage dürfte daher eine Näherungsmethode zur Bestimmung der Schmelzkurve, 
die der raschen Veränderliehkeit von x in der Umgebung des Schmelzpunktes wenigstens 
dadurch Reehnung trägt, daß den beiden Aggregatzuständen entsprechend zwei verschiedene 
konstante Mittelwerte für z gesetzt werden, nicht ohne Interesse sein, wenn auch ihre Durch- 
führung in einem konkreten Fall ziemlich viel Arbeit erfordert. 


I. Ansatz zur strengen Lösung für den unbegrenzten Leiter. Um die großen Schwierig- 
keiten, die einer exakten Lösung des in Frage stehenden Problemes im Wege stehen, deutlich 
hervortreten zu lassen, wollen wir für den Fall des beiderseits unbegrenzten Leiters zeigen, 
wie sie sich bei Benutzung der bekannten Ansätze gestalten würde. Durch die Indizes 1 
und 2 unterscheiden wir den festen und den flüssigen Zustand und nehmen an, daß der Leiter 
zur Zeit t=0 für 2<V fest und für 2>0 Nlüssıe sei. Dann handelt es sich um die Be- 
stimmung je eines Integrales der beiden Gleichungen 

04; 7, 


dA; 


— 4=132), 
of ’ 0) 30" 


die den folgenden Bedingungen genügen: 


a)d (er, Ve, le) 0 für O0 und d,(e, 0) p,(e) 0 für z=0, wobei p,(0)=g,(0)=0 
die Schmelztemperatur bedeuten und 9,(a#), 9,(#) stetig differenzierbar sein sollen, und es 
möge sich ein eventuell ursprünglich an der Trennungsstelle vorhanden gewesener Temperatur- 
sprung bereits ausgeglichen haben. 


b) Für 720 seien mit === s(f) die noch unbekannte Gleichung der Schmelzkurve 
bezeichnet, mit #; die als konstant angenommenen Wärmeleitzahlen, mit o die Dichte und 
mit / die Schmelzwärme, dann soll: 


a) YArZUmIE KATI M IE Tim 

Or, Q | i ds ®) f 

/ Hr” Ko A° . . 
P) | ) V.r ©; Or ® {N O dt | 


Ks seien nun We, N zwei den Anfanesbedineungen a) eenürende Inteerale und 
{ h few) tes) i 


! 


J; Is (hd, #] vi), 
dann haben wir je ein weiteres Integral 9; [.r, f; s(D] so zu bestimmen, daß 8 
9, [2,0;:s(d]=0 und 9’lsd,t,sd)=y;(i). ! 
Dazu setzt man bekanntlich 
/ e 
| (7) Ix—s(n]? 
[5 7% m 
gr ap 
); le. t: sh] - u ur u ze 7 
d;] N _ (f ı ud 
u 
woraus sieh zur Bestimmung der beiden Funktionen @;(r) die folgenden Integralgleichungen 
ergeben, deren Kerne aber von der ja noch unbekannten Funktion s (f) abhängen: 
ö Is (f) s(r)]? 
| s(f) — s{r) ; 
ta;2(t —ı 1} 
| [): od; if) Pa Pr 319 .e& N .dtr=wy;f(t. 
2a yn. (E — 7)" ? 
v 3 
& 


Für die Temperaturverteilung in den beiden Teilen des Leiters erhält man somit 
la, t; s dl= 9,’ (a, N + I, le,.t;s(d], 
wobei identisch für jedes s (f): 


9;,[2,0:s(N]=g;(&) und 9[s(d, t:;s(d]=0, 

5) ,andolt-Börnstein: Phys.-chem. Tabellen, 5. Aufl., 2. Bd. (1923), S. 1289. 

6, db man hier für © die Diehte des festen oder des flüssigen Zustandes zu wählen hat, hängt, strenge genommen, 
von der Riehtung ab, nach der die Sehmelzgrenze fortsehreitet. Man müßte aber dann, wenn man diesen Umstand 
hericksiehtigen wollte, folgeriehtig aueh sehon bei der Bildung der Wärmeleitungsgleichung die Wärmeausdehnung 
des Leiters beachten. Es soll also o einen Mittelwert der im allgemeinen ohnehin nieht stark mit der Temperatur 
variierenden Dichte bedeuten. 
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während sich zur Bestimmung der Schmelzkurve &=s (f) die gewöhnliche Differentialgleichung 
erster Ordnung (1) ergibt. Um mit diesem Ansatze s (ft) näherungsweise bestimmen zu können. 
wäre es wohl das Nächstliegende, s(f) durch ein Polygon zu approximieren, aber die dann 
nicht zu vermeidende Lösung der beiden Integralgzleichungen würde dabei erhebliche Schwierig- 
keiten verursachen, die sich für einen beiderseits begrenzten Leiter noch bedeutend steigern 
würden. Deshalb war es unser Ziel, die nun zu schildernde Näherungsmethode zur Bestimmung 
der Schmelzkurve von vorneherein schon so einzurichten, daß dabei wenigstens explizite 
die Auflösung von Integralgleichungen umgangen wird. 


2. Ansatz zur näherungsweisen Bestimmung der Schmelzkurve. Um unnötige Weitläufig- 
keiten zu ersparen, wollen wir nun unsere Näherungsmethode an einem besonderen Fall aus- 
einandersetzen, der aber doch wieder so allgemein ist, daß ihre wesentlichen Züge deutlich 
hervortreten. Wir betrachten dazu einen linearen Leiter, dessen Enden die Abszissen £ —0 
und «=1 haben mögen und der zur Zeit £=0 im Intervall (0, 4,) fest, in (4,,9,) Nüssig und 
in (9, 1) wieder fest sein soll. Natürlich setzen wir dabei voraus, daß sich auch im flüssigen 
Teil die Wärme nur durch Leitung und nieht auch durch Konvektion fortpflanzen soll und 
daß die Wärme- und die Temperaturleitzahlen für jeden der beiden Aggregatzustände konstant 
bleibende Werte haben mögen. Bezeichnen wir also die Schmelztemperatur wieder mit 9=0, 
dann soll: 

(a) <0 für 0<e: 1 


(x, 0) DE ee ; 5. ers TR 
Eh Me at 


wobeı die Funktionen y;(x) als stetig differenzierbar vorausgesetzt werden und an den Inter- 
vallgrenzen stetig ineinander übergehen sollen, also 9, (q,)= 9, (y,)= 9, (9)=95(9%)=0. An 
den Enden z=0 und ===! sollen die Temperaturen etwa als Funktionen der Zeit vor- 
eeschrieben sein: 


‚en u Peru : .  errerererı i 


wobei noch 9,(0) = f(0), 9, ()=g(0) und der Einfachheit halber stets fi) <O und g(d)<O 
bleiben sollen, damit wir uns bei der Bestimmung der Schmelzkurve vorläufig auf einen die 
Punkte #2—=g, und @=q, verbindenden Bogen beschränken können. Man wird später leicht 
bemerken, daß diese Einschränkungen keineswegs wesentlich sind, sowie auch, daß man die 
Randbedingungen (3), ohne die Anwendbarkeit der Methode zu beeinträchtigen, durch andere 
ersetzen könnte. 

Die Schmelzkurve geht also jedenfalls dureh die Punkte #=—=qg, und ===, hindureh, 
und die Richtungen ihrer Tangenten in diesen Punkten ergeben sich sofort aus (1), die ja 
auch für t=0 gelten muß: 


i 2 a ds\ 
KH Po (I) TR" P, )=4o|,, 
/’4ı 
(4). 
’ r . ds 
A, "Ya (4) IE (45) _ A:0' it 
N 42 





E a . (ds h 
Bei unseren Annahmen über die Funktionen 9;(#) werden die 3 ’ sicherlich endliche 
\Ar)qi 


Größen sein, die eventuell auch verschwinden können, so daß also die Schmelzkurve die x-Achse 
im allgemeinen unter einem von Null verschiedenen Winkel schneiden wird. Der Fall, daß die 
Schmelzkurve die &-Achse berührt, was dann eintreten würde, wenn die eine oder die andere 
der auf den linken Seiten von (4) stehenden Ableitungen unendlich werden würde, bedarf einer 
besonderen Erörterung, weil ja dann die Konstruktion der Schmelzkurve gleich zu Beginn 
versagen würde. Dies wird sich praktisch sogar sehr häufig einstellen, nämlich immer dann, 
wenn der feste Teil des Leiters mit dem über dem Schmelzpunkt erwärmten flüssigen Teil 
plötzlich zusammenstößt, so daß im Augenblicke der Berührung eine Unstetigkeit der Temperatur 
selber an der Berührungsstelle auftritt, doch soll vorläufig dieser Fall bei den folgenden Über- 
legungen ausdrücklich ausgeschlossen sein, was Ja auch schon in den Annahmen über die 
Funktionen 9;(x#) zum Ausdruck gebracht worden ist. 

Wir denken uns nun in den Punkten 2<=g, und z=g, die Tangenten der Schmelz- 
kurve konstruiert, deren Richtungen ja dureh (4) gegeben sind, und schneiden sie mit einer 
Parallelen zur x-Achse im Abstande At>0O in den Punkten 9,” bzw. 9,. Wenn wir sodann 
unter der Annahme, daß die Schmelzkurve für ein hinreichend kleines At durch die soeben 
konstruierten Tangenten genügend genau approximiert wird, die Temperaturverteilung im 

je 
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eanzen Leiter, also von z2=0 bis z==l, für { I! den gegebenen Anfanes- und Rand- 

bedingungen gemäß berechnen können, dann lassen sich wieder mittels der jetzt aber für 
ds " (ds ’ - \ fl 

EI} er . At bzw. für eg, +! 4) - It analog wie (4) zu bildenden Relationen die 
( 4 y i 12 


Tangenten in den Punkten 9,’ und @, konstruieren. Man kann also durch wiederholte An- 
wendung dieses Verfahrens die gesuchte Schmelzkurve durch ein Polygon um so genauer 
approximieren, je kleiner man die Zeitintervalle 4 wählt. 


04 ‚0; ; \ 
Wegen | N 1 won kann man mittels der aus (1) folgenden Relation: 
VE Oi; 
(2 J, ds N, | k u ), ds 0° ), ) ß d’s 
Be, 3 a * m Ha es pe | nee“ - 
o XL ft Vr of en Q d ! Or of x s(f) zu d f° 


auch leicht eine parabolische Approximation der Schmelzkurve erreichen. Wir werden in der 
Tat später sehen, daß die Bildung der Werte der Ableitungen höherer Ordnung von 9 nach 
e allein, auf die es Ja letzten Endes ankommt, auch praktisch leicht durchgeführt werden 
kann, so daß man für die Funktionen s«4;(f), oder besser noch für ihre Umkehrungen auch 
in einem größeren Intervall brauchbare Darstellungen finden könnte. Diese Überlegungen 
ließen sich übrigens ohne erhebliche Schwierigkeiten zu einem „Existenzbeweis* für die 
Funktionen s,(M) ausbauen, doch ist hier nicht der passende Ort hierzu. Wir werden aber 
später, sobald wir die hierzu nötigen Formeln aufgestellt haben werden, leicht zeigen können, 
daß unser Verfahren das sehr vorteilhafte Bestreben besitzt, die infolge der Approximation 
der Schmelzkurve durch ihre Tangenten entstehenden Fehler von selber auszugleichen, so daß 
man praktisch wohl in den meisten Fällen mit einer Polygonapproximation der Schmelzkurve 
auskommen wird, zumal da es ja angesichts der großen Unsicherheit, mit der die Angaben 
der hier in Frage kommenden thermischen Konstanten naturgemäß behaftet sind, überhaupt 
nicht viel Sinn hat, an die Genauigkeit der Reehnung hohe Ansprüche zu stellen. 

Unser Problem erscheint somit auf die wiederholte Lösung der folgenden Randwert- 
aufgabe zurückgeführt: es ist der räumlich-zeitliche Temperaturverlauf in einem Leiter zu 
bestimmen, dessen Anfangstemperatur als differenzierbare Funktion seiner Abszisse gegeben 
ist und von dem beide Enden auf der Temperatur J=0 gehalten werden und sich mit 
bekannten konstanten Geschwindigekeiten bewegen. Es seı hierzu noch bemerkt, daß man 
die für die beweglichen Enden gemachten Annahmen #F=0O auch durch allgemeinere ersetzen 
kann, ohne die Anwendbarkeit der nun zu entwickelnden Methode zu beeinträchtigen, so daß 
sie auch zur näherungsweisen Lösung von Randwertaufgaben für ein Gebiet dienen kann, das 
unten von einer zur .«-Achse parallelen Strecke und seitwärts von zwei Kurven begrenzt wird, 
die von einer Parallelen zur x-Achse in je einem Punkte bloß getroffen werden’). Unsere 
Methode versagt natürlich nieht nur in dem bereits erwähnten Falle eines Temperatursprunges 
an der Schmelzgrenze, sondern auch dann, wenn der Schmelz- bzw. Erstarrungsvorgang an 
einem Ende des Leiters beginnt, da ja dann die Relation (4) nicht verwendet werden kann. 
Wir werden später noch ausführlich darauf zurückkommen und wenden uns jetzt zur 


3. Lösung der Randwertauigabe für einen Leiter, dessen Länge eine lineare Funktion der 
Zeit ist. Wir bezeichnen jetzt mit X die Abszisse eines beliebigen Punktes unseres Leiters, 
mit 7 die Zeit und mit 
(); Pi "+ fi 5 A a Se 
die Abszissen der beweglichen Enden. Dann soll dasjenige Integral der Gleichung 
0,09 
ns — 
O1 0X? 
gesucht werden, das den folgenden Bedingungen genügt: 
a) (X, V)=PI(X) für ZSÄSQ@ | 
i | | (6). 
b» (pP T+4.NM=dV, wobeı Pig,) Diqg)—=0)| 
Um das Ergebnis gleich in einer für die numerische Auswertung bequemen Form zu erhalten, 
führen wir dimensionslose Variable ein, indem wir setzen: 
N ( 7 u : na 
am hi : f .. J(X,T) (ae. : D(X) D (9, 7 (9 4,) ge) == ya) (D). 
1. =4: 1. = i 


') Vortrag bei der Tagung der Deutschen Mathematiker-Vereinigung in Bad Kreuznach am 21. Sept. 1957. 
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Ks ist also nun jenes Integral der Gleichung 


RE), 0’u 


Of 0x” 
zu suchen, für welches: 
a) Kr, O)=ypla) für Ozr—l: b)Y(pyt.D=dl(p tl, O0 167, 
’ . 4(.— 1 
wobeı noch p; . : . -p; gesetzt wurde. 
X 


Wir zerleren die gesuchte Lösung in der üblichen Weise, indem wir durch den Ansatz 
) "0", wo 


Du (X, !) 


die Anfangstemperatur zunächst auf Null reduzieren. Wegen 9 O)=p()=V ıst für 0 
N ij! 


zwar noch 9 selber, aber im allgemeinen nicht mehr auch bei 20 und x = 1 stetig, 


en 
0) A 
doch ist dieser Nachteil, der sich übrigens leicht beheben ließe, für unsere Zwecke belanglos. 
Nun ist 9” so zu bestimmen, daß 9’ (#,0) =0 und 


(pt) (pp, tdey (dd) und 9"’(kpt+Hh,de - 9 (pt+l,d=y,(th) (6°), 


und diese Aufgabe läßt sich, wie wir gleich sehen werden, durch eine einfache Transformation 
auf eine analoge für einen Leiter mit festen Enden zurückführen. 

Schon vor längerer Zeit hat nämlich Appell’) gezeigt, daß durch die folgende Trans- 
formation die Wärmeleitungseleichung in sieh übergelührt wird: 

> Bi) (' (A P)? 
u Pr „r ‚ ii ‚ gr tt oe - 3 ry = 
N y a: u ! = u ı) = :6& (ce, tt) ‚ . (4). 

t—a '' T—a vi—a 
Die Transformation der beiden unabhängig Veränderlichen bedeutet eine Kollineation, bei der 
die Gerade t=a der (#,#)-Ebene in die unendlich ferne Gerade der (x’ f)-Ebene übergeht. 
Nun sieht man aber leicht, daß man über die Parameter a,«', P,p' und y so verfügen kann, 
daß diese Kollineation eine zentrische wird mit der #-Achse als Kollineationsachse, wobei 
noch die beiden Geraden z=p,-t und z2=p,-t-+1 bzw. in die Geraden «’=0 und x’ =] 
übergehen. Damit nämlich zunächst die »-Achse zur Kollineationsachse wird, muß 0 
aus t=0 folgen, also muß gemäß der zweiten Gl. (7) 


> ’ 


Au = A . «A ö . ° . . . . . . . . . . (S), 


‘ 


und die erste Gl. (7) für x’ = x identisch in = erfüllt sein: 
(Yy + a) yD Da u. 


woraus wegen (8) folgt: 


so daß auch 


Damit nun noch z= p, tin @=0 und »—py t-+1 in „«’=1 übergehen, müssen die aus der 
ersten Gl. (7) sich ergebenden Identitäten in f bestehen: 


BP-t—a)—a-(p -t—- B)=0 und (P-1)-(t-a)-a-(py-t+1-P) 0. 
Dazu muß offenbar 5=ap, und % -1l=a:p,, womit nach (S°) und (8°) folgt, daß: 


Pi | | 
P=P et a a=—y P ; RILTO BEN EEK OASE 
Pi Ps» 
Man rechnet leicht nach, daß « und 5 Ordinate und Abszisse des Schnittpunktes der beiden 
Geraden x =p, -t und «= p, f-+H1 bedeuten, in dessen senkreehte Projektion auf die «-Achse 
das Kollineationszentrum fällt, weil ja aus der zweiten Gl. (7) folgt, daß !=«"= —- a für 
0?2 82 


8, Journ. de math. (4), 8 (1892), S. 187 bis 216, Sur l’&quation _ et la Theorie de la ehaleur, 
8x2 OY 
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=, Mit den ın (9%) angegebenen Werten für « und 5 lassen sich nun die ersten zwei 


(il. (7) so sehreiben: 


3 aAvx h- / y n+ f | 
4 in = ; ; . ö . 0 
A ! A / ); 
oder 
ıL° r’ t— )* rt’ (1 f’ 
7: n / er R ’ } i : . . . u" . 
I / (dl { (1 ) 


Die für 9” vorgeschriebenen Anfanes- und Randwerte erhalten schließlich zemäß 
der dritten Gl. (7) für Pt’, #’) die foleende Form, wenn wir noch das ohnehin herausfallende 





E | setzen: 
ir. VO) (0) für 0) r" | 
f’ a“ tt’ | af f’ 
JO.) | PL; .y 5 ME 
a1 "\a+r) A ie .°° 
(] h) 
| ‚ (ı” t ıf’ ) | af j 
| Ei | - { ! - ( } 
v1 ' A ! | L 
Damit ergibt sich 9 (.@’, F) als Summe der beiden „Faltungsintegrale*“ ?): 
Be, Darth ei ihre 
wobeı 
(2 'ın) 
_ (fr = x’ 12 nn) 
inte A, Wh \ | 4m 
Gie,®#) z 7% (ve +2n).e 
GR ee Vıt BUNT” aa 
S FT 


Für die Anwendune unserer Methode ist es nun von eroßer Wiehtiekeit. daß man die 
in (12) vorkommenden Doppelintegrale auf einfache reduzieren kann, indem sich nämlich die 
Inteeration nach der Zeit elementar ausführen läßt. Es wird nämlich zunächst nach (6”\). 


| pt s)? ) (pDer’t+ll S))” 
tt h “ tt ) 
m / \. - 3 > r £ | ’ / - - I p’ > . 1 
ur If) Ple)-c ad& 111€ us. IT) Die)» € Gun, 
2yat\‘ >2yat)9 


ı ı 


also nach (9). (10°) und (11): 


In entsprechender Weise ergibt sich, wenn man 5 >(1- pP) und 2>(1 —£) mit Ausnahme 
des Arzumentes von « 


f 


Zur Bildung des ersten Faltungsintegrales in (12) setzen wir noch der Kürze wegen y„=."+2n 


und finden 
f’ 
RE > Dr 


FH 
£ 4 S 
Y Une er 
, = MA, L ); f’ „ ” 
y„.sirie) .Pp . \ Da z FR ERS ee -dE, 
FT. 


D".yi T . 
da ja dıe zur Vertaüschung von Summation und Integration nötigen Voraussetzungen erfüllt 
sind. Daraus entsteht nach für »„==0 offenbar erlaubter Vertauschung der Integrationsfolge 
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Führen wir im zweiten Inteeral statt r die neue Integrationsvariable o "_ ein, so wird 


nach einfacher Rechnung '") 


t’ % 
Yu“ L&2 /* Yn?® ] 2 ] ( Yı 2)2 
„ f (fr E y fr ‘(t6Ö T fr 
\ —e le er sen (n)-2 oe 
f "Yale / 
0 () 
wobei sgn (n)= + 1, je nachdem » =O und sgn (O)=1. Es wird somit 
| 
| 4 , (> S)* (Yua'+>)° 
A v ’ 
* L ; he 1 / £ ., 
Yı* Ga (#’) z ur de \ sen (n)\p (de dz 5 5 
> 7 
N F u 
und entsprechend, wenn y„=1-—.x"+2n gesetzt wird: 
| 
| - (3 =)? ( YUn +1 s)* 
4 Y 
j ‚ ' \ hi t 1’ ‘ y‚ 
JG 4) > or. \ sgn (m)\ 4 (£).e FE u; ' 
= 2 u . 
[71 FT. () 
v .. .. fi dj! art _ . A j A 
Damit erhält man für 9 — 9 + 9’ weeen (7) und (12) und weil jr 
y ’ f’ | 
WIE. } te) 
Zynt, ! 
ı 
| (X s)? X 7) 2 (3 &)° 4% | (Yn L. £)2 ( Yı 3 =)? 
4 t a) & , 2 ‚ 
Ye t po ) ! , \ sen (n)|e hi Le It 


| de, 
a  r 
[7 F- 


wo die rechts vorkommenden x und f noch nach (10°) dureh x’ und FF auszudrücken sind. 
Führt man dies aus, indem man besonders die folgenden Relationen 


l H | on p ae P. / y’ 3 
r f’ ns / y ; (a )‚ırtie)=Ma 


beachtet, so erhält man nach einigen Umformungen die für die numerische Auswertung wohl 
bequemste Darstellung, wenn man noch setzt 


yiä)-e di wi£): 
| 
(x’ 3,2 pr 
t(@ +’) Be; (2n+x’— 5)? 2n+x’+s)? 
e Y 
( “ ir ir . 
la’, t)- i ER \ le P BEE er 
aynt i 
nz 
0 n % 


verfügt, die für genügend kleine Schritte des Argumentes f berechnet sind, kann man davon 
zur Auswertung von (14) vorteilhaft Gebrauch machen, doch kommt man auch mit einer um- 
fangreicheren Tafel für e”* aus, da man nur wenige Glieder der in (14) vorkommenden Summe 
benötigt, indem ja das !’ in der Regel sehr klein sein wird. Bei Beachtung der Formeln (10°) 
sieht man sofort, daß für «>x und endliches 5 die Darstellung (14) übergeht in die für die 
Temperaturverteilung eines Leiters mit festen Enden. Die Lösung (14) ist natürlich nicht 
nur für den mittleren flüssigen Teil unseres Leiters brauchbar, sondern nach entsprechender 
Berücksichtieung der Randwerte auch für die beiden äußeren festen. Für den linken Teil 


| 
wird nämlich p, =g,==0, also auch p,’==0 und daher nach (9): A=0 und «a a lÜbenso 
wird für den rechten Teil p,=0 und 9, =1, also auch p,’ = 0 und daher 5=1 und a 2, 


10, Vgl. z. B. Meyver-Dirichlet: Bestimmte Integrale, Leipzig (1871), 3 99. 
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Diese Überleruneen bedürfen jedoeh noch einer Ergänzung in dem bisher stillschweigend 
ausgeschlossenen Fall, daß p, == p,, wo ja die oben betrachtete Kollineation in eine Affinität 
übergeht. Da wir die dann sich ergebenden Formeln für den unendlichen Leiter später be- 
nötieen werden, wollen wir kurz die Grenzübereänee andeuten, durch die man aus (9) die 
den Formeln (7) entsprechenden erhält. Es wırd zunächst: 





r ) P I’ / A r 
lım Pi 1: z I p y i 
m’>) d ( z 
wobei p,=ps° p. Damit wird (107): 
(. ( . Ze 
x Hr p' / E 55 I: er I: pP E: / f ; i r ; ö ; i . (15 ). 


Um schließlich die der dritten Gl. (7) ın unserem Falle entsprechende zu bekommen, setzen 
wir zunächst 


Wird nun |a sehr groß, dann ergibt sich für den Exponenten der e-Potenz in (7): 


" PL > 2 f 2 A? n 2 .3 
(.r » | & .» [1 n | ’ P By Bi A 
7 f\ | 7 7 A MM a” lu la \ 7 p 
n 
und somit weren Im p 
‘1 
la p’t) 
i,P Pr ( ’ - r 
J (ar. e t IX p nn, j i BER A i ee (15 ), 


was man natürlich aueh dureh einen naheliegenden Ansatz direkt hätte finden können. Die 
dureh (15°) dargestellte perspektive Affinität bedeutet offenbar eine Scherung parallel zur 
r-Achse., die selber punktweise festbleıbt. 


4. Der Beginn der Schmelzkurve bei stetiger Randbedingung. Die soeben auseinander- 
vesetzte Methode zur Bestimmung der Schmelzkurve ist natürlich nur dann anwendbar, wenn 
der Erstarrungsvorgang schon ins Innere des Leiters fortgeschritten ist. Wenn aber etwa 
ein flüssiger Leiter an einem seiner freien Enden zu erstarren beginnt, dann kann man die 
(4. (1) zur Berechnung der Tangentenrichtung der Schmelzkurve an diesem Ende des Leiters 
offenbar nieht mehr verwenden und muß daher zu diesem Zwecke andere Überlegungen an- 
stellen, die, je nachdem die für das betreffende Leiterende vorgeschriebene Randbedingung 
stetig oder unstetig ist, wesentlich verschieden ausfallen. 

Wir beschäftigen uns zunächst mit dem ersten Fall und betrachten einen flüssigen Leiter, 
dessen Enden stetie abzekühlt werden sollen. Da es sıch für uns nur darum handelt, die 
Richtung der Tangente an die Schmelzkurve in ihren Endpunkten zu bestimmen, so genügt 
es, wenn wir uns auf den mit geringerem Rechenaufwand zu behandelnden Fall eines bloß 
einseitig begrenzten Leiters beschränken, da ja trotz der theoretisch unendlich großen Aus- 
breitungsgeschwindigkeit der Wärme praktisch die Temperaturverteilung in einer hinreichend 
kleinen Umgebung eines Endpunktes während einer kurzen Zeitspanne At nur von der für 
diesen Endpunkt vorgeschriebenen Randbedingung und den für seine Nachbarschaft an- 
verebenen Anfangswerten abhängen wird. Die Größe von It bleibt dabei natürlich bis zu 
einem gewissen Grade willkürlich, für seine Wahl wird jedoch besonders die für die Rechnung 
zeforderte Genauigkeit maßgebend sein, 

Ks seı also ",(r, 0) p(e) für x 0) und ’,(0, 8) f(t), wobeı noch p(V): fO)=V0, 
indem wir wieder mit #=0 die Schmelztemperatur bezeichnen. Ferner seien 9’ (0) und f’(0) 
vorhanden, wobei unseren Annahmen entsprechend g’(0) >O und f’(0)<0O sein sollen. Ist 
nun @=p+f mit noch unbekanntem p — die Tangentengleichung der Schmelzkurve 2 =s (f) 
im Punkte #—=7--0, dann handelt es sieh doch darum, für hinreichend kleines t je ein 
Integral der Gleichungen 

04; 0’ v7; 


i 3 i 5) 
da; _ (i 1.2) 
0 dr 


zu finden, die den folzenden Bedinzeungen zenügen: 


— 


a) a, (0,9) fit): I (pt) 0) 
EN a te 
h) I, (ar, V Wink ",(p i{.#) 0) 


7 
l 
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wobeı noch p so bestimmt werden muß, daß 


0%, 04, . 
lım (x, - z PER en ee 
I>\ (x VE /ı—= pt y 


Das den Bedingungen b) von (16) genügende Integral läßt sich sofort angeben, wenn 


PA 


man in (15°) und (15°) 4, = 0, q, = 1, a,’ t an Stelle von f und 9, (x) = (x) e?“” setzt, nämlich : 


F- 
ı PN , 
(2% pt) (x — pt —5)° (x—pt +5)? 
) 9, B2 | han? t u. | ! | 
(„pt 9,(e, 8): \o,(&)-1e a e e BE, „_. TB 
! a 3a,ynt r & | 
M) 
Für den Gradienten finden wir zunächst 
s 
D ji 
(2% pt) (x / -)2\ Es Et -+<)? 
j to? # D Ss / 
Ri) ), p g Ä N) e Fo | (£ | Ö | 4 as? Ö | t as? 1 } 1: 
= 0 0, r ( e)- 40 .. I\Ie . ade 
DE Ba: 2a, yat Jr rs oE&\ | 
0 
und daraus durch partielle Integration wegen y, (0) 0: 
s 
E_ (2x — pt) [ (x ‚Al 2 (x t 4-5)? 
07, p ini h | | vr 21 4 er | 
: 55'9,(x,)+ > -\o, (£)-je ’ un. j dE ..d(19). 
Öx au za,yrt J | | 
E m 
Wegen 9, (pt, b) =0 folgt daraus für z>p-t: 
N 
pet v 
04 — 
—? | EM! WE Gi... : GP - ;, ' 
UA ‘X pt d,] rt « 


0 
woraus, wie es Ja sein muß, wegen 9,'(0) = 1’ (0) folgt 


N, 
® ( ! “ ’ 
im —— vd). 
!>\0 X a pt 
Zur Bestimmung von %,(x,#) denken wir uns zunächst den Leiter ım festen Aggregat- 
zustand nach links ins Unendliche fortgesetzt und bezeichnen mit 9 (#) =0 seine Anfangs- 
temperatur, wobei g ())=0 und 9’ (0) vorhanden sein sollen. Dann ergibt sich entsprechend 


PA 


wie vorher, wenn noch 9, (#) = (#)  e?"'“ gesetzt wird: 


D a / . 
pP as“ un | (xo-pt 5)? ix pt- 1 
( FR 21 N 2t e: 
vl) > -\p(d)-je u WdE... . (8), 
za,yrt . | 
SF. 
woraus wieder wie vorher folgt: 
Nr) 
. 03 Br 
lim 7’ (0), 


I>VN UF x—pft 


und dieses 9’ (0) wollen wir nun bestimmen. Dazu bedienen wir uns am zweekmäßigsten 
einer geometrischen Ausdrucksweise, indem wir die gesuchte Integrallläche 9= 9, (x, t) in 
dem von der (— x)-Achse und der Geraden z=pt für t=0 gebildeten Winkel als regulär 
voraussetzen, so daß ihre durch 2=t=0 gehende Tangentialebene die Tangenten an die 
Kurven d=g (x) und d=f(t) und auch die Gerade z=pt enthalten muß. Ist also 

2.9 (0) +1: f (0) -9—=0 


f’(0) 
pP 


die Gleichung dieser Tangentialebene, so muß 9’ (0) sein, wenn sie auch die Gerade 


x» —=pt enthalten soll. Damit wird aber (17) 


(0) 
I Pr %s'0Q (0) /. OP, 


%; 
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und diese quadratische Gleichung hat wegen f’(0)<0 und 9’ (0) >O nur eine positive Wurzel 


p: 2 2,’ (0) + Y (x, 9’ (0))? t2ox,f' ()). 

womit wir die zur Ausführung des ersten Schrittes erforderliche Richtung der Tangente der 
Schmelzkurve kennen. 

In ähnlicher Weise ließen sich die Fälle behandeln, wo eine der anderen üblichen Rand- 
bedingungen vorgeschrieben ist, doch wollen wir darauf nieht eingehen, sondern untersuchen den 

3. Beginn der Schmelzkurve bei unstetiger Randbedingung. Wir haben schon früher be- 
merkt, daß ın diesem Fall, den wir gleich näher präzisieren werden und der in der Wirk- 
lichkeit sogar meistens die Regel bilden wird, unser Näherungsverfahren zur Konstruktion 
der Schmelzkurve beim ersten Schritt versagen muß. Da hier nieht der Ort ist, um auf alle 
mit diesem Fall verknüpften Schwierigkeiten einzugehen und Existenz- sowie Eindeutigkeits- 
[ragen zu erörtern ''), soll dermalen nur versucht werden, die klassischen Ansätze zu ver- 
allgemeinern. Es sei also etwa ein Nüssiger Leiter mit bekannter Anfangstemperatur gegeben, 
dessen Enden plötzlich unter den Schmelzpunkt abgekühlt werden, so daß also dort, „momen- 
tane* Erstarrung eintreten muß, wobei wir an der auch sonst üblichen Annahme festhalten, 
daß eine anfänglieb vorhandene Unstetigkeit der Temperatur nicht bestehenbleiben kann. 
Offenbar können wir uns aus denselben Gründen wie im vorhergehenden Abschnitt auf die 
Betrachtung eines einseitige begrenzten Leiters beschränken, dessen Ende sich bei @=0 be- 
finden und der sieh in der Richtung der (+ x), bzw. (— x)-Achse erstrecken möge. 

Betrachten wir zunächst den folgenden einfachen Fall, der sich mittels des von H. Weber '?) 
gemachten Ansatzes für das schon von J. Stefan’) in Angriff genommene „Eisproblem“ 
erledieen läßt. Es seı 


P.IED= U <R ur Bei 9,(2,0)=C,>0 für 2 >0 
» (s(d. N = I, ls(Hd.N=O 
| (20), 
2 ), | (° vr) ds 
H Ku AO» 
'\or/x—s are IP - di 





und wır setzen 


I (a, ld) = A; + Bi Po. ) CU 5 IR ‚ EEE EEE 2° 


wobei P (x) 


— 
u 
J 
. 
_ 
— 
r 


<. Die ersten zwei Bedingungen (20) liefern dann 


A,—-B=(C; 4,+B,=0 


und die dritte, wenn s (f) = 9-y t angenommen wird: 


A+B.:bl, )=0; 4,+B,-®|, 
woraus folet 


# i ) -) 
I; " p| A | ( 1,2). 
( ‚2 “«D 

2a; 
Schließlich ergibt die letzte Gl. (20) nach Multiplikation mit yrt die folgende transzendente 
Gleichung für 9: 


) 02 
u. ( pP tar? u: P tits? 
“2 Br ee 2 2 ] T N .).) 
A0tz ( a te A 
Ad, | L p| d | (dl, p| Ö “ 
“ö «DD e ()) 
2a, \>u, 


von der man leicht zeieen kann, daß sie nur eine einzige reelle Wurzel hat, wenn, wie an- 
genommen, ©, <O und €, >0. Während nämlich die rechte Seite mit öd wächst, nimmt die 
lınke für << von + nach © monoton ab, da man sich ohne Mühe überzeugen 

tt) Vgl. z. B. das Referat von G. Doetsceh: Les equations aux derivees partielles du type paraboliqne, 
l,’enseign. math. 35 (1936), S. 42 bis 87, bes. S. 4s ff. 


2), Riemann-Weber: Die part. Diff.-Gl d. math. Physik, 5. Aufl., Braunschweig (1912), 2. Bd., S. 118 ff. 
13) Über die Theorie der Eisbildung, Monatsh. f. M. u. Ph., 1 (1890), S. 1 bis 6. 
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kann, daß ihre Ableitung nach Ö nie positiv werden kann. Man braucht hierzu nur bemerken, 
daß für & = stets 


e—Xx? | ++ [1 P(x)] e—%° 2x |) ai \e .dE 2\(€ z\)e-:°’.dE& + 


Wenn. wie z. B. beim Wasser, ein erheblicher Unterschied in der Diehte der beiden 
Aggreratzustände besteht, dann hat man aus (22) zunächst das Vorzeichen von Ö zu ermitteln 
und zur Berechnung von Ö sodann o, oder o, für o zu nehmen, je nachdem $=Z0. Dazu ge: 
x.C, %0, 


nürt offenbar die Bestimmung des Vorzeichens von Verschwindet dieser Aus- 


a, d, 
druck, dann ist ö&=0, also x 0 und die Schmelzgrenze bleibt fest'°). 


Bei der nun folgenden Betrachtung des allgemeinen Falles, bei dem also die Anfangs- 
temperaturen des festen und des flüssigen Leiterteiles nicht konstant sein sollen, machen wir 
die durch das soeben erledigte Beispiel nahegelegte Annahme, daß man näherungsweise für die 
Schmelzkurve e=s(t)=6:y t setzen darf, solange f hinlänglich klein bleibt. Es wird dabei 
insbesonders auch die Frage zu prüfen sen, ob dann ö etwa auch von den an der Schmelz- 
erenze gebildeten Gradienten der Anfangstemperaturen abhängen wird oder nicht. Wir setzen 


also: 

,(r,0)=p,(r) für -<O und 9,(#,0)=9,(r) für 2 >00, 
wo aber nun 9, (0) ==, (0) und 9; (0) vorhanden sein sollen, und bestimmen zunächst je ein 
Integral der Gleichungen 


( ); 2 Pol ); 
Gi 


; 9 N p I “ 
of 92» 


so daß 
I;(#,0)=p; (a) für (-M@>O0 und 9;(Ö vti,)=V. 


Dafür machen wir den üblichen Ansatz 





| 2 | (x — s)? x 4 7] 
i ta;2t ‚ 4ay2t 
3; (a, 8) savzie ()|e ! f ! (ds 
{ 23), 
1)‘ x Ö yı (x-dVr)? 
r2a] 7 A (t 1)°!? .p taj2 dt ). dr 
( 
wobei 
Dr) . ee 
px) = P 0 je nachdem („1x0 


und die y;(r) die zum Bestehen der bekannten Diskontinuitätsformel '%) notwendigen Be- 
dingungen erfüllen mögen. Setzen wir der Kürze wegen noch 


es 


| 
vd) . \ ß (2a; +0) yt|—pl2a;y —d)- yYi |. E ı En 34° ° 


I. 
so erhalten die zur Bestimmung der 7; (r) dienenden Integralgleiehungen die folgende Form: 
! 
: de(Vt-—-Vr) 
] T ta;2 (Vi {Vr) 


r)*? u ; - (dt - y;(t) . . j j . (2). 


(— Did | 4 
a 7 4 Kr) (f 


1%) Die beim Kisproblem sich ergebende transzendente Gleichung für d unterseheidet sieh nieht wesentlieh von (22), 


aber da Weber von einer komplizierteren Rechnung sprieht, die zum Nachweis der Unität der Wurzel nötig sei, und 


Herr Fürth an der in Fußnote?) angeführten Stelle darauf auch nicht eingeht, so hielt ich diese Andeutung nicht 
fiir unangebracht. 


15) Vgl. das auf S. 192 ff. des in Fußnote?) a. ©. ausgeführte Beispiel, insbesonders die Gl. (18) auf S. 194. 
16) Vgl. z. B. Goursat: Cours d’analyse, »!PMesd., T. III (1923). S. 308. Die dort geforderten Voraussetzungen 


können für unseren Fall etwas gemildert werden. 
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12 Hube Uber das Fortschreiten der Schmelzgrenze in einem linearen Leiter 
Führt man hier statt 7 die neue Inteerationsvariable 
() ı yı d | / ’ 
t) 
2a vt—ı 3a; | VA, Yır 
ein. so wird 
0 f f 7 7 0 
do dr und - 
be; 7 (f T) f Q$ 
womit (25) übereeht ın 
/ A I)’ \ r Lıua,o d” La;0° 
YjıE) "1 Yıll-I ’x n'’e 
/ | T N ()” | (dl; Ö e)- > | dl,” 0° 
Da weıter nach (24) 
() 
y; (0) | I.4 (0). 25.) 
nn I 
so folgt, wenn Iım 4; (N#=7;(0), aus (25) 
ET i 
| I) ed)” bt; Ö 
(lo ( L)’ lg 


Um das auf der linken stehende Interral auszuwerten. 


zunächst 
| | 


*.do 


vr; (#) (25’). 
(24'), 

d\ I 
(O0) Pl, | (6) 

= 4; 

an 
setzen wir h und finden 

-(; 


23h (1-#& Ih le-®- 
-\ ee .d: \ „eds PD (h). 
| 7 | & } I. ı+E 
Setzen wır ferner 
p \ 
(ı (.r) \ . (dd & 
\1 & 
so wird 
| ] 
(ı (.r) (1 (.r) di =; 
: 5 = 
Derselben Differentialgleichung genügt aber auch 
” JT 
(tr (») i e | (Pi x )) 
Y T . 
und da G (0) ( (VO) = so Ist 
„T » 
ı f eN. | $D ] e-))-1 7 | 
| 
Differenziert man nun (25) nach F und ist 
hm |p # +; (dh eh 
> \ 
so folet aus (24) für t>0: 
.) ; ’ > IE 
| ee‘ I)’ . lı,o | . N) i 
Ne gara) ee" de=z pi )) 


‚ralG ia). fü 


17) Dieses Inte er: 


also dur: h 


S, 141), läßt sich „bekannte* Funktionen darstellen. 





(26’). 


r welches schon Raabe eine Reihenentwieklung angegeben hat (Journ. f. Math. 48 (1854) 
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woraus man ersehen kann, daß man für kleines f näherungsweise 4; (f) +4 (0) +e;y f setzen 


darf, wenn lım 9; («) 
x>+ 0 


—Yp; (V)=F0, 


Y’';: 


Wäre dagegen erst lim |[(# «) 


oo +-( 


) 


{ 


y; ()]=F0 und endlich für 


doch wollen wir von 


;>0, so würde dies die folgenden Überlegungen erheblich erschweren, 
diesem praktisch ja ohnehin bedeutungslosem Falle absehen. 

Nach dieser knappen Rechtfertigung des für die Schmelzkurve 
d-Y t haben wir nun die zur Bestimmung von d erforderliche Gleichung nach (1) her- 
zustellen, die trotz ihrer sehr komplizierten Form sieh mit der Gl. (22) identisch erweisen 


eemachten Ansatzes 


x 


\9. 
| | ' r OU... j 
wird, indem die z;(0) herausfallen werden. Zur Berechnung von lım \ „bilden wir zu- 
> ‚Vi + 0 ou 
nächst '°) 
4 
\ | ) ? (A \” \ - 
hd 0 |, at tat y 
\ > "I9pAK5)'ız|e e Id 
(U 5 d;] nt . O&| | 
SF 
! ”. Vi 
x 
| 14 ‚» 4ajtlt—ı) 
( [) |,;(0) kaj2t | + yi(r) r Ö] 7 | le dt; 
. . pP H\l z . rz (tr)! - dr : 
ayrz \yt lda’yr i—ı yi—ı | 
u 


0, für das zweite Integral 
Setzt man die so gefundenen Aus- 


wo man wieder leicht zeigen kann, daß für t>0, wenn 2>odyt 7 
reehts die Diskontinuitätsformel angewendet werden darf. 








u wi 7, 0%, : -. DE 
drücke für lım = und lım ® in (1) ein. so erhält man nach Multiplikation 
xm0Vt „08 v>dvVt+ro‘ L 
mit y ft für >00: 
te | 02 ) 
5 ta,? 
%, Ir. (0) Lay? e 1 () N) | 
€ =. IM. i [1 D 
a,\ya Kl) Van a, | 2a) | 
(28). 
N 
0)? 
I (ls 
Ka ‘ > (0) Ida? e u N) | N) | P ) 
erg y,(0)- 1+-P 10: 
a,\ya Ka) 20, | 2u,)) 2 


.)».) 


Da hier jedoch die 9; (0) gar nicht vorkommen, so muß diese Gleichung für 6 mit (22) identisch 
sein, wovon man sich auf Grund der Gleichungen (26) und (27) sofort überzeugen kann, wenn 
man noch 9, (0) = €, und 9,(0) = (, setzt, so daß also bei stetigen 9; (x) der Faktor $ nur 
von den 9;(0) allein abhängt und nicht auch von den 9; (0). 

Hat man also aus (285) bzw. (22) den Faktor $ berechnet, dann kann man nach (26) 
und (26°) die z;(0) bzw. &; bestimmen und damit nach (23) die nach einem hinreichend 
kleinen f vorhandenen Temperaturverteilungen im festen und flüssigen Teile Leiters, 
worauf man den folgenden Schritt nach der allgemeinen Näherungsmethode ausführen kann. 
Ks sei noch besonders darauf hingewiesen, daß die e; nur dann verschwinden würden, also 
in (23) die ;(r) nur dann in erster Näherung durch die Konstanten y; (0) ersetzt werden 
dürfen, wenn die 9; (0) =. 

Auf Grund dieser Ergebnisse kann man nun auch die Verallgemeinerung des Stefan- 
schen Eisproblemes erledigen, bei dem die Anfangstemperatur 9, (#,0)=y,(x) mit 9, (0) >O 
eines einseitig begrenzten flüssigen Leiters gegeben ist, dessen linkes bei #0 befindliches 
Ende plötzlich abgekühlt wird, so daß 9, (0, =f(t), wobei f(V) <O und lim |y t -f’(t)| vor- 

>( 
handen sein sollen. Wir führen diese Aufgabe auf den früher betrachteten Fall zurück, indem 
wir an das Einde unseres Leiters einen nach links unbegrenzten im festen Aggregatzustande 
befindlichen Leiter aus demselben Material anlegen, dessen uns unbekannte Anfangstemperatur 
),(2,0)=y, (x) sein soll. Setzen wir für genügend kleines # für die Schmelzkurve wieder 


des 


x==0-y t, so können wir aber jetzt nach (22) den Faktor 9 nieht sofort berechnen, da wir 


7, (0) noch nicht kennen, denn es wird im allgemeinen ja 9, (0) ==f(0) ausfallen. Es folgt 


2) . 
nämlich aus (23), wenn wir h z setzen, für t>0: 
10; 
! 
N » Ih? 
/ . 
—: r -dı 
(0, f) We 1 ‚ 
| -_ 39 A (f 4 
0 
18) Ve]. M. Gevrey: Sur les equations aux derivces partielles du type parabolique, Journ. de math. (VD, 


9 (1913), S: 320. 
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P} * [2 * T . ”.. 
oder. wenn wiır die neue Inteerationsvariable o einführen : 
1 
N » h? 
! 1 ] ) . do 
3.0, 8) .\ v.(t:0)-e 5 (29). 
l | 7 \ /ı (] ( ? 
il 
Daraus folgt aber, wenn lim z, (F-0) =, (0): 
>\ 
| 
h \e-:- 
r ’ > 1 ( dd ) * 
lım 9, (0,8) = y,(0)- \e Bu ” /() 


a0 ] JI : (] 


und wenn man die neue Integrationsvariable =: E = einführt, erhält man daraus nach 
kurzer Rechnung 
vd 1 —h-ya Hi -2lMEO) : . . 2 rer. 
Anderseits ergibt sich aus (26) mit Berücksichtigung von (27) nach einigen einfachen Um- 
formungen 
1,0): [I +-Pc(h))- [| —h-ya et. - Plh)]=y,(0)-P Ch), 
so daß man nach (29) für 9, (0) den folgenden Ausdruck erhält: 


Q | -P(h) (0 
( ) 7 0077, 
ni Pih) / 


Setzt man dies für 9, (0) an Stelle von ©, in (22) ein, so bekommt man zur Bestimmung von 
d die foleende transzendente Gleichung: 


J2 \2 
1laı? I its? 
H e : Fa e ” n JI ’ 
-f(0)- 9,10) * iu3 .d Fee 4; ° 
d, ) di, . R) 9 
Do, | er | 


2a, 
die genau mit der für das Stefansche Eisproblem sich ergebenden Gleichung übereinstimmt '”). 
Hat man 9 aus (22°) berechnet, so findet man z, (0) aus (29), aber wie wir schon oben be- 
merkt hatten, genügt es im allgemeinen nicht, in (23) die z7,(r) einfach durch z;(0) zu er- 
setzen, wir müssen also auch noch die &; berechnen. Da uns 9, (0) und auch d bekannt 
sind, ist &, durch (26°) bestimmt: um aber e, zu finden, müßten wir zuerst 9, (0) kennen. 
Deshalb bilden wir zunächst aus (29): 


N 9 J h . 
Or, / \ ’ ] ‚[ 6 Er ‘) 
-\y, (l-o).e ae ee 
of \ () | JT i Ai |; 2 
( 
Nimmt man nun wieder wie oben an, daß lim |y # »z,’(B|=e,, so folgt aus (30), wenn man 
t>\ j 
mit Y £ multipliziert, für >V0: 
N) J | u 312 
. OU, ! O a . ’ BT o3N’ 
lim /yt E, .\e y le see IV FU} 500 
>| ( f \ ) N JT i | Ö > \ 
u 
Ist also lim | £ -f(H 0, dann ist auch 8, = O0 und man kann für kleines # in (23) näherungs- 
'>( 


weise 7, (r) = 2, (0) setzen, während man sonst für 7, (r)=%,(0)+e,-Y r nehmen muß. 

Hat man nun e, aus (30°) bestimmt, so ergibt sich g,’ (0) aus (26), und damit können 
wir für 9, (SD) > 9,(Ö)+E-9,' (0) nach (23) für hinreichend kleines t und VDzExr=6:.y t die 
in dem während dieser kurzen Zeitspanne erstarrten Teile vorhandene Temperaturverteilung 
bereehnen und dann den nächsten Schritt nach der im Abschnitte 3 angegebenen Methode 
ausführen. 

19) (sl. (52) auf S. 22 des in Fußnote?) a. ©. 


20), Für das in (30) stehende Integral ließe sich sowie für das in (26) eine Darstellung dureh bekannte Funktionen 
leieht finden, aber da sieh unbequeme Ausdrücke ergeben, so sehen wir davon ab. 


STE En 
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Der Fall, wo ein flüssiger Leiter plötzlich an einen festen stößt, läßt sich noch dahin 
verallgemeinern, daß die beiden Leiter nicht aus dem gleichen Material bestehen und ins- 
besonders verschiedene Schmelzpunkte besitzen, doch wollen wir trotz des praktischen Inter- 
esses, das dieser Fall besitzt, dermalen uns nicht mit ihm beschäftigen. 


6. Ausgleich der durch die Tangentenapproximation der Schmelzkurve entstehenden Fehler. 
Es wurde bereits bemerkt, daß unser Näherungsverfahren das Bestreben besitzt, die infolge 
des Ersatzes der Schmelzkurve durch ihre Tangenten entstehenden Fehler von selber teilweise 
auszugleichen, und dies wollen wir nun nachweisen, wobei wir uns aus denselben Gründen 
wie in den beiden vorangehenden Abschnitten auf einen unendlichen Leiter beschränken 
können. Es befinde sich die Schmelzgrenze zur Zeit t=0 bei ==0 und es sei 
,(2,0)=y,(r)<O für -<O und (2,0) p,(X) >O für «>00, wobei die 9;(#) zweimal 
differenzierbar sein sollen. Die „wahre“ Schmelzkurve == ({f) 


(Bild I) sc ln ice die | arallı le zur x-Ac lıse im Abstande At 0 Im 2esft)ıste-ä)/ I 
Punkte ©, ihre Tangente in O dagegen dieselbe Parallele in einem 4, ya 
' EN 4 t est) / 5 na 
weiter rechts liegenden Punkt 0, den also unser Näherungs- | VAN „u 
verfahren als einen Punkt der angenäherten Schmelzkurve statt | fs" 37 
1 f N 





Tangente 0 9 durch ihre Sehne 0 Q approximiert, so wäre die Lage 
der Schmelzgrenze zur Zeit Jt natürlich richtig bestimmt, aber ' N > 
die nach der in Abschnitt 3 entwickelten Methode berechnete 
Temperaturverteilung würde von der „wahren“ doch ein wenig 


() liefern würde. Hätten wir die Schmelzkurve anstatt durch ihre | 
| 


abweichen, wenn statt des Bogens 00 die Sehne 0 0 genommen 49 

würde. Es seien dureh die Kurven %, und %, die Temperatur- | 3 . 
verteilungen dargestellt, die man etwa nach (15°) und (18’) ge- | 277 e 
funden hätte, wenn man das dort stehende p durch das p in der "777677475 ee 
Gleichung «=pt der Sehne 09 ersetzt hätte. Da uns aber p | / 72 
unbekannt ist, so setzen wir für das p in (18°) und (18°) das in | 3 7 

der Gleichung z=p-t der Tangente 0Q vorkommende p und er- / 2 

halten so die Temperaturverteilungen, die durch die Kurven %, Eos: 

und 9, dargestellt sind. Während wir aus 9, und 9, nach (1) mit 


sehr großer Genauigkeit die Tangente der „wahren“ Schmelzkurve 
bestimmen könnten”), erhalten wir aus 9, und d, nach (1) natürlich überhaupt keine Tangente 


der „wahren“ Schmelzkurve ==s({f), sondern eine Gerade: z==p-AIt+p*-(t It), wobei 
p* aus 
u, O7, ö . " 
974 — Mn A0o»D 
.” Ö K £ Ö L x p-J1 JA 


zu berechnen ist. Aber diese Gerade ist steiler als die Tangente der „wahren“ Schmelzkurve, 
so daß wir beim folgenden Schritte doch wieder näher an die „wahre* Schmelzkurve heran- 
rücken, ja sie unter Umständen sogar kreuzen können, worauf beim folgenden Schritte der 
Fehlerausgleich natürlich im entgegengesetzten Sinne stattfinden würde. Um dies einzusehen, 
genügt es offenbar, wenn wir zeigen, daß für ein hinreichend kleines £ mit wachsendem p 


— wenn also die Gerade 2=pt einen spitzeren Winkel mit der »-Achse bildet! der 
I», 0%, i 
Gradient abnimmt, “ dagegen zunimmt. 
Or Va or R 
PX 
Setzt man nämlich in Übereinstimmung mit Bild 1 4 ;(@)=e"'".g;,(#) für - bi >0 
» .. ‚ 
so wird zunächst nach (18°) 
ı 
r \ e 
p ja“? pt | (X pi 5)° IX pi 5° | 
ı (. f) - @,lE&- E taz? 2. ta? :-dE, 
da,yat J* 
TI. 
und man findet nach einfacher Rechnung 
ı) ’ 
= QO,l&) tr o,.lö)i:e ae: 7 
x ix—pt a,yrt 2a? Hi ni | 
SF. 


” 


21) Der nur geringfügige Unterschied rührt davon her, daß an Stelle des Bogens 0Q die Sehne 0Q als rechter 


bzw. linker Bereicehrand bei der Bereehnung’von %, bzw. %s benützt würde. 
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Damit wird nun 









gs A £ pi s)* (pt )? 
\: I. 23 \ 
( J), | ıYıle?, ta,2t / ( \. ta,2t | yp (2) BR z 
. E . > I ee Dr + ! = | en 
Oro Pa pi ad, ] at 2a; 0x 2a, a fi I 
F 
woraus durch partielle Integration wegen q, (0) =0 folgt 
u 
; (pt 
\2 
SE j \ (&)-e (ls 
O2OPx=pt 2a,’ ıt)”° ” 5 
JH 
(u 
J ! y (pi )2 
| ! ; ta; pP N /=- ta,zt | 
N (O)-.« 1 | er (& ( (&)e / 
1 ı»*) 2 2 > ı » 
dt, 2 „led, | 


Setzt man im ersten Integral der rechten Seite £&=-pt+2ay tn, so sieht man, daß es 
für #*>0 ebenfalls wie Y f verschwindet. Zieht man also den Faktor y ? aus allen Gliedern 
der reehten Seite heraus, dann werden die für 2>0 endlich bleibenden Terme negativ, also niımmt 
9 
(v4 


‚ „ für genügend kleines # mit wachsendem p ab. In ganz entsprechender Weise zeigt man, 


0 x 


\ ( 


BR, | ; ö 
daß | z für kleines £ mit wachsendem p zunimmt, womit aber das Bestehen der erwähnten 
( “ 


vorteilhaften Eigenschaft unseres Näherungsverfahrens nachgewiesen ist. Übrigens ist der 
soeben bewiesene Sachverhalt ja auch anschaulich sofort einleuchtend, indem ein zu schnelles 


\9 
. oR 
Fortschreiten der Schmelzgrenze nach rechts offenbar eine Verkleinerung von = und eine 
0x1 
22 Or, ’ 
Vererößerunge von \„ nach sich zieben muß. 
(X 


\ 


7. Ein Beispiel. Da sich das von Herrn Lachmann durchgerechnete Beispiel für Kupfer 
infolge eines dabei unterlaufenen Versehens’”’) zu einem Vergleich der von ihm entwickelten 
Methode mit der unsrigen leider nicht heranziehen läßt und ich mir auch die Wärme- und 
Temperaturleitzahl für Kupferschmelze nicht verschaffen konnte, so wollen wir unsere Methode 
auf einen /=10em langen aus geschmolzenem Zink bestehenden Leiter anwenden. Seine 
Anfangstemperatur sei 520° C und seine Enden mögen plötzlich auf 20°C abgekühlt werden 
Setzen wir die Schmelztemperatur des Zinkes gleich 420°C, so wird mit der in (22) benützten 
Bezeichnung 


rn 100: C,— 10. 
\ . ® x .. a )4°® * cal nn ® md f} ur . ® 
"erner sind die Schmelzwärme 4 = 26,6 er die Dichten o,=7,1l und 0,= 6,9, wofür wir im 
C O3 


Mittel o = 7,0 wählen, so daß 40 = 156 genommen werden kann. Für die spezifischen Wärmen 
setzen wir 


: cal cal 
2, = —, | Rt ae 
Y 4 
und für die Wärmeleitzahlen 
a eal a cal 
%, 0.249 3 BE a 0.15% 0» 
Chr Set i CH See 


womit sich für die Temperaturleitzahlen die folgenden Werte ergeben: 


5 ann em’ H P ‚em’ 
dt,“ 0.360 ’ u. — () [64 
C, oO, a hr B. O, = S0G 


Pd 


Die Endstücke der in unserem zwar für die Rechnung aber nicht im Wesen durch 
Symmetrie sehr vereinfachten Problem auftretenden Schmelzkurve wollen wir durch die 
Parabeln X=d-y T bzw. X==10- d:.y T approximieren. Wir nehmen dabei an, daß sich 
der Leiter nach der jeweils von dein betreffenden Ende ins Innere weisenden Richtung ins 
Unendliche erstreckt, so daß wir den folgenden Ansatz’) für die im festen und flüssigen 
Teile sich einstellenden Temperaturverteilungen am linken Ende machen können: 


Vgl. meine Zuscehrift an den Herausgeber der ZAMM 17 (1937), 8. 379 f. 
Vel. Gl. (40), (1) und (2) auf S. 221 des in Fußnote?) a. ©. 
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X Y 
) A+B:®D — 1: A A,+B,:-P — 
% Ba Fp® r) : chen APR r) 
’ n | h) | 
„.® er ME 
r\ Ö, 0, +. 
A: ER RB, - A. B. ” 
N) n h) 1 / 0 \ 
( «DD «)D) 
Pau) 3a.) 24.) 
Dabei bedeutet 8 die einzige positive Wurzel der Gleichung 
‚2 d2 
x, Ge u PT nz _,yı7 
— oO 7 FE Er — . ' 
dt, p| ) (Od. | N), - > 
()D) 2 “D 
3a) (&a,) 


für die man hinreichend genau d=0,7144 findet. Damit erhält man für die für kleines T 
vorhandenen 'Temperaturverteilungen die folgenden Ausdrücke: 


X 
(X<öd-yVT) 9(X,T) 100 + 666,56: D| 121 #) 


N 
0.S1 | I 





(X>d.-yT) 9,(X,T) 370,81 + 470,81: D| 


Um zu prüfen, bis zu welchem Werte von T wir (32) benützen dürfen, wenn der Fehler 
von d, bei X=5cm kleiner als 0,1°C sein soll, berechnen wir die Temperaturverteilung in 
einem Leiter mit a?=030, für den 9 (X,0)=0"C, dessen Länge | Idem und dessen Enden 
plötzlich auf 500°’C abgekühlt werden, aber ohne daß sich sein Aggreratzustand dabei 
ändert. Da also hier keine Schmelzwärme frei wird und obendrein a?’ >a,, so muß 

PR m: 
| 5 r)- 9, =, r\ I00. Nun ist aber bekanntlich 
EN 4 u. X+-2nl I X—2nl] 
J(X,T) 500. 12 KL nn )t+ PP 
| | \2ayt) | >ayıl | 


nn el 
„ SI. 


Damit also | =; | — 0,1 wird, muß T aus der folgenden Gleichung bestimmt werden: 


[fl „ | "u. 31 z„f 5 | 
dran da Mrz; + Pliaı 7 Rz, 7)" 22 01. 
oder hinreichend genau: 

0 


| derer 7 10°". 


v 


Daraus folgt, daß bis T 25 see die Gl. (32) sicher noch benützbar sind, da 


.) 
0 _ 


10 
D|5; 15) 


D (2.78) = 0.I9I99156. 


Wir haben somit insbesonders für 7, 2,25 see und 9 = 0-y T,= 1,0716 cm die folgenden 
Temperaturverteilungen: 
VO <X<10716: EIWSA<XN 10.0) 9, (= 00 + 666,56 + PD (055556 - X) .. 
(1,0716 = X <= 89284) 9, (X) = — 370,51 + 470,51 - P (0,52304 X) I en) 


Zur Ausführung des folgenden Sehrittes setzen wir zunächst für den erstarrten Teil 
10,9, =dyT, — nach (5°) die dimensionslosen Größen x und f an Stelle von X und T: 


4,.\° 

Ze, 3=1B- =: a FT: Fr -t=3,183 -t, 
: a, 

so daß 


I, (2) = p (x) 100 + 666,56: PD OSB) V! al). 


Mit dieser Funktion 9 (#) haben wir nun die folgende Randwertaufgabe nach (6°) zu lösen: 


a) da0)=pla) ar drei; b9l0,N wald +p/tN 0, 
wobei 
‚ ‘fs ‘fa Ö I /o\ car 
Azul 35 Passir"r 35.777 mar a 0.70885. 
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Kerner wird hier nach (9) a ; 1,4109 und 5=0. Wegen (0, f)==0 ist aber nun 
l’2 
die Formel (14) noch zu ergänzen, indem hier zum Faltungsintegral z,* @ (’) in (13°) auf der 


reehten Seite noch die Summe 


[’ j 


’ 


oe \ = -(r(t.)) I. \ / 


hınzuzufüeen ıst Setzt man ın 


statt 7 die neue Inteerationsvarıable 





On - at+t T In+x A r | / 
; und 2, | ; - int): 5 
> r (a ) 2 I (a f at 
so erhält man 
% Sir w7 
{ ’ \ 
I , I ( I’) 
/ ( d = 
| re ü 
und somit an Stelle von (14) 
4 senin 
J / 
SO) ’ n.(m 2X"): 
te ft 8 
np] Eu ’\ . \ & . € (de 
| T CT —n 0 
) J N 
| 
A 
y.? 
ww pr, . ( 7 x? ) ( 7 x? ] ) 
‘ \ WL It’ 
\ wie) ( E ds 
=} r/ | 
ı 7 I 


für die Temperaturverteilung im erstarrten Teile des Leiters, also für EX = 9,+ ps: Al, 


I ; °% « ur uf 
bzw. Or I, wo die Werte von f und F gegeben sind durch: 47 | -1+$ | .L ; 
A, dl, a—+1t 
Für die Temperaturverteillung in der Schmelze = yI, u = B—q, lindet man 
nach Einführung der dimensionslosen Variablen x, vermöge 
R m en Fo 0 2 d,4\° PFAOR 
\ 10 —2q4)c+q, .S5568 : x -+1.0716: 17 | | - t= 376,41 - 1 
dl, 
zunächst zur Zeit T, = 2,25 see: 
9. (.r) ı (.r) 370.S1 + 170.51: PD (6.1665 r + 0.585196) (Or 1). 
Damit ist nun die folgende Randwertaufgabe zu lösen: 
Aa) dar. O) g(a) (0 r I): h) "(p, #,#) Hl HP, 4, 0. 
. ( ’ I0- Ai .. » ’ ’ ) 
wobei wegen p, P=o/T und 9; = ‚> Pi zunächst folgt p, Ps 11,408 und 
‚3 i a, , 
| »).)*+) ) 4 » Er m 
nach (9) a =, ,70,04353 und = ,, und damit kann nach (14) sofort für dasselbe 17 
sp, 2 
wie vorher #9, (@’, #) berechnet werden, wobei q,+p,:ITEX ZI —- 1,—p,:AT, bzw. 
| |  (10—2gq,\2 0-2a.\ a-# 
0.’ 1 und fund # bestimmt sind dureh: AT! Ei 2 2 A). ‚, Weeen 
di, (, a+t 
2 ’ ’ . . . ’ | 
der Symmetrie unseres Problemes braueht man 9, (@’, f) wirklich nur für 0= 0" 5 aus- 





zuwerten. 
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= 2 u | ds. 
Nun hat man nach (1) den für den nächsten Schritt erforderlichen Wert von yp für 
( 
nr; 9 5 04, Or, 
I !,+AT zu bestimmen, wozu man \Y und \ kennen muß. Dazu beachten wir, daß 
Ö. 0. 
0%, 
04, AH, oO 2 0x Or 


ON or or ON (+ Ps° Wi 


und entsprechend, wenn nun 4, und p, durch g, und p, ersetzt werden: 


( ), 0 ), dx 0x dır' 
oX 9x dxz dA W—2q,—2p, AT 


wodureh einfach die Verlängerung des festen bzw. die Verkürzung des Nlüssigen Leiterstückes 
Ö ); . . . ‚ 

Jr MUS den Formeln, die sich aus (14’) bzw. (14) 
dureh Differentiation nach x’ bilden lassen, für #’=1 bzw. x’ —0 würde zwar im Vergleich 
zu der ohnehin zu leistenden Rechenarbeit keine bedeutende Mühe erfordern, doch findet 
man sie mit der hier angemessenen Genauigkeit bequemer durch numerische Differentiation 
der 9; (#’,f’) selber, wobei die relativen Fehler nie '/; vH erreichen, wenn man die 9; auf 
Zehntel genau berechnet hat. 


zum Ausdruck kommt. Die Berechnung der 








7 5 | | | } / 
40: ! / 
! 
oO 30 | | A 
& 
+ 
& ß 
20 = | 
Q 
| vr. | 
'Xcm En | | | |Xcm 
A 7 2 3 4 5 
Bild 2. Bild >. 


In der angedeuteten Weise kann ein einigermaßen geübter Rechner einen Schritt an 
einem Tage ausführen, wenn man es nicht auf eine besonders große Genauigkeit abgesehen 
hat. In Bild 2 ist ein Teil der von mir berechneten Temperaturverteilungen eingezeichnet. 


Nach etwa 33 see ist die Temperatur der Schmelze der Schmelztemperatur schon so nahe, 


daß ihre örtliche Anderung graphisch nicht mehr dargestellt werden kann. Da sich zeigt, 
07, 


daß - 
Ö x \x ) 


von dort an nahezu in geometrischer Progression abnimmt, so braucht man zur 


ds 


Bildung von IT 


das 9, (’) überhaupt nieht mehr berechnen, indem man sich bloß mit einer 


( 


7, 
0x 

Die nach dem geschilderten Verfahren punktweise berechnete „angenäherte“ Schmelz- 
kurve ist in Bild 3 (ausgezogen) dargestellt. Sie weicht erst gegen Ende des Erstarrungs- 


entsprechenden Korrektur begnügt, die unter einem '/» vH des von , , gelieferten Beitrages liegt. 


vorganges von der strichliert eingetragenen Parabel X =0,7144:y T merklich ab, und es ist 
sogar zu erkennen, daß sie die zwei Wendepunkte aufweist, die ıhr nach der genaueren 


‘)# 
y4 
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Berechnung des Herrn Laehmann zukommen sollten, da ja unser Beispiel einem der von 
ihm behandelten ziemlich nahe kommt °*). Man kann aus ihr entnehmen, daß nach 41.4 sec 
der eanze Leiter erstarrt ıst, und dieser Wert dürfte um höchstens eine Sekunde zu klein 
sein, wie eine erobe Aufsummierune der beı den einzelnen Sehritten zu erwartenden Fehler 
ereibt, wenn man annımmt, daß die „wahre“ Schmelzkurve etwas links von der eineezeichneten 
verläuft. Hätte man die Parabel X=9-y T als „Schmelzkurve* genommen, so ergäbe sich 
Fe Yo 
” u. s u de) ! u 
für die Erstarrungsdauer die Zeit 7 2 I9 sec. also um rund 20vH zu eroß. so daß man 
S2 
also, wenn man etwa für technische Zwecke nur eine obere Schranke für die Dauer des Er- 


starrunesvoreangees zu kennen wünscht. bloß die Gl. (227) aufzulösen hat. 


Wir wollen noch feststellen, wie sich in unserem Beispiel das Ö ändert, wenn man 
z»,— 2,0245 und ebenso a,’= a,’ —=0,36 annehmen würde. Die Gl. (22’) bekommt dann 
eine etwas einfachere Gestalt und ihre Wurzel wird 6 0,7169, womit sich für die Erstarrungs- 


dauer T= 148,6 sec ergibt. Daß sich hier d nur um 0,3 vH ändert, kommt davon, daß wir 
GC, IC, haben. Es ergibt sich auch tatsächlich eine bedeutend größere Änderung von 6, 
wenn (Ü, (,, wie man sieht, wenn man &, = — 100 und ,—=400 in unserem Beispiel 
setzen würde. Man erhält dann für #,==%, bzw. #z,=z, die folgenden Wurzeln von (22’): 


5 — 0204 bzw. d—= 0,1%, also eine Änderung von 7 vH. so daß man also besonders dann die 
thermischen Konstanten der Schmelze nicht durch die für den festen Zustand ersetzen darf. 
wenn (€, ein großes Vielfaches von (-- €,) ist, wie ja auch aus (22”) sofort ersichtlich ist. 
In diesem Falle sind also die in der Einleitung diesbezüglich geäußerten Bedenken wohl 
hbeeründet. 


8. Schlußbemerkungen. Es ist einleuchtend, daß sich unsere Methode mit geringfügigen 
Abänderungen der Formeln auch auf dıe homogene Kugel übertragen läßt, wie ich hoffentlich 
bald an der Berechnung der Dicke der festen Erdkruste zeigen werde. Auf dieses Problem. 
das sehon die Aufmerksamkeit von Clapevyron und Lame°°) erregt hatte, weist auch Herr 
M. Brillouin im zweiten seiner im Jahre 1929 im Institut Henri Poinecare gehaltenen 
Vorträge hin?®). Nach Erörterung der bisherigen Lösungsversuche in speziellen Fällen, wobei 
insbesonders der fehlende Existenz- und Unitätsbeweis für die Wurzel der unserer Gl]. (22) 
entsprechenden transzendenten Gleichung beanstandet und ein graphisches Verfahren zur Auf- 
findune eines Näherungswertes derselben angegeben wird, führt Herr Brillouin unser Pro- 
blem auf die Auflösung eines Systemes von Integralgleichungen zurück, aber „la complieation 
des noveaux.... nous laissent completement desarmes“. Die von ihm dazu benützte Grund- 
lösune der Wärmeleitungesgzleiehung ist allerdings eine andere, als wir sie im Abschnitt 1 
benutzt haben, aber die dabei sich ergebenden Integralgleichungen sind einer numerischen 
Behandlung genau so unzugänglich, wie die von uns angegebenen für den Fall des unbegrenzten 
lweiters. Bei der Ausdehnung seiner zunächst für die ebene Platte erhaltenen Ergebnisse 
auf die Kugel ist aber Herrn Brillouin eın Fehler unterlaufen, auf den ich hier hinweisen 
möchte. Die von ihm mit M, ın seiner Gl. (XXIII) bezeichnete Funktion genügt nämlich 
ear nicht der Wärmeleitungsgleichung für den linearen Leiter, und daher sind die in AXIV 
und XXV dargestellten Ausdrücke keine Integrale der mit XIX bezeichneten Differential- 
eleichung. 


Schließlieh möchte ich noch einige Ergebnisse der noch unveröffentlichten Wiener 
Dissertation des Herrn Robert Riecek aus dem Jahre 1924 erwähnen, auf die mieh Herr 
Hofrat Wirtingzer aufmerksam gemacht hat’). Es werden dort zunächst solche Integrale 
der Gleichung 


——_ 


.„ [du m ou\ ou 
Or: r or of’ 


die für m=0,1,2 die Ausbreitung der Wärme in einem Stab, einem Zylinder und in einer 


Kugel beherrscht, gesucht, die bloß von 2 abhängen, wobei r=9(f) die Gleichung 


) 
‚ Q (#) 
der „Schmelzkurve* sein soll. Aus der an der veränderlichen Schmelzgrenze vorgeschriebenen 


>24) Kurve % 0,8; % 0,1 der Abb. 7 anf S. 35? des in Fußnote?) a. ©. 


5) Memoire sur la solidifieation par refroidissement d’un globe liquide. Ann. d. Ch. et de Phys. 47 (1831), 


nn 0) hbıs Ib. 


6) Sur quelgues problemes non resolues de la Physique Mathematique elassiqne, Propagation de la fusion. 
\nn. de IV’Institut H. Poineare 1 (1931), S. 285 bis 308. 
7) „Die Wärmeleitungsgleichung bei Änderung des Aggregatzustandes. (Schmelzen u. Erstarren.) Lösungen 


Kugel und Zvrlinder.* Ein Teil dieser Arbeit soll demnächst in den Monatsh. f. Math. u. Phys. erscheinen, 
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Bedingung folgt dann zunächst, daß 9=y€C,-t+ €, sein muß, wobei €, und €, willkürliche 
Konstanten bedeuten. Dann muß aber «(z) der folgenden gewöhnlichen Diflerentialgleichung 


veenügen: 
d’u m Ü, du 0 
de 1% u 1 
woraus sich für C,=1 und (,=a ergibt 
/ 
y. (dl /t + 
r dr 
u(r,)=B,+B,:\e Es „Mm 


Von den noch verfügbaren Konstanten BD, B, und «a ıst B, durch die Temperatur an der 
Schmelzgrenze bestimmt, während DB, und a willkürlich bleiben. Die auf diese Weise ge- 
wonnenen neuen partikulären Integrale lassen sich für m =1,2 durch den Integrallogarithmus 
bzw. das Gaußsche Fehlerintegral darstellen. Natürlich sind sowohl bei diesen sowie auch 
bei einigen anderen noch aufgestellten neuen partikulären Integralen die Rand- und Anfangs- 
bedingungen beinahe schon vollständig von vornherein festgelegt. S64 


A Non-linear Wing Theory and its Application to 
Rectangular Wings of Small Aspect Ratio 


by William Bollay'), Cambridge, Mass. U.S.A. 


Zusammenfassung. Die vorliegende Arbeit befaßt sich mit der Berechnung der Kräfte, 
die auf einen rechteckigen Tragflügel mit kleinem Seitenverhältnis wirken. Erforderlich dafür 
ist eine nichtlineare Tragflügeltheorie, die nicht die übliche Annahme macht, daß z. B. der 
Anstellwinkel klein ist und daß die abgehende Wirbelbahn in der Tragflügelebene liegt. Der 
wesentliche Unterschied der folgenden Entwicklung gegenüber den von Prandtl und Blenk 
eegebenen Theorien der tragenden Linie, bzw. der tragenden Fläche, besteht ın der Annahme, 
daß die abzehenden Wirbel die Tragfläche unter einem gewissen Winkel zu ihr dort verlassen, 
wo sie gebildet werden. Deshalb erzeugen die abgehenden Wirbel eine Geschwindigkeit 
tangentiell zur Tragfläche und dies ändert die resultierende Kraft an der Tragfläche. Da eine 
eleichförmige Auftriebsverteilung angenommen ist, liegen die Wirbel in den Endflächen des 
Flügels. Die Resultate zeigen gute Übereinstimmung mit den Versuchen von Winter. Im 
Grenzfalle des verschwindenden Seitenverhältnisses gibt die Theorie einen Normalkraftbeiwert 
U’ =2sin’a, während im Falle eines unendlich großen Seitenteiles sich das bekannte Resultat 


UN =27zsina-cosa ergibt. 


Introduction. The theory for flow about a wing of large aspeect ratio has been worked 
out quite completely along the lines suggested by Prandtl. This ıs the flow which is of 
prineipal interest for the main lifting surfaces of present day airplanes. However, for some 
unconventional airplane wings and especially for the control surfaces of dirigibles, the lifting 
surfaces have generally a very small aspeet ratio and the existing theories have so far been 
inadequate to give either the magnitude or the distribution of the forces correctly. 

The Prandtl theory which replaces the wing by a lifting vortex line holds very elosely 
for large aspect ratios k of the order k>4. For smaller aspeet ratios, Birnbaum [1] and 
Blenk[2] carried out an extension of the Prandtl theory by replacıng the wing by a 
distribution of vortex lines, thus giving the so-called hfting surface theory. This theory holds 
approximately down to aspeet ratios of the order 1. However, here the experimental lift 
curves indicate a eurvature concavely upward, and Blenk stated already ın 1925 that some 
additional effeet causing this eurvature becomes important at small aspeet ratios. He realized 
that this effeet could not be obtained by means of the existing wing theory, because this 
latter is essentially a linear theory. Since this time no theoretical work has been done on 
the wing of small aspeet ratio?). In 1955, however, Winter [3] reported on a very extensive: 


1) This paper is essentially an abstract from the author’s thesis-work earried out at the California Institute of 


Technology, under the supervision of Prof. Dr. v. Kärmän. lLacknowledge gratefully on this oceasion Dr. v. Kärmän’s 
« > . 


suggestion of the subjeet of my thesis and the very effective help during the work. 

2) Since the writing of this paper W. Kinner has treated the case of the eireular wing in the Ingenieur- 
Archiv, 1937 (p. 47), and F. Weinig has discussed the application of the lattice theory to small aspeet ratio wings in 
Luftfahrtforschung, 1936 (p. 405). 














.. a. ' ' Rn . 4. angew. Math. Mech. 
3) Ola \ E \pplieation to Reetaneular Wines ol Small Aspect Ratio Bd 19 Nr | Febr 1020 





experimental investigation carried out at Danzig on wings of small aspect ratio. Also at 
about the same -time, interest arose in America in the forces and pressure distribution on 
dirieible fins in conneetion with the ıinvestigation of the Macon aceident. Dr. von Kärmän 
then suggested to the author that he carry out an investigation for the wing of zero aspect 
ratio, assuming trailing vortices leaving at some angle to the wing, and to try to extend these 
results to the wing with finite aspeet ratio. The present paper represents a report on the 
results of this investigation. It represents a first step in the establishment of a non-linear 
wing theory. Since such a theory ıs principally useful for very small aspeect ratios (smaller 
than about one), the approximations which were made to simplify the ealeulations were made 
so as to fit especially for this region of aspect ratios. As a result, the theory which gives 
the eorreet result for both the limiting cases of zero and infinite aspeet ratios holds best for 
aspeet ratıos below about one. 


In the present paper, we shall restriet ourselves to finding the normal force on a 
reetangular flat plate wing. In addition, in the appendix we shall further investigate the 
limiting case of the wing with zero aspeet ratio, that ıs finite span but infinite chord. Work 
is also under way to apply the same theory to the dynamie lift of dirigibles, which exhibit 
a behavior very similar to wings of extremely small aspect ratio; however, these calceulations 
are not vet completed. 


A. Theory for Rectangular Wings of Small Aspect Ratio. 


1. Assumptions. Let us assume that a flat plate wing is inelined at the angle of attack 
a with respect to the free stream velocity U. We can resolve this veloeity U into a component 
alone the wing, U-cosa, and a component normal to the wing, U-sina. Now it can be 
shown that for purposes of a mathematical analysıs, the solid boundary of the wing can be 
replaced by an equivalent system of bound vortices. The strength of these vortices is determined 
by the boundary condition that there be no flow through the plate. This means that the 
normal component of the free stream velocity is cancelled by the normal component of the 
induced velocities at every point on the plate. A theoretical analysis of the wing of zero 
aspeet ratio shows that the vortex strength and the induced velocity are constant along the 
span. Experiments show that this ıs still a fairly good approximation at an aspect ratio of 
about 1 (see Fig. 14). Thus for the analysıs of wings of small aspeet ratio, we shall simplify 
our analysis by the assumptions that 


a) the vortex strength is constant across the span, 


b) the downwash along the span ıs eonstant and equal to the value at the center of 
the span. 


let us now consider the vortex distribution corresponding to a rectangular wing of 
small aspeet ratio. We shall denote the variables along the chord and span by x and y 
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respeetively, and eurrent variables in the same direetions by Sand n. The eoordinate normal 
to the plate is denoted by z, the positive direetion being chosen upward. The origin O0 is 
ehosen in the center of the wing. Then an equivalent vortex system replacing the wing 
consists of a distribution of bound vortices which is assumed constant along the span and of 
variable intensitv y(£) along the chord. At the ends of the span each bound vortex is as- 
sumed to leave as a straight trailing vortex at some angle © to the wing. These trailing 
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vortices are assumed to lie in the two vertical planes passed through the sides ofthe wing. 
Thus we get the system of horse-shoe vortices shown in Figure 1. 


Our first problem ıs to caleulate the induced veloeities along the center line of the plate 
for any point « along the chord due to the bound and trailing vortices respeetively. "These 
are obtained by the Biot-Savart Law. 


2. Induced Veloeity Due to Bound Vortices. A single bound vortex of strength y gives 
an induced velocity which is normal at every point to the plane of the wing and which has 
a magnitude 


l 
.), 
/ J) - () 1 Ä _ 
11° . :- (COS I), + cos (7,) 6Y I, Er 
In h 1 . ’ Ir a > l,\2 
| .) | oe Ss) 
We have to integrate the effeet of such vortices of variable strength y(£) from < to 


.) 


Sri 


— +, and obtain as the induced velocity normal to the wing due to the bound vortices 


In the plane of the vortices the tangential induced veloeity ww, is zero. However, an infinitesimal 


| (a) EEE | v (.#) 
distance above the plate wx, -",°, and an infinitesimal distance below urx, , 
’ x & .b2 h | 

Introdueing dimensionless eoordınates 5 > 35 >58, 40 > k:=aspect ratio, we thus get for 
the induced veloeities at the plate 

= tag 

| k-{ 7 ds 
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3. Induced Velocity Due to Trailing Vortiees. A pair of trailing vortiees of strength 
give an induced velocity mw, normal to their plane of magnitude 


v1 
m" ,, 20;-sıind 2). —— -(cosi+1)- sın d 
I h 
where 
, b \? Sn. wi 
h“ | > r (2 — &)”-sın? 9 
8 hi 
Sin « 
hı 
R Pr =)- COS (-) 6 &) cos () 
cos 4ı 


\ (x E)”.cos?9) + h? | | h\? 
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As the integral of the downwash due to the trailing vortices 
coming from both ends of the span, we get 
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A The normal component w;,, and the tangential component 
x,of the induced velocity are respeetively 


























W>. I, COS (-) 





and 


rs = mi sind WUx, Un sın O: Wr,‘ tan. 


ig. 3 Introdueine dimensionless coordinates we can therefore write 
| FE EEE EFERSE 
| nei v(£E)-dE (.r E)-cos I) | 
>, r \9 9 | Fr 
27 )R®+(2— £&)-sinO |yk’+ (ec — £)® Mi 
| 
Wr ;,. tan 7) 








4. Derivation and Solution of Integral Equation. Expressing our boundary condition now 
that the normal eomponent ofthe induced velocity at the plate is equal and opposite to the 
normal component of the free stream velocity, we get an integral equation for the determination 


(X) Tr N: (X) !/’.sın a 
or 

"U v(& d= k f v(E)-cos@Y).d:E 
: 2? + (x — &)?- sin? © 


I; \ v(E)-cos’(@): (x )-dE 
27) |k’+(& — 8° - sin? 0]: yk?+ (ec — 8)? 


In order to get the pressure distribution over the wing, it is necessary to solve this 
integral equation for y(£). In the present investigation, we shall restriet ourselves to finding 
the total force alone. "This makes it necessary to satisfy the integral equation only in the 
mean. We assume some reasonable distribution for the vortieity over the chord such as 








- 


y(£) - (whieh holds for large aspeet ratios), and then use the integral equation 


0 
‘ I 


2 - 
[2+8 


to determine the constant y,. 


KExpressing the fact that the mean value ofthe induced velocity over the chord is equal 
to the normal component of the free-stream velocity, we get 
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or Iinterehaneing the order of integration of x and & 
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where 





1 (£ da | k+yk+(l+5” 1-58 
(6) \( Dykr(z— 5 % Pk +yR®+i—H: 1+E 
l 
dr | (1+-5)-sın (1 &)-sın | 
# (S) \x: + (x &)?. sin? (-) k.sın® E vos l; | tan ; u 
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(‚x E)- dr 
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Now substituting these values in the above equation and letting 
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and solving for y,: 
Ir-U/.sına 
de F\, (k) -+- Got () e F, (de, ()) ot () P F, (k. ()) 


The problem is thus reduced to solving the three integrals FF, (k), F,(k, 9), and F,(k,o). 


5. Discussion of FF, (I). 


The second part of the above integral is integrable by elementary methods. However, 


the first cannot be integrated by any of the elementary methods nor does ıt adapt itself to 
solution by contour integration or by a differential equation. Some method of expansion was 
thus the only alternative. The first attempt was to expand the radieal in the domains I + &)<k 
and (1 +&)>k and then expand the logarıthm and integrate. This was carried out. However, 
the result was very long and the eonvergence of the series rather doubtful. A better method 
of expansion was found to be the following: As a first approximation, since k <1, we can 
write (letting wu =1+%£) 


log(k+yk’+ wu)=log(k + u). 
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This approximation is good if k<u over most 


ar a of the domain, but it does not hold near u —0®, 
/u) } ° ” . .. . 

u .e. £ | (the leading edge). Since the weieht 

ee [unetion becomes infinite in that vicinity, it is 


15 > Nhecessary to have a good approximation there. 
A better approximation can be obtained by letting 


exact: f{u)-JoglkrVk“u“) 


InnrY Jagd Fru)*r ge 
pror: f{u)=kglKru)’ae io 


log(k+ yk’+w)=log(k+u)+a:e 





where constants a and / are determined to fit 


Fir. 4. the difference eurve of 


log(k+yk’+ u?) -log(k +u)jatu=0 and u=k. 


Ihe exact and approximate formulae are plotted up below for the case k —=.2 and show a 
very close fit indeed (Fig. 4). It was found possible to integrate the new funetions. The 
lirst came out in terms of elementary funetions, the second in terms ofthe Bessel funetion 
with imaginary argument ofthe order of 1,1,(4). The result of this rather long analysis is: 








kr ) 
Ik I I; E+23'” 1.08 
j’ 1.435 e h | 
F(d= gt: re .+435 Au |. 
Kyar! 
6. Discussion of F,(k, ©}: 
| a/I-EI (1-+&): sin ©] 1-9-sinO]\ |. 
”..(@)) = A | —-t; . (dl E 
F,(k, \| 4 z ran 1: | tan | 1: e. 


-] 
This integral can be evaluated exactly using the method of contour integration. The result, 


not derived here is 

















| | f ‚je sın 2 
sınl „tan! 
| (2: sın O 2 k | 
Fk) =2xrx-tan tı -tan"’i: E | 
’ I 12sin@Q)\: | ‚|2sin 
I+| 1; t cos| > tan l; 
ai, | 2.-sınQ 
his expression can be reduced so as to involve only one parameter, namely x 1. "hen 
‚f{1 
| sin | 5 tan ) | 
| FF, (k,)=2n-tan”! u—4n-tan”! — Ä 
j l-+ u -cos| „tan u) 
7. Diseussion of FF, (k, ©): 
| S tan Oo i Zu tan oO n „N 
Fk, | „tan + (1-+ 5] —tan:'| —+(1- E)7|,-dE. 


This integral, like the first one, cannot be exactly ealeulated by any of the conventional 
methods. However, ıt ıs again adaptable to the same method of expansion as F, (k), by writing 


‚tan o 1: _ A RR tan \ 
an "+ u?) =tan ur +a:e* 
ie |" 
which (1 + &)= u, the eonstants a and 4 being again adjusted to fit the eurves at u —=0 
The exact and approximate formulae again showed very elose agreement when 
m ’ ‚ftanY-u i 
plotted for the case k == .2. The term involving tan a u could again be evaluated bv 


In 
and u k. 


the method of contour integration similarly to F,(k, ©. The term involving a-e=-*" again 


came out in terms of the Bessel funetion /,(4). The result of the analysis, again not 


dderived here, 18 
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| ‚tan 2tanQY . [1 ‚2tane 
Eu | n | cos (5 tan I; + J; sin | > tan 1; 
(k, 9 I 
| 2 .tan | Gun dach 1 
"7 
H2n-0-.e-*.1,(4) 
where 
ti 2 tan ©) 
/ log | u | | ıı 
I: \ (-) | 
Fun e ee ztanY 
he first half of this ean be reduced in terms of the parameter » Be hen 
| ee 
| cos | > tan 'r) sın | > tan '» 
F.tk,)=4n - z - +-2n70-.e-*.1,(4) 
v ” Pi-rr vl-+»’ 
Now 


Irl/.sına 
» F,(k)+ cot 9: F,(k, O9) — cot 9: F,(k, 9)’ 


Since F,(k, ©) can be expressed as the sum of two terms which are funetions of the simpler 
variables » and 2 respeetively, it is convenient to define new parameters B(r) and C (A), in 
which we can at the same time absorb the faetor 47. Thus we write 


u F'(k) F,(k, ©) | Rn F,(k. ©) 
F\ (k) Ada) . Bier) +9. CA) h 
ia | kn I 
where 
| i | 
COS > tan ") sın | > tan” | 
’(r) " 
) ) | y? | y? | 


| 
Ci=—ze-'L, (4) 


then 
!!’.sın a 


© F (k)+ecot 9|A(ı) + Bir) 9-C(A]' 


1 
8. Caleulation of Induced Velocity Tangential to Plate. In the ordinary wing theory 
which assumes small angles of attack, the effeet of the induced veloecity tangential to the 
plate is neglected. Since the angles of attack"whieh wings of small aspeet ratio reach before 
stalling may run up to 45°, it is necessary to consider this velocity in the ealeulation of the 
. ryv . . . ö . j v(x) 
forces acting. The bound vortices induce a veloeity tangential to the plate which is +, 


v (&) 


-) 


_ 


at the top of the wing and at the bottom. Their mean value at the wing is therefore 


zero. The trailing vortices, however, induce a veloeity which has a component along the 
wing equal to w, == mw; -sin$ = w.-tan®©. The mean value of w, over the chord is then 


eiven as 
m, | 5 IF, (k,9) - F,(k,O)=y,|A(ı) + Bir) -9-C(A)] 


and since 


v[Ale)+ Bir) 9. C(A)|= U-sina-tanY - Fi’ (k)-tan 9» y, 


| W. U ’ sın A: tan (-) R (k) . tan (-) ® Yo 
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9. Determination of Inclination of Trailing Vortices. The inclination of the trailing vortices 
with respeet to the plate which we have denoted by © is determined by the Helmholtz 
vortex law at every point in space. This law states that a free vortex follows the fluid 
partieles. 'Thus the vortex will lie along the resultant velocity vector. The inelination of 
this resultant velocity veetor has been assumed constant all during this development. If it 
should be variable, we should therefore desire an effective mean value for the inclination. 
However, this is quite diffieult to obtain exactly, and as an approximation we shall therefore 
choose, for the effective inclination ©, the angle at which the vortices leave the plate. The 
justifiecation for doing this is that the vortices right near the plate are the most effective in 
determining the downwash there’). 

With these approximations, the determination of the angle of the resultant velocity 
veetor at the side of the plate ıs a simple matter in terms of our previous results. For it 
is clear that the induced veloecities at the sides of the plate are exactly half as much as the 
induced veloeities would be at the center of a wing of twice the span and thus twice the 
aspeet ratio. This is due to the fact that ın the plane of the side of the plate only the 
opposite vortex sheet can induce any velocities in the x and z directions. We obtain for the 
mean values of the induced velocities: 


Y. Filz) F(2h) 
Iz, In \ Yo’ ) 
0. cos [F,(2k,O9)— F,(2k,O)|  Yy.-c0tO | (u\ „[r A| 
"2 4.sinQ| 2 | 2 \ala)+2lo) (5) 
I t/2 
| 1(y@& 
x, 2 n \ 2 dg 8/0 
m 


Wx,=—Mm;, tan Co 


and the components of the free-stream velocity are U-sina and ÜU.cosa respectively. The 
inclination of the resultant velocity is thus 
{ () U u sın A >, Wzs 
an € — _ 
l[.cosa—+ wx, + W x, 














® ’o | u] ' u v 4 / \ 
sın «a 1 .—+IF, >k) + cot 9.14(5)+B(7,) () . ( (> I} 
| tan © RR = — = 
| Yo | T | hh > ’ Y / u 
| COS ar I > | . Als)+B| 5) (). ( 7 h 
De I'm 
Since 7 ıs a funetion of © this equation for © Al 
nce etlon oO 1 ’ Ss [U2? 
ean only he solved by iteration. As a first S/ 





a & 
approximation, we can assume =, calculate gr 
-_ [A863.5] 
1A9909.>) Oros u 
0 > . R > m 
T and the funetions A, B and € and then find SE 2, 
uw 
a new expression for tan ©. A second appro- Pig: 9. 


ximation will, with a little experience, give very good agreement between assumed and 
resultant value for ©. 


10. Caleulation of Normal Force. By applying the formula for the force on a vortex, 
we can now ealeulate the normal force N acting on the wing as 


N=o(U-cosa-+tm,)- Tb 


where 


- 8 
- 





—\y(&)-d: "\| is e 


') Another approximation which we shall use is that we assume the trailing vortices to lie in the two vertical 
planes passed through the sides of the wing. This amounts to negleeting the 4-component of the induced velocity 
when applying the Helmholtz condition. 
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N=o[U-.cosa+ÜU-sina-tanQ -U.F/’(k)-tan 9.5 \- 17. 7 z hi 
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3 de U?.bt 
a rs x! 
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This formula for Cy represents the final result of our analysis. We shall determine its 
asymptotie behavior at aspeet ratio k=0 and k=x. Then we shall compare the theoretical 
curves of Üy vs. a with some experiments made by Winter on wings of small aspeet ratio. 

For convenience, let us summarize the important formulas which are required for the 
ealeulation of the normal force eoefficients. We have 


nr | no ] 
‚vo 
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43 | 
sin ( „tan' “) 
Au)=z-tan”! u—tan-' — (Fig. 7a, b) 
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The last five funetions have been plotted ın terms of the parameters k, O,/, u,» in the figures 
dlesignated above. 'T'he parameters « and » have been defined as 
2.sın O 2 tan O 
TE 1; and Pa 

11. First Limiting Case — Infinite Aspect Ratio (k=x). If the approximations which 
have been made in evaluating the integrals are still valıd, we should get for this case 
the wellknown formula given by the Prandtl theory, namely C,=2rzsina and thus 
UN =2nzsina-cosa. We have for this case 


k=o F'’(k) = .530 y 
ı WM) 2 —9.sina 
Ü 
uU Aluı)=V0 
r ö } ER R tan ©) \ 
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„—— Bir) —0 \ > 








u Ci)=Pd or|Uy 2a: sina-cosa . 
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Ihus we find that the normal force coeffiecient for the wing of infinite aspeet ratio is given 
correetly. It is found to be independent of the angle 9. This suggests that probably for 
Iinite but large aspect ratios the angle of the vortices © is not very important and that the 
’randtl theory which assumes 9=a should therefore be a gooll approximation. 


12. Second Limiting Case — Zero Aspect Ratio (" 0). This case of zero aspeet ratio 


corresponds to a wing of finite span but infinite chord. The leading edee is at r—= x 
the traılıngz edee at -, For thıs case our formulas reduce to 
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This last equation ean be solved for tan © and we obtain 


cot a | A 
tan 9 5 ty Veol’a-+2 

and for the normal force eoeffhieient 

Bu | 9 .| e si Q | sın a | 

IT * sın a» tan (9 -lcos a + sın a - tan € sın a-tan (9 + | | 

\ Rn |2tan 9) 
or 

Un > sın? dA | 


- 
This simple result is very interesting, for it indieates that the wing of very small 

aspeet ratio behaves very much like the elassieal Newtonian flat plate which experiences a 

normal force proportional to the time rate of change of momentum of fluid partieles hitting 

against it. The normal component of the veloeity of these fluid partieles would be V U 

and thus if we assume inelastie impaet, we get for the force N 

| 


-9sın’ a- > oU’:.8S 


-SINn GA, 


N= o V:-V.sS 
C\ >2sın? a. 
Experiments made by Winter have also shown that a sin? law holds for wings of very 
small aspeet ratio. Flügel [5] extrapolated the experimental values to zero aspeet ratio and 
stated that the eoefficient of sin’ a1s2. Thus the present theory seems to check experimental 
values also for the second limiting case of k—=V0. 
It is elearly seen that the result in this case depends very much upon the angle ©. 
This suggests, therefore, that for wings of small aspeet ratio the assumption that the vortices 
leave at some angle to the wing is of basie importance. The result is that a nonlinear 
theory such as the present one Is required to get at the fundamental behavior of such wings. 
The results for the case of k=0 can also be derived from very simple direet con- 
siderations. These are given in the appendix. 


13. Results and Comparison with Experiment: We have already seen that the present 
theory gives the correet result for the normal force eoeffieient in the lrmitinz cases of zero 
and of infinite aspeet ratio. For intermediate aspeet ratios, the agreement of the present 
theory with the experimental results of Winter [3] and Zimmerman [6] is shown in 
igs. Il and 12. "The eurves obtained by the Iifting-line theory are also plotted for comparison, 
though it is, of course, understood that this theory was never expected to be applied to the 
very small aspeet ratio wings. Uonversely, the results obtained using the present theory are 
shown for aspeet ratios as large as k=6, even though it is elear that some of the basie 
assumptions do not hold any longer. 
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or the very small aspeet ratios of k= 1/30 and k— .134, the theory gives the shape 
of the Cy vs. a eurve very nicely. The theoretical values agree within about 5° with the 
experimental ones. For the larger aspeet ratios of k—.20, .35 and .50, the present theory 
gives at least the order of magnitude, the maximum diserepaneies with experiment being of 
the order of 20°. This is probably due to the approximation used for the angle © at 


which the trailing vortices leave. ‘As Fig. 13 shows, there is a peculiar increase of the value 
oT © as a function of k, with a maximum value somewhere near k—=.35. Perhaps a better 
method of determination of the mean value for © will improve the agreement with experiment. 
It is to be noted that even for these very small aspeet ratios, the initial tangent already 
begins to agree with the lifting-line theory. This tendeney increases with increasing aspect 
ratıos kb ,66, I, and up. At k=,66 the theoretical ceurve begins to be too high rather than 
too low as it was for the smaller aspeet ratios. At k=1, it is uniformly about 10°Jo too 
large, while the lifting-line theory gives the initial tangent quite elosely but then departs and 
eives results as much as 20% too low. For k =1.25 and 1.5, the experimental points lie 
somewhere between the present theory and the lifting-line theory results. At k=2, the 
Iılting theory fits already pretty well, and of course very well at k=6. For these larger 
aspeet ratios, the present theory gives values about 25° too high. It would not be expected 
to hold, since the assumption of a constant lift distribution across the span breaks down. 
Also the approximation used in expanding the integrals F,(k) and F,(k, ©) become less 
satısfactory at large aspeet ratios. 

It is interesting to consider the validity of the fundamental assumptions used in the 
present analysis. The assumption of a constant lift and constant downwash across the span 
is verified exactly for the case k=0 in the appendix. However, pressure measurements made 
by Winter (Fig. 14) show that even at k=1 the assumption of a constant lift across the 
span is still a good approximation. Thus it should be good at least in the interval O<k<1. 
2 — x 
2 + 


The assumption for the vortieity distribution along the chord y=y, \ also appears 
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reasonable from Winter’s measurements at k—=1. The actual pressure distributions showed 
very great suction peaks at the leading edge as postulated by the theory. In addition, they 
showed, however, smaller peaks near the sides and at the trailing edge, which have not been 
included in the theory. The measured pressure distributions indicated a considerable variation 
with the angle of attack. This has not been taken into account in our theory. It indicates 
that the true solution of the integral equation for y(£) is probably quite complicated. For this 
.) » 
reason, the assumption of some mean distribution funetion such as y: Yo| rs is probably 
as good an approach as can be made theoretically. Rn: 
The assumption that for small aspeet ratios the vortices leave at some angle to the 
plate seems also to be verified by flow pietures taken by Winter. He finds that for very 
small angles of attack, a smooth Prandtl theory type flow takes place. This checks the 
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fact that the lifting-line theory gives the initial tangent of the Ey vs. a eurves. At larger 
angles, however, a strong tip flow oceurs with vortices leaving the plate at the sides. These 
vortices actually roll up immediately into a single vortex core, and according to Winter 
leave at approximately the angle of attack. The two vortex sheets which we assumed in 
our theory are apparently unstable; their effects are, however, approximately equivalent to 
the single vortex core, as Is seen by the agreement between theory and experiment. 


14. Conclusion. T'he present discussion seems to show definitely that the low about 
wings of small aspeet ratio is not purely a stalling phenomenon, but is amenable to treatment 
by the perfeet fluid theories. A non-linear theory is required, however, whose basie difference 
[rom previous wing theories consists in the fact that the trailing vortices are assumed to 
leave the plate at an angle © which is determined approximately from the Helmholtz 
vortex laws. For the limiting case of zero aspeet ratio, this theory gives a „Newtonian*“ 
type law for the normal force coeffieient, namely Cy=2sin’a. For the limiting case of 
infinite aspeet ratio, it gives the result of the Prandtl lifting-line theory, Uüy=27sina-cosa. 
For intermediate aspeet ratios, the theory gives the transition from the first type flow to the 
second. The agreement with experiments is quite good for very small aspeet ratios and fair 
for the larger ones. 


B. Appendix: Rectangular Wing of Zero Aspect Ratio (finite span & infinite chord). 


Considering a wing of very large chord t where t>x, the vortieity distribution 
2 HK e | x |tj2 oives v>v., as t>x For a wine of verv small 
m+ 2 "Tr ie Tee 5 ! 

aspect ratio (k=0), the bound vortieity has then essentially a constant distribution along 
the chord. We may assume, however, that the region is still singly eonneeted since we 
still have a definite leadine and trailinge edzee. Thus in ealeulating the force on the wing, 
the Bernoulli constants of the region inside and outside the vortex wake are still the 
same and the formula for the force on a vortex F=oYVJI’' may still be used. 


For our vortex system replacing the wing, \ 
we assume again a system of bound vortices lyıng ei; 5 


[unetion y = y, 


in the wing, with straisht trailing vortices coming AN 0’ 

off at each end ofthe span. The bound vortices | / " 

are assumed to have a constant strength along d e 

the span and the same value y, at every point | > — 
on the chord. "The trailing vortices agaın lie at g’ wg 
an undetermined inelination © with respect to 8’ 


the plate. 
The veloeity induced by this vortex system 
can be written down at once.  Evidently the 
eontributions of the bound vortices normal to the ar % 


plate cancel at every point. However, they induce «— 
. ag 
a velocity tangential to the plate, namely 5 Pa 
n "e: er x 
on the upper surface and , on the lower.- aan > 
2 —— 
The induced veloeity due to the trailing vortices er 
Is the same as ifthey were arranged as starting . 
from a vertical plane through the origin and 
had a strength per unit height of .” .. Forthe 1 > - 
sın eo 
eontribution normal to the wine of a vortex ee. 
segment AB is exactly the same as if it were a u 
placed along A’b’. Thus we get the equivalent | ee < 
trailing vortex system shown, consisting of two De 
semi-infinite walls of vortices. The downwaslı Zu Bez = 
due to this semi-infinite system of vortices at any [ARssıs] a 
point along O0’ is Just half that ofa doubly infinite 
system. This latter has a uniform downwash ig. 15. 
Yo 


rn everywhere between the two walls of vortices, and zero outside. Thus the trailine vortex 
N ( 


system of our wing has a uniform downwash at every point of the span. This downwash is 
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perpendieular to the vortices and thus the component of the induced veloeity along the plate 
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EI r- -sin O the eomponent normal to the plate is w -c08 9. 


2. sın 0) 2 2.sın O 
If the wing moved through the air with a velocity U at an angle of attack a, the 
component of the veloeity normal to the wing is Usina. Our boundary condition is now 


as before, that the velocity induced by the vortices normal to the wing is just equal to 
Usina, 1. e. 


5 .recot I !’.sın a 


or 





2U.sına-tan®. 


In order to determine the angle © as a funetion of a, it is necessary to have expressions 
for the velocıty field due to the vortex system at an arbitrary point. This is determined by 
the use of the Biot-Savart Law, similarly as for the case of the finite wing. Letting 


ı = veloceitv along chord «x ( ); = veloeity inside vortex wake 
"== velocity along span y ( )„ = veloeity outside vortex wake 


MR veloceitv along normal to wing z 


the results obtained are 








| h 
M; Yo for Y , 
h 
" 0 7 z 
ER 
eot Oo) | | ) Y) ee r f 
’o -108 or allı 
! Se Fr Ö , | = „ 
| ) | „) +? 
cot I | 2 2 I 
; "pi 7 tan ] tan 7 | fo1 yI<5 
| 319 29 
',.cot 9 SEE. ... BE . 
Ta 5 tan u’ tan 7 for y|>5 
Jr ) „ ) 





The angle © is determined by the condition that the vortex lines follow the streamlines. 
As an approximation, we shall assume that the vortex sheets remain in vertical planes through 
the sides of the wing, i. e. we negleet the »-component of the veloeity normal to the vortex 
sheet. The inelination of the vortex filaments is then given by the equation 


mean lv b]2 


tan oO a 
> Hmean |» DE 
where 
WW U.sın«a| 2 2] 
W mean |y — bi? HU. sına 1 -tan | 
y y | JT I | 
u; u 


Hincan | bir 9 4 l’.cos a l/. cos a 1-l/.sın a -tan@. 
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Thus 
U.sına 2 2 | 
” | --+tan' 
7 h| 
tan 


l/.cosa-+ l-sina-tan © 


or 


« 
Enz 


eota , eota f 2 
tan = un ae I+-2tan’a!l : tan 


> 2 b|' 

The trailing vortices thus are not actualiy straight lines but lie along eurved paths which, at 
a great distance behind the wing, are parallel to the wing. This contradietion with our 
original assumptions arıses out of the fact that we have in part violated the Helmholtz 
condition that the vortices follow the fluid partieles. It will be remembered that the Prandt| 
wing theory which assumes the trailing vortices to be carried along the direetion of the 
undisturbed low sımilarly violates the Helmholtz conditions. It thus seems justifiable to 
retain the assumption that our trailing vortices follow straight lines at some efleetive 
inelination © with respect to the plate. For this effeetive inelination we shall choose the 
inelination near z=0, since the vortices near the plate are the most effeetive in indueing 
the normal veloeity at the plate. In other words, we shall assume that for our wing of zero 
aspeet ratio 

cota , cota 
tan ©) ar - vyı + zta’, 


-_ 


The normal force actinge on the wing can now be determined from the formula for the 
force on a vortex. We obtain for the normal force per unit length along the chord 


N o(lU-cosa ru,)'Yy DB 


jo ’ 


where 


tt, Induced u-velocity at plate = 


y„_ 2U-sina-tan 
and 
N 
( ' | 
=o0U°.6 





I/\ —=2mnal,. 





Thus we again obtain the same result as before for the normal force eoeffieient. In 
addition, we now have a pieture of the velocity field corresponding to the wing of zero aspeect 
ratio. This seems to indieate that our violation of the Helmholtz condition for the cross: 
low v» has resulted in an internal contradietion eoncerning the shape of the vortex lines. It 
seems hopeless to try to caleulate the induced veloeities corresponding to vortices Iying along 
a curved surface as is actuallv the case. The most we might expeet to do in a refinement 
of this theory is to obtain a better approximation for the effective angle © at which a system 
of equivalent straight trailing vortices leaves the wing. 363 
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Über die zahlenmäßsige Beherrschung und Anwendung einiger 
den Besselschen verwandten Funktionen nebst Bemerkungen 


zum Gebiet der Besselfunktionen '. 
Von Reinhard Miller in Berlin. 


Übersicht. Die vorliereende Arbeit bringt Zahlentafeln für die Integrale der Besselschen 
und Hankelschen Funktionen erster Art, nullter Ordnung für rein imagınäre Argumente 
dy 


i.w, sowie für eine Funktion 9 (#*) nebst Ableitung 5, die in der Theorie der Trag- 


’ da 
Nürelschwingungen auftritt. 

Ks wird eine erößere Anzahl von Funktionen ın Inteeralform anzereben, die mit Hilfe 
dieser Zahlentafeln zugänglich geworden ist. So ergibt sich u. a. die Klasse der Struveschen 
Funktionen für halbzahlige Parameter 2r =0,1,2... 

erner wird die Lösung der inhomogenen Besselschen Differentialgeleichung für eine 
Störungsfunktion in Form einer Potenzreihe mit ganzen positiven und negativen Exponenten 
auf die Zahlentafeln zurückgeführt und gezeigt, daß dies für die praktische Durchrechnung 
des Problems der erzwungenen Schwingungen elastischer Körper von zylindrischer Gestalt 
von Nutzen Ist, 

Schließlich werden einige Typen von Differentialgleichungen aufgestellt, die sich auf 
die Besselsche Gleichung zurückführen lassen. 


Zusammenstellung der verwendeten Bezeichnungen. 


(a \ sine. : er Dr Er 

nv“ u . (2) \ sin“ (af) cos“ (P #) sind (z sin + 0, f)ecos’:(zsint+Öd,N)dt . 1V, Gl. (17), 
fr. (2) Fa BEI RE. Evo re EI 

hu) ia te init, ER 2 u er EA er a 


wis) = 2,(i09)=0,J,19) +0, 2,” (e) 


Il. Zurückführung eines Integrals aus der Theorie der Tragflügelschwingungen auf Besselsche 
Funktionen. In einer im Flugtechnischen Institut der Technischen Hochschule Berlin dureh- 
veführten Untersuchung über Tragflügelschwingungen (v. Borbely*) trat das Integral 


JH 
Sn @ 
Det) \ gr „de. Hr v) 
J ur 5) he 
tl 
auf. Ergänzt man g,(a@) dureh 
7 
COS (It) 
(rt) Ir tt r) \ „do 


] ze? 4 m? 


und integriert die Funktion 9.) -+-i gs (a) über den ersten Kreisquadranten der komplexen 


ie!) 


bene. so findet män die für = 0 eeltende Beziehung 


Dissertation, genehmigt von der Fakultät für Allgemeine Wissenschaften der Techn. Hochschule Berlin, 
reeranreen aus einer im Flugtechn. Institut (Prof. Dr.-Ing. H. Wagner) durchgeführten Arbeit. 


\ \ er]. diese YA itsehın ıft Bd IS (1938), He to 
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Wir betrachten etwas allgemeiner 


plz2)—|\e-zsint di Re = 5 | Be er ER 


y(z) befriedigt die Differentialgleichung 


j 
m (U) 
RT: I dw | Q 2 
| ’- 1) 
dz? 2 dz u F® 
da | I 
Durch „Variation der Konstanten“ erhält man als Lösung von (2) 
i Pr. | -7 gr T (1) =. | Ft dt. =. . .) 
(2) li) \i— H, Sa, PtzH, iz)\J, id): En > A 
2 MM 5 u m 
Hieraus folet durch Differentiation 
di n |, TREE 
y ide), \’, TAGORE 9 H"’(iz)\J GH) dE  . k4). 
m m 
Für die Funktionen (3) und (4) findet man ım Anhang Tafeln für 2° x: fürz=i. gıbt es 


bereits Zahlentafeln ’?). 


Il. Mehrfache Quadraturen über die Besselschen Funktionen. Setzt man ın der Bessel- 
schen Differentialgleichung 


mw I dm | y" 
da? 2 dz > | ) (.)) 
den Parameter r=0, so bekommt man 
di dw\. 
14 rzw=VU DE DR SP DSB an RG De 
dz 


dz 


Durch zunächst zweimaliges Integrieren von (6) erhält man 


Im r 
2 z|wdz ((wdzdz ar Ge er Ger eo u; ° 
dz 
zw (wdz z((wdzdz2+2[(((wdzdz2dz2 ....:.:...68 
USW. 
Die Lösung von (6) entsteht dureh lineare Kombination von ww ==.) (2), IH” iz). Für 
. y be er d m (1 4 r » er) 
die Funktionen ww und für N , J, (2), - HH," (2) sınd Zahlentafeln für z2=.x und 
(di 2 
z=i.x vorhanden. Die Funktionen (wdz=/(J,()dz, ({H”(z2)dz sind für z=.x bei 
Watson tabelliert’), bzw. auf andere Tafeln dieses Werkes zurückgeführt; für 2. ver: 


öffentliche ich Zahlentafeln am Schluß der Arbeit. 


Wiederholte Quadraturen über die Besselschen und Hankelschen Funktionen nullter 
und erster Ordnung kann man also nach (7) und (8) auf tabellierte Funktionen zurückführen. 
Allgemein lassen sich sogar ı-fache Quadraturen über Z,(z) durch tabellierte Funktionen 
ausdrücken. 


_=.. DR. u sa: 2=x oder 2=i&, w= 2,(2)=c,J,(2) +c, H,“' (2) 


befriedigt nämlich die Differentialgeleichung 


d | dm y? 
Maid ERFRUN . gie Er u 
dz dz | g° 
2) Watson: A Treatise on the theory of Bessel Funetions, Cambridge, 1922; im folgenden als „Watson 
hbezeiehnet. 
T d T ‚[7 

r(ix) (Ju (x) -iH,(Xx)), 7 WÜR)) Jı (x) -+il H, (x) 1). 
2 AEX 2 2 


H,ix) und H, (x) sind bei Watson tabelliert, Tafel I, S. 666/7 MT. 


3) Watson, 8. 722. 
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Durch Integrration von (9) erhält man 







d 1 


z p 2 \ w(dz rn \ \ wdzdz-»? \ = dz. | (10). 


7 
































I. ? 
Man kann aus (10) \ _ dz für »—=1,2,.. berechnen, falls man (w,dz und ((w,dzdz 





schon kennt 








Ferner eilt bekanntlich 


oder ınteeriert 


er tn, 
\w,_,dz F\w,+.,dz 2» \ Eier DD an Wa 


Aus (11) erhält man dureh Multiplikation mit z und nachfolgendes Integrieren 


‚dz—((w,_,dzdz+z([w,.,dz2— [[w,.,dz2dz=2v(w,dz. . . (13). 
Wiederholte Quadraturen über »*, und w, waren auf Zahlentafeln zurückgeführt. Für v = | 
5 | mw R | 

liefert daher Gl. (10) \  dz. Darauf erhält man aus (12) (w,dz und aus (13) ([w,dzdz. 


Setzt man v=1,2,.., so bekommt man wiederholte Quadraturen über »,(z).. Um eine 
-malige Quadratur über w,„(z) zu erhalten, multipliziere man (11) mit 2”! und integriere 
danach mehrmals partiell, wie das für «==2 in Gl. (13) geschah. 


Nerative v erledigen sich durch die Beziehung 


. n v R er. 
w_„(2)=(— 1)" w,„(2). 


Ill. Das Integral | 27 0,(2)dz und die Integrale [cosz.J,(2)dz, [sınzJ,(z)dz für 
3% 1.0.1 
w„(2)=c,J,(2)+ c, H,'" (2), =, ,-ELRT,., s= a md sie 


Für positive ganze 4 kann man | 24 w,dz durch 4-malige partielle Integration?) zurückführen 
auf wiederholte Quadraturen über »v, (2), die in II behandelt wurden. Man kann für /=1,2,. 
auch zunächst der Rekursionsformel’) 


Ar'w,dz=(t —I)fAtTwdz+tztiw,ı, tA—r)Aw, ... 0. (1) 


folgen, die das linksstehende Integral für ungerade / auf [ıw,dz, für gerade / auf 


(zw,dz2=2 |w,„de— ((w,d2dz 
zurückführt und erst jetzt nach II weiter verfahren. Im Falle A + v» ist (14) zu benutzen. 


Der Fall positiver ganzzahliger 4 ıst hiermit erledigt. 


ser I, u. * . \ . .. r. 
(12) führt \  dz für v-1,2,.. auf tabellierte Funktionen zurück. Wir nehmen vor- 
übergehend an, die Funktion | _" dz sei tabelliert. Division von (11) durch z und nachfolgendes 


. . » Z 4 Do. .. ® da \ > 
Integrieren liefert \ _ dz für »—1,2,.., während sich | _?dz aus der Gleichung 


az a. eh 


/ Fl; gr: 


für = 1 ergibt. Dividiert man jetzt (11) durch 2° und integriert danach, so ergibt sich 
1° En ‚ , . Li - . . .. 
\ , dz für v»—=1,2,.., während man dieses Integral für v=0 aus (15) mit 4=2 erhält. 











I", 


Fährt man so fort, so zeigt sich, daß sich das Integral \ -ds für k=1,2,3,... auf die 





> 


tw 










. w, 


Tabellierungz von \ dz reduziert. 


') Die partielle Integration ist so auszuführen, daß sieh der Exponent der Potenz bei jedem Schritt um l erniedrigt. 
Siehe z. B. Watson, S. 132 bis 133. 
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Diese Tabelle ist aber zur Zeit noch nicht berechnet. Daher werde noch untersucht. was 


h zo. 1, en 
die vorhandenen Zahlentafeln leisten. (11) zeigt, daß sich \  dz für vu, a«-+l,.., d.h. 
jr rp\ ei ie o . ; wu 
für u» (u,» ganz) durch Tabellen ausdrücken läßt. Insbesondere ist dies also bei \ ——- dz 
y 1] nd 4 (d \ er p) m, 
der Fall. Setzt man hier «= 1, so liefert (15) für A=1 das Integral \ !dz. Man setze 


jetzt in (11) »=1, dividiere durch 2° und integriere danach, sodann setze man » = 2, dividiere 


dureh 2° und integriere, bis man die »-te derartige Gleichung nach Division durch z”*+' inte- 
. ie m, . oo. =, R i 
eriert. Man erhält so \ zu + „dz. Setzt man hier v 1, so liefert (15) für 1=5 das Integral 


m FRE 2. ’ 
\— da. Dividiert man nun (11) für v | dureh z* und integriert dann,..., bis man die 


iw 


(IR 


v-te derartige Gleichung nach Division durch z”*+* integriert, so erhält man \ 57:48. Die 
a, mw, Th 
Integrale \ = 73 „dz und nach früherem auch \ - y dz für k=0,1,2,.. sind somit auf heute 
berechnete Zahlentafeln zurückgeführt. 
Integrale der Form 
nt s+if Ir 
\ zi \ w,dzdz 2: \ w,dz , \ zu tiw,dz 


beherrscht man ebenfalls. (4 —+1=-0.) Unsere Integrale sind auch bei jeder Funktion f(z), 
die durch eine Reihe positiver und negativer Potenzen approximiert werden kann, für das Integral 


(f(z) w,(z2) dz 


von Nutzen. Bisweilen führt aber ein Kunstgriff schneller zum Ziel. Hierfür als Beispiel 


F(t)=\cosz J, (2) dz. 
0 
Mit 
G(t)=|\sınz J,(2)dz 
0 
wird 
! 
F+iG=eit\cos(2 HM) J,(2)dz+eiti\sin(—tH)J,(2)dz=et (fh) High). 
" 0" 
f und y genügen dem System f’+-f = J,, g+f=0 


ft)=C,sint+(,cost+tJ, (N. 


Wegen f(0) =0 und f’(0) = J,(0) ıst C, = (0,0, 


{ ! tt 
(th) \tJ, (Hal t\J, dBat+\\J,(Hatdt. 


0 0" u 


Hiermit ist alles auf Tabellen zurückgeführt. "Man findet schließlich 


! 
\cos2J,(@)d2e=t[costJ,(d)+sintJ,()). . .». .» . .  « ABB 
0 
2 ! ! 
(Übrigens sind auch \cos2.J,(z)dz==F,(tf) und \sinzJ,(@)dz=- G,() für v=..- 1,0,1.. 
U ( 
auf bekannte Funktionen zurückführbar. Man kennt F',, @, und erhält wegen .J, : J,’ dureh 


partielle Integration auch F,, @,. Aus der Funktionalgleichung J,_,(2) - J,4,(2)=2.J, (2) 
gewinnt man leicht die Beziehungen 


F,., d)=F,_,() +2(J, (0) — @, (fl) — cstJ,(t)), 
G,., d)=G,-,(d +2 (#,(t) -sıntJ,()). 
Hier hat man nur der Reihe nach »=1,2.. zu setzen. Die Integrale 
[cosz H," (z)dz, (sinz H,"” (z)dz 


erledigt man ganz entsprechend, jedoch ist die Singularität von H,“” bei z2=0 zu beachten.) 
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IV. Weitere Funktionen in Integralform, die mit Hilfe der Zahlentafeln des Anhangs zu- 
gänglich geworden sind. Zunächst sei angeführt 








F sa „®, (2) \ sin“ (af) cos": (PB sind: (zsint+ 9, cos’ (zsint+0,h) dt . . (17), 


Ad a v, en ıı.D. ) gi R Am 


/ N 2 .. . u iw 


N Iy. 


Man forme F in eine Summe von Ausdrücken 





piihrt rt sınt) fl f (eos ft N isın PH" etmzsint ff 


ı () 


um. Ersetzt man © mz dureh eine neue Variable, so hat man ledielich 


fr, u (2) = \ sın" t cos" teesint dt EAU Fe ee 


u 
zu betrachten ®). 
Ist in (15) « gerade, so erhält man durch Ersetzen von eos?* dureh (1 — sin?)* die Form f,.. (2). 
von der später zu sprechen sein wird. Ist « ungerade, so kommt man durch die gleiche Um- 
formung auf f,.,(2); in diesem Falle führen einige partielle Integrationen nach # zum Ziel. 
Somit ist (15) zurückgeführt auf die Behandlung von f,.,. (2): 
Unter Berücksichtigung von 
ı.J,(i2) EEE 2 DE a Ve A we 


bekommt man 


eitidgt=n Ju —- ol) .: - 2 2. 5 Bon ar. 


Vel. I, Gl. (1) und Fußnote 2. 
Hieraus erhält man dureh Differentiation 


| | do 
\ sin fe:sint di id, (iz) ae 
dz 
und weıter 
| I dag | 
\sin’te sıntidf 1), (iz) ae iJ, (iz)) (2) + r a |} 
zZ 2 G% 2 


ı 


Man kann nämlich bei den auf den rechten Seiten der Gleichungen stehenden Funktionen die 
höheren Ableitungen mit Hilfe ihrer Differentialgleiehungen (2) und (5) ausmerzen. Durch 
forteesetztes Differenzieren von (20) läßt sich also 


Lei. ee res 
il 
und somit aueh (17) und (18) für 2=.r und z==i.x auf tabellierte Funktionen zurückführen. 


Wegen der Periodizität des Sinus könnte auch in (23), bzw. in (17), (18) die obere 


(‚renze des Integrals an Stelle von, ein ganzzahliges Vielfaches von sein, ohne daß sich 


DNS 


etwas Wesentliches ändert. 
Wir notieren noch die sieh durch partielle Integration ergebende Beziehung 
| | | 
\t/coste sint .]f r 4 (2) e | ; - { ; : i A (24). 


M 
Man beherrscht die Integrale (23) auch für negative ganzzahlige v. Man bekommt sie durch 
wiederholtes Integrieren von (20). Auf der rechten Seite der Gleichung entstehen immer 
tabellierte Funktionen, vorausgesetzt allerdings, daß noch eine Tabelle für 








Piz) 









‚befriedigt die Differentialeleiehung (ze) + (vr +4 I) or (z w)’ vu vr, seien reell> 0, 
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gerechnet wird. Mehrfache Quadraturen über ® (z) wären ja durch Integration der mit z 
Salz. ei dq es 
multiplizierten Gl. (2) auf P (z), 9 (2), re zurückführbar, wie schon in II entsprechend an der 
; ( 
Besselschen Gl. (6), (7), (8) gezeigt wurde. So gilt beispielsweise 


| > -sint 
I dt--D(z) 


u 


Me 


\ zsint Sincz sin) 


| 
dt = \[P () HDi 2)| dz. 


sın? / 2. 
u 0 
is Ist 
71 Dr ML. z 
PD (z) Fe), \i = N dHdE +2\ FISHER UDBE . 2. 5» 
ü 5 
mit 
F(e)=\J,t Hd: 
Plz) -—P, (—z)=nF kl) (26). 
Aus (26) oder (20) folgt 
o(2) Tp|(—2) tJ,(t2) (2) 
po (2) -p (-2) niJ,(i2) (28) 
Ferner ist bekanntlich 
-Wai=sisetnhtt+r di --; » 2 ei MM. 
(u 
Setzt man nach (17) 
N. \cos (zsint+rt)dt=y,(z), 
0 
so gilt für r=..—1,0,1, 
v,„@)+(— 1’ yv_,@)=(-1N”’aJ,le -» -» » 2 202020. (BD) 


Wir werfen nun einen Blick auf die Struveschen Funktionen S,(2). Sie sind definiert 
durch die Gleichung 


> v HE 
RM |, 
J,(@)+iS,(): lesimteossstat ...2.2.20.0. 080. 
Tv Ft, 
0 
Unter Benutzung von (18) wird j 
|» Y 
2(5) 
J,(2) +tS,(2) U ı % | Ps Er 5  ° 


| 
Val r +5] 
Die Klasse der Struveschen Funktionen’) ist somit jetzt tabellenmäßig erfaßt für z=.r, 
ir and Br => 1 2,.. 

Es sind jetzt einige Betrachtungen im Reellen durchzuführen. Aus (15) folgt mit dem 
Moivreschen Lehrsatz, daß man die Funktion 


ci} 
\sın (vr #) eXsint dt 
IE 
”) Näheres über die Struve schen Funktionen findet man z.B. bei Watson, 8. 328 1V Wevrich:Die Zylinder- 
funktionen und ihre Anwendungen, Teubner, 1937, S. 89. — Über die Bedeutung der Struveschen Funktionen in der 
Physik, vel. Watson, 8.417. — Struve: Annalen der Physik und Chemie, 17 (1882), S. 1010 bis 1011. Mae Lachlan: 


Bessel Funetions for Engineers, Oxford, 1934, S. 65ff. — Lord Rayleigh: Theory of sound, (1896), Bd. 11, S. 169. 
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für ganzzahlige » tabellenmäßig beherrscht. Nun bestehen die Gleichungen 


j BE 1 sen 
\ cos(actgi)tg?st!|, Kae \e-esintsin[2s+l)fldt . . .. (3), 
Y | er i 7] 
\ reg \e asint sın (2s t I)tldt: | BI ee 


Sie sind abgeleitet bei Theisinger, Monatshefte für Mathematik und Physik, 24 (1913) 
S. 337 bis 311. 


Kerner eilt für Inteeration länes der reellen Achse 


I 


(“ k\’ 
| Be | | da ya 2 a. 
cos| 5 vr), (ax) sin | 5 v7) S, (a re 1: | | fuvt-ah). . (BB). 
Pie 
| Ba 
Vorauss. Re (k)>V, -a<RBRe(r)<2, a>—, 
| (3) 
J,lax)dx - ) h 
N - . Ba te 7 er a 
7 Ai u k° &L k: . | y 
" I ( 2 | 
' | 
Vorauss. d, Re (k) 0, Re (») > 1 
. 16 IR d 
N, iO) d oc rt“ ” > 


Y‚,„(ak)+i]|J,(ak) nee (37). 


+ 2 tk’+t'cos(en). 














Vorauss. Re (v)« 


Hierbei ıst 
T 
> 


jer at JPG) =K,(E) 


vesetzt. Die Gl. (35) bis (37) findet man bei Watson, S. 425 bis 26, abgeleitet. Dort sind 
noch mehr derartige Integrale angegeben, die auf die Funktionen (15) führen. Dasselbe 
eilt von folgenden Arbeiten: Meijer, ©. S.: Integraldarstellungen für Lommelsche und 
Struvesche Funktionen. Akad. Wetensch. Amsterdam, Proc. 38 (1935), S. 628 bis 634, 744 
bis 749; Proc. 39 (1936), 519 bis 527. Mae Lachlan and Meyers: Phil. Mag. VII, s. 21 (1936), 
125 bis 448. | 

Schließlich integrieren wir noch die Funktion 


. "T * % e ur PR a 
175 HrGO=\ = dv &,v reell, E>VO Br (39), 
2 Jyv’—1 4 
| 
mit der Substitution » von OÖ bis x und erhalten 
sın 
X | . 
Be HN.’ (iA)dAE ” 7 / » = f) 0) 
1  NKkö)de = \e Ed, A Be «Er Th 489). 
N! Mm 


Wiederholtes Differenzieren bzw. Integrieren‘) von Gl. (39) nach .e ergibt das Integral 


c 


|») 


x 
’ En, : n 
h,()=\sin”t-e sintdt, > Se 
M 
sS) Auf dieselbe Art lassen sieh aus Gl. (29) weitere Integrale gewinnen, deren Anzahl noch dureh Benutzung des 
Systems (26) bis (28) und (30) vergrößert werden kann, 
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denn links stehen tabellierte Funktionen. Außerdem zeigt (39), daß 


X 
Er =. 7 
im iS Hrüddl=Z....... 2 A 
xo-n, _ Pit 
S u 


ist. Von dieser Beziehung werden wir später Gebrauch zu machen haben. 


V. Die Funktionen f,,.(z2) und %(z) in der komplexen Ebene. Aus ler Definitions- 
eleichung für fi,« (2) X, «—=0,1,2,..) erhält man die Reihenentwicklung dieser ganzen tran- 
szendenten Funktion. 


v4 (Ab 
\ u. 
fr, (2) ‚\sin” + ae: ee ) (42). 
u Ir ” . 
m 0 0 


Durch Auswertung”) des Integrals in (42) folgt 





im 
77 u! \ (v -- m)! | 24 | (2 
Iv,u(2) Pratly, ‚mv +utrm, m! (0) 
9 -m=V .) . .) | 
Demnach gilt 
m 
| a ylr) 
ols)=f,..il%) = en ; ee . (444) 
E m ,\' 
im 0 (5 
oder 
\ev 
x | L 
7 r\ ) Y 22’ + 1 M 
G (2) 4) \ \2 \ .) - oO f 5} (4»). 
2 a KH" PAAR I1-3-5... (Zr -+1)]’ 


T 
5J,(tz). Setzt man 


_ 


Die erste Summe in (45) ist die Reihenentwicklung für die Funktion 


7 
p (2) ss, lu): : ee ER, 
so wird 
T r 
oi- Des t+el : 5 


(Die Funktion 6 (z) wird bei Watson mit ©. H,(- %z) bezeichnet.) Ferner erhält man aus 


der Definitionsgleichung für fr, .« (2) die Entwicklung 


es 
Y ‚m 


» N D A » . " 
f5, 7 (x 7 l Y) _— \ N  fv 4 m, HH (? Y) R j i ; ; . . ) IS). 


Setzt man 


so ergibt Gl. (48) in Verbindung mit der Definitionsgleichung für die Funktion F die 
Beziehungen | 


on 
vw.,Mm 
\'a 1,1,0,0 19 
Ur, ul, 9Y): ; tr rm) - RER 9 Me 
ee Z ® 
m 0 
SL 
vw_,Mm 
\ 4 a... 3 
Ur, u(X, 9%) Eur BD 2 Bhurl 2.5  A 00 Se 
TEE ° 
m 0 


», Vegl.z.B. Jahnke-Emde: Funktionentafeln, 2. Auflage, (1933), S. 95. 
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y =180° 








Bild 2. 


Bild des Halbstrahls ı 





NITH 


in der w-Ebene 
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: | ’ j ER f 
Für die unendlichen Reihen (48), (49), (50) ist die Potenzentwicklung von —, eX eine Majorante. 


In IV wurde gezeigt, daß die Funktionen f und F unter den Summenzeichen auf den rechten 
Seiten dieser Gleiehungen für reelle und für rein imaginäre Argumente tabellenmäßig zu- 
eängelich sind. Aus (48) erhält man noch 


(— 2) , . 
7 (.r rımy) Ft K ıı) \ y Ivo! ı ı/) . . u . . (>1) 
rn: . 
oder 
\ ”( ie)” dd’ Q (0 „) (59) 
. 3 ar AR 
AA ıy) y N dy’ 


) u 
Die Abbildung wg (2) zeigen die Bilder la, Ib (2=reiY; w=u-+ir). Durch Spiegelung 


an der reellen Achse entsteht ein volles Exemplar der ı-Ebene, die man sich längs der 
nerativen reellen Halbachse aufgeschnitten denken mag. 




















v 
h 
50 
V 
| | 
-100 00 200, u | 
1/ 
/ 
/ 
v 955 
05 sy 
/ 
/ 
a5 
-100 | 05 ol / 0° et 
14 gi mi 
80° 
>» 
-05 
Adsi 
nn 
Bild des Kreises r =6 in der ıw-Ebene 
Bild 3a. Bild 3b. 


Den Halbstrahlen y = const entsprechen Spiralen mit unendlich vielen Windungen. Für 
0° <= 90° werden die Windungen mit wachsendem r kleiner und kleiner, für 90° < y < 180° 
erößer und größer. Ein Beispiel gibt Bild 2 (Zahlenwerte aus Watson entnommen). Bilder 
der Vollgeraden durch 2=0 sind Spiralen, die sich um O0 und & winden (Ausnahmen: 2° ., 
2=tYy). i 

Die Kreise r = const verlaufen für kleine Werte von r (etwa für O<r<4) nur im 
Il. Exemplar der w-Ebene. Für r > 4,3 tritt bereits Spiralbildung ein; mit wachsendem r nimmt 
die Windungszahl der Spiralen mehr und mehr zu, bleibt aber stets endlich. Ein Beispiel zeigt 
Bild 3a, 3b. (Bei Bild 3a ıst der Kurventeil im ersten Quadranten der w-Ebene zwecks 
besserer Veranschaulichung etwas übertrieben gezeichnet.) 


VI. Anwendung der inhomogenen Besselschen Differentialgleichung bei dem Problem der 
erzwungenen Schwingungen elastischer Körper von zylinderartiger Gestalt. Bekanntlich '") löst 
man mit Besselschen Funktionen die Aufgabe, die freien Schwingungen zylinderartiger 
Körper zu berechnen. Erwähnt seien Kreismembran, kreisförmige oder kreisringförmige Platte, 
Zylinder endlicher bzw. unendlieher Länge usw. Für die Behandlung der entsprechenden 
erzwungenen Schwingungen findet man in der Literatur die Entwieklung der äußeren Kraft 
(x, y,.#)'') sowie der Lösung V (x, 9,#) nach den Eigenfunktionen des freien Schwingungs- 
problems angegeben. Es sei hier ein anderer Weg angedeutet, der in manchen Fällen die 
Rechenarbeit vermindert. 


10, Vel. z. B. Handbuch der Physik von Geiger-Scheel, Bd. V1. 
11) Wir betrachten hier der Einfachheit halber das zweidimensionale Problem. 
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Man entwickle © in eine Reihe 





y ı ! 
DVI. Y, f) > 4, (X, y)e I 
(k) 


und setze für V an 


iv, t 


V ".(a,y)e Ft 


ı(w,M, f) Vi, Y, !) Bo 


(k) 


Hierbei sei die Lösung des freien Schwingungsproblems V*(#, », f) bereits bekannt. Geht man 


nun zu Polarkoordinaten (r,g) über, so ist die Entwicklung 


a 


4, (#, 9) 4x (r,Y) yulld (vr) ein“ 


durchzuführen, während man für ”, den Ansatz 


Apr P} PR) (r) ei 
(1) 
macht. Die Funktion Ru’) (r) muß der inhomogenen Besselschen Differentialgleichung 
genügen, wobei, abgesehen von Konstanten, die Funktion 9.) (r) als Störungsfunktion auf- 
tritt. Man approximiere gu") (r) durch eine Summe von Potenzen mit ganzzahligen negativen 


m 
und positiven Exponenten a; ri. Die Lösung der inhomogenen Besselschen Gleichung 
f n 


läßt sieh dann aus Besselschen Funktionen und Inteeralen der Form 
j 


’ır ra “en r we 
\r*Z,G0*)dr, r i Er N > du, 2 


) 
N 
sehr leicht aufbauen '?). (z und v sind eanze Zahlen.) _Derartiee Inteerale konnten in III auf 
tabellierte Funktionen'’) zurückgeführt werden. 


VII. Einige auf die Besselsche Gleichung zurückführbare Differentialgleichungen. Die all- 
gemeinste dureh Punkttransformation in die Besselsche Gleichung überführbare Differential- 
eleichung 2. Ordnung haf Neumer angegeben’). Indessen ist die explizite Auffindung der 
Transformation sehr schwierig. Daher erscheint es nicht zwecklos, einige auf die Besselsche 
Gleichung zurückführbare Differentialgleichungen nebst der zugehörigen Transformation an- 
zugeben. Wir betrachten zuerst 


v+y’oy)+yfa)tyu)y(a)=0 .:. . 2.20.2002. (54). 
Die Bedingungen 
dy 
dy i 
(m) ee ee rt re Sa er 
yviy) 
.) ” F R . 15) +) Y » 
ga) + e@)‘ FR er 
wober F(r)= | f(w) de ıst und v» einen Parameter bedeutet, sind hinreichend dafür, daß sich 


(Dh) ın die Besselsche Gleichung transformieren läßt. Die Substitution ist dann 


1 
Fix) . \ d 7 
& f Aa : ) — p 
/ 
. > du /, (8) obere , ' , ' 
Die alleemeine Lösunz von (54) ist e z . . für das Ja Vorzeichen in (2”). Bei- 
iyitce) tere 
ji - 
spielsweise hat die Differentialgleichung 
a : Pr . er ä 
7 y (ige + cig x) + cig” zz (sın? x v)y In y 0 
. 7 “ « « 


die Form (54) und erfüllt die Bedingungen (1’), (2”). Dasselbe gilt von 
Y (+ p2AX y*) /2 Y —{) \. Y — eONnSt. 


el. z.B. Wevrich: Die Zylinderfunktionen und ihre Anwendungen, (1937), S. 87. 
13) Hinsichtlich negativer * sind die in III erwähnten Einsehränkungen zu beachten. 


14) Journal für die reine und angewandte Mathematik, Bd. 173, Heft 3, (1935), S. 131 bis 132, Satz (4), S. 128, 


’ 





7. angew. Math. Mech. Pr . 20: R ' = 
19 Müller, Zahlenmäßige Beherrschung und Anwendung Besselscher Funktionen 47 


3d. 19 Nr. 1 Febr. 1939 . En ee 





Ferner sei die Differentialgleichung '’) 
’ Ka Ei 2.92 ei 4 
"tHtW+nl.t+Fr@)t+RW-k-lt1- Gt F,@)=0 .... (65) 


aneeführt. Unter den Voraussetzungen 


» » / 67, ) | ’ 

f, () f: (y) Ep (k,» Parameter; k=-=-0) . or 
ist die Lösung von (99) 

nn E - \F,(x)dx 
e iM —const-e ° [7, (e)]". 
E obere 4 f b be d er 
Für das Vorzeichen in 55) gilt z =... 
untere iX 


Es seien jetzt in der homogenen Differentialgleichung 2. Ordnung 


n"+n’o,@)wo.l@)=0. . . .: .: 2.2.2.2...) 
y, und g, vorgegebene Funktionen, 9,320. Setzt man 

a, » +»... Tu Ge ra ee ;;° 

EEE ARTE a. ‚79% f ; ;; ° 


wobei F noch zu definieren ist, so handelt es sich also um die Differentialeleichung 
"+ f )+wf, 2) twF(f,,f)=0 .:. . .: 20.2 0.20.(56) 


Für das Funktional 


\ al f; fe\,. | 'G fs, 1 ) 
2 Ay 2\(ı —2 z )da — \(hı —2 7 )d: ) 
Fe m 
1 f; | | 
wobei » einen Parameter (» ==0) bedeutet, ist die Lösung von (56) 
Zu | 
wie, 1,\&:' 
Hierbei ist 
\ 24 fr, L \dz 
us ’ı Dez ° 1 dz 
Für das Funktional 
fr fa fa E 
f. 2 f: d \lrı- 2, )d: N. fr 2 5)a: 
1 1 


ist die Lösung von (56) 
w=(,7,(;), 


wobeı 


fo Er: f; 
—-dz I —ıdz 
C, pP hi ; R - ee Hi n)‘ 


dz 

gesetzt ist. Den Übergang von (5’) zu (56) vermitteln die Gl. (6’), (7°). Hierbei ist (7’) bei 
gezebenen y, (2), 9; (2) eine Integrodifferentialgleichung für f, (2), die selbst nach willkürlieh 
gewähltem » immer lösbar, dabei aber recht kompliziert ist. Wir beschränken uns daher 
darauf, Gl. (56) nun unter der Annahme zu behandeln, daß f, und f, vorgegeben sind; der 
Fall f, 0, f,320 werde ausgeschlossen. Dabei ist uns bereits bekannt, daß (56) sich mit 


Besselschen Funktionen integrieren läßt, wenn F mit (57) oder (58) übereinstimmt '®). 
15) Der Fall fs (y y,„k=1, &=x findet sieh schon bei Jahnke-Emde -Funktionentafeln. 


16) Vgl. hierzu Rommel: Mathemat. Annalen XIV, 1879, S. 510 bis 536. Lom mel arbeitete mit der Differential 


eleichung 


ei, Y 
ae(‘) 24 alx zu r 
dv u dy x 


welehe die Lösung 
y(z2)=X (z) | (z) |’ Z, (v (2)) 


hat Dabei sind X (z) und W (z) beliebige Funktionen. Wir geben hier dem Lo m m e|]schen Problem eine neue Fassung. 
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Kür f, a) 0, (2) Fid)-flz) geht (56) über in 


’‚ 


ww -wfF=#6 . oo: ae (9). 










Aus (57) entsteht dann 


(60). 


Vorauss ı 0. 





(seht man von einer beliebigen Funktion f(z) aus, so ist durch (60) ein F(z) definiert. Dann 
gehört dazu eine Differentialgleichung (59), die sich leicht mit Besselfunktionen lösen läßt. 



























Beispielsweise gehen für f(z) ae: (39), (60) über ın die Gleichung 
2 
| 
z; r 2) 
w + (a 2) 50 ı u 


v,a, 4 Konstanten: v0. 


Für f,=f, f. „ entsteht aus (56), (97) 


2 E: 
we rote! wi - Alz+B)' ' we. er a a 00 006,7 © 


A. B,vr=const: A.v==0. 


| 


Wir zeigen jetzt, daß die partielle Differentialgleichung 


x’ +y) 2 ee BT re a ae 


für jede reguläre Funktion y’ auf die Besselsche Gleichung zurückführbar ist. Mit z=rcos y, 
yrsing und dem Ansatz «= ei’’ y(r) erhält man aus (63) die gewöhnliche Differential- 
eleichung 


(l- 7 | dl 7 L 
I = vw(r) + :;|]? 0 DT (64) 
dı“ r dr | ‚U 
Setzt man andererseits ın (D6), (D7) f, (2) .,2-f,(2) f(z), so findet man, daß sich 
7) m de 
"+ erg . rear ae 


leicht auf die Besselsche Gleichung zurückführen läßt, wenn die Lösung f(z) der Inte- 
erodifferentialeleichung 


(z) ” | | D del, —° u 
(2) A I’+lf + e /dz|2y\e? [d: 1 ee re ae ee 
Faß | } = | 
(vr beliebige '”) == 0) 
bekannt ıst. 
Kennt man also eine Funktion f(z), die mit 
u — 
4(2) ur (2) + -— } er 9 Reg oe an Ar Rs A 


die Bedinzung (66) erfüllt, so ıst hiermit die Reduktion von (63) auf die Besselsche Gl. 
erreicht. Nimmt man f beliebig an, so ist durch (66), (67) eine Funktion y definiert, für 
die uns die Lösung von (65) mit Hilfe von Besselschen Funktionen bekannt ist. Als 
Beispiel sei angeführt, daß 


ete) j | „sin’r | - 
v(r) -ete’r + 5; = - (65). 
| r sın’r pr r’ 


$ econst ==0 


Für » ı würde an die Stelle der Bedingung (66), die aus (57) gewonnen wurde, eine entsprechende aus (58 
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die Bedineune (66). (67) mit ı 3: \fOo)dr=Insinr+1In € erfüllt. Der Fall 
Ad dl. 
(ir) a, 
) 


®\, anzeereeben. (64) zeeht dann ın die sorenannte 


l,ommelsche Gleichung über mit der Lösung y(r) = r"ı Zu (b, rs). Es bedeuten die a; vor- 


ist in dem Buch von Gray und Mathews 


vegebene, die 5b, geeignet zu bestimmende Konstanten. (Voraussetzg. a,-a, 0). In der 
Fußnote'”) ist Literatur zu VII angegeben. 


de 


VIll. Rechnungsgang und Fehlerbestimmung bei den Zahlentafeln. 4 (.) und F brauchen 
U x 


weeen (27). 28) nur für 20 berechnet zu werden, wenn man sieh auf reelle Argumente 
beschränken will. Ebenso eünstie steht es beı den Interralen der Besselschen Funktionen, 
denn es ist 


JtHd)AdAE+H\J(kHD)dE=OV | | (69). 


ı . ı) 


JI i ; PRREN . . r : : . 
Setzt man? „ I,” Ga)=K,(#), so folgt aus der "Theorie dieser Funktion für = >0 die 


Beziehung ?”) 


Daher eilt für #0 
SJ)AEHIK,(S)AE in\J,lieE)dE . | (70). 


Alle Tabellen benötiet man also nur für = 0. 


Die Intererale der Besselschen Funktionen wurden im wesentlichen nach der Methode 
der Differenzenreehnune eewonnen. Hierbei wurde die Stirlınesche Formel 


Vo hı 
P > N 5 
h A 


) [2 | - n : | R 


dr 2 - — E 
\yyı ‘ [9% 6 ISO I51> nn (‘ 


XV) hı 


benutzt. Die benötigten Funktionswerte wurden den siebenstelligen Tafeln des Watsonschen 
Werkes entnommen. Die Quadraturen wurden über Teilintervalle von 2h=0,1 bzw. 0,2 berechnet 
und danach alle Teilintegrale addiert. Fast stets war das Glied der Integrationsformel mit 1°, 
manchmal bereits das mit A,’ vernachlässiebar klein. Es wurde zunächst mit der Watsonschen 


0 
Stellenzahl gerechnet; als Ergebnis der Fehlerabschätzung ergab sich nachher die Notwendieg- 
2 ’ 
keit, einige Stellen abzustreiehen. Uber die Funktion \J,(T&)dE ıst nur zu bemerken, daß sie 
ı 


für -=6,8, 10 und 16 dureh Vergleich mit den- sich aus der gliedweise integrierten Reihe 
von J,(i.x) ergebenden Werten kontrolliert wurde. 


. x . r JA e . . . . .. vs . 
Die Funktion © „ I,” (ix) hat bei @2=0 eine logarithmische Singularität. Hier war 


also die Differenzenreehnune in der Umeebune von #0 nieht anwendbar. Für #04 bis 
x —=06 wurde die Reihe 


a i 4 te) ” 
ne J,(tx)ln. | \ —— I > ee el Fan 


a yıX 7 ARE -) 


j 


15) (ray and Mathews: A Treatise on Bessel Funetions. 


19), Jahnke-Emde: Funktionentafeln, (1909), S. 166 bis 168, (1933), S. 214 bis 215. Gray and Mathews: 
A Treatise on Bessel Funetions and their applieation to physies, (1895), S. 226 bis 240. ,ommel: Studien über die 
Besselschen Funktionen, (1868), S. 98 bis 120. Mathemat. Annalen 3, (1871), S. 475 bis 487: 14, (1879), S. >10 bis 536. 
Nielsen: Hdbeh. der Theorie der Cylinderfunktionen, (1904), S. 129 bis 149. Watson: Theory of Bessel Funetions, 
S. 95 bis 110. Weyrieh: Die Zylinderfunktionen und ihre Anwendungen, (1937), S. 86 bis 8». Schafheitlin: 
Besselsche Funktionen, (1908). Mac Lachlan: Bessel Funetions for Engineers, Oxford 1934. Partielle Differential 


gleichungen und Besselsche Funktionen: Kepinski: Mathemat. Annalen 61, (1906), S. 397 bis 4. 
20) Vel. z. B. Watson,S. 75, Gl. (5) a. a. O. 
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nach Potenzen von .r entwickelt und danach gliedweise integriert, woraus sich ergibt: 





17 





















Hierbei ist r’ 











Im Intervall #=06 bis 16 wurde wieder nach der Stirlineschen Formel (71) zerechnet. 
Man erhält so die Teilintegrale über die Intervalle ?2h recht zenau: indessen kann sich bei 
der Addition der einzelnen Teilergebnisse natürlich der Fehler summieren. Deshalb wurde 
das Intereral 


\ 
GE TEEN . 
\i-— H. BE, ee 
& . tl 
für ır IVO u. .r I6 noch auf anderem Wege berechnet. So ereaben sieh in der Reihe der 


Summation der Teilergebnisse zwei feste Kontrollwerte. Nach Gl. (39) ist nämlich 
\ö- YAuttir \e u > (1D) 
.) -) 0 - - kn . . . . . . . . . ° 


Dieser Ausdruck ıst für große Werte von . zur Berechnung geeignet. 


\e "Fdge=\+\. | . nr. 176). 


Für kleine = und eroße .r eilt näherunesweise 


\e sds=e | 5. Plel 5 En. 


Hierbei bedeutet P das Gaussche Fehlerinteeral. 





"ür große .r ist es möglich, ein e zwischen O und —- so zu wählen, daß (77) noch gilt und (75) 


sehon vernachlässiebar. klein wird. Als brauchbarer Wert ereab sich bei x = 16 





10 
13° 








£e = ATCCOS 













Auf diese Weise wurde (7D) für «= 16 bereehnet. Da (75) für «= 16 schon sehr klein ist, bedeutet 
selbst eine Fehlermöglichkeit von etwa 4°/, dennoch einen sehr genauen Wert für (74). Der 
Fehler fällt in Dezimalstellen, die später unterdrückt wurden. Das Intervall » I) bis 2 = 16 


wurde nun durch fortgesetzte Subtraktion der Teilintegrale von dem so berechneten Wert 
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erhalten. Bei x — 10 entstand die Möglichkeit einer Kontrolle durch Vergleich mit demjenigen 


Wert. der dureh Addition der Teilintegrale der von 0,6 aufsteigenden Reihe enstanden war. 
Der Unterschied der beiden Werte fiel außerhalb der als zeltend vermuteten Stellenzahl. 


Der Wert für £—= 10 wurde noch auf andere Weise kontrolliert, wie gleieh erläutert wird. 
, dq rn i "m 
Die Funktionen y (x) und sind nach (3) und (4) berechnet. Wegen (41) führte das 
{ ET, 
Auftreten der Funktion (75) in diesen Gleichungen dazu, daß die geltende Stellenzahl beı 
dq : i a , BZ j 
ya) und kleiner und kleiner wurde. Zur Kontrolle wurde (x) für Id nach Gl]. (Dh) 
Yen 
durch graphische Auswertung des Integrals mit dem Integraphen ermittelt. Es ergab sich, 
daß für I0O sorar noch die in der Zahlentafel als unsicher angegebene dritte Stelle richtig 
ist. Ferner war dieses Ereebnis auch eine Stützune für die Riehtiekeit der übrigen zur 
Reehnunze benutzten Funktionswerte an der Stelle Io, 


An Hand der Tafel für y (x) wurde sodann interpoliert mittels 


dq 
dx 


Iy, - 40 | " rI"y .. 
Jı Y (r,) Yo „ 0.1666 j — « (79). 


« 
w 
- 


Beträgt die Unsicherheit in der letzten angegebenen °') Stelle bei den y-Werten eine Einheit, 
so kann sie dureh die Interpolation bei hy’ auf vier Einheiten anwachsen. Der Vergleich 
dag 


mit der Tafel für gr 


ergab einen Fehler, der nie mehr als drei Einheiten betrug *”). 

Die Fehlerabschätzunge wurde bei den Besselintegralen in der Weise durchgeführt, daß 
der Fehler e,, jeder Teilintegration nach (71) bestimmt und angenommen wurde, daß sich 
alle Fehler im gleiehen Sinne addieren, was sicher nieht der Fall ist. Bedeutet 10” die 
aus den Watsonschen Tafeln zu erhaltende Genauigkeit, so wurde für &,, angesetzt: 


" | 16 | 
e,,=2hlW 1 ar an ie BR: 2 2 % » 0. MO 
) ( 


Für # wurde der maximale Wert des ersten vernachlässigten Gliedes der Stirlingschen 
Formel in dem betrachteten Intervall eingesetzt. Es ergab sıch stets 


ey; <pPhl0 >h RE En 


Wesentlich für den Verlauf der Fehlerabschätzunz war, daß sieh » im Intervall # 0 bis 16 
erheblich ändert. 


Nachdem die Funktionswerte entsprechend der Fehlerrechnung abgerundet waren, wurden 
sie dureh Bildung der Differenzenstaffel kontrolliert. 

aerr) v > but .. ® \ " . on dq . 

Für die Besselintegrale beträgt die Fehlerschranke eine, für g (@), j,. Im Intervall = 0 

Ü x 

bis 6 eine halbe, zwischen x = 6 und S eine Einheit der letzten angezebenen Stelle. Im Inter- 

ch ur do \. ce h 
vall S bis 16 ist bei g und j„ die letzte Stelle mit noch etwas größerer Unsicherheit behaftet 

2 


A 


und daher klein gedruckt. Eine überstrichene 5 ist durch Erhöhung beim Abrunden ent- 
standen und bleibt bei weiterer Abrundung fort. 


Für die Interpolation benutze man bei Vernachlässigung der dritten bzw. vierten Diffe- 
renzen die Formeln von Bessel bzw. l,aplace-Everett; man benötigt dann nur die 
zweiten Differenzen. S61 


21) Die klein gedruckten Zahlen in Tafel 3 gelten nur dann als „angegeben“, wenn weitere Ziffern nieht von 


R d 4 R U 
handen sind. | für x = 15—16.] 
Ada 
92 oT nis dy 
>2, Allerdings war es im Intervall a 6,2 bis a Ss nieht möglich, dureh Interpolation genau genug zu er 
( 


halten, was nach dem Verlauf der Differenzen auch verständlich ist. 








«e) 
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\ A‘ 
Ji d: 1? r AGD: 1? 
0 0 

0.0 0) 5.2 37,5146 1.1736 
0.1 V.LOOOSZ 0.000502 2.4 44.9541 1.4114 
(> V.IOOHHN 0.001005 2.6 33.204995 1.6979 
03 0.30? 0.001517 I. 63.7538 > 0438 
0.4 0.105369 0,002043 6,0 76.0464 2.4609 
5 0.510515 V.UOO25S0 6,2 30.7999 > I640 
0.6 0.618245 0.003139 6.4 108.5174 3.9715 
0,7 0,729114 0.003721 6,6 129,064 1.305 
ILS ().S843704 0.000433 6.5 155.400 2.190 
00 0.062625 0.004975 7.0 156.184 6,261 
1.0 1.086521 0,0065 1.2 223.229 1.51 
1.1 1.2 1607 V.OOHBS 1.4 657.825 9.11» 
1.2 1.35200 0.00715 7.6 321.536 11.001 
1.3 1.149508 000798 1.8 386.248 15.284 
1.4 1.641614 0.0DOSST Su 464,244 16.042 
1.» 1.0607 0.0098? 8,2 558.282 19.37 
1,6 1.937582 0.01085 8,4 671.701 23.416 
1. >.15642 0.011098 S.6 SOS.536 28.300 
1.S > 34900 0.01317 8.8 973.671 4.210 
1.0 > 55475 0.01450 930 1173,016 41.355 
() > 77500 001592 92 1113.72 50.03 
a 3.01117 0.01746 9.4 1704.45 650.49 
2.3 3.206480 0.011915 36 2059.67 73.19 
>23 HRaHN EN. 002100 08 2480,08 8855 
2.4 3.831536 0.022099 10.0 293.045 107.18 
In 1.148513 0.02519 10.2 3613.175 129.70 
I 1.419009 0.02757 10.4 1363.02 157.02 
7 1.8596? 0.0301 10.6 5269.88 190.11 
2,8 3.250933 0.035303 10,5 6366.85 >3),22 
0) 5.659207 0.053615 11.0 7694.04 278.80 
0 6.165096 0,03956 11.2 9300.03 3377 
3.1 6.669415 0.04330 11.4 11243.7 409.2 
32 7.22116 0.047357 11.6 135096.6 195.6 
3.3 T.S?0>8 0.0184 11.8 16449. 1 600.6 
3.4 8.47124 0.05675 12,0 19804,2 121.9 
1.) 3.178095 0.062209 12.2 24071,2 882.5 
6 494875 VUHTOS 12.4 29130,5 1069.4 
1.1 10,786 007441 12.6 35259,.2 12964 
N 11.659872 008146 12,8 42684,3 1972.0 
0 12.609237 0,08919 13.0 »1681.4 1906.0 
[.0 3.779521 009766 13.2 62584.9 2311.4 
t.1 1195571 0,10695 13,4 757990 2803.3 
1.2 16.245316 0,11714 13.6 91816,5 3400 
1.3 17.647755 0.1283? 13.8 111235 4125 
1.4 19.180655 0.140595 14.0 134778 3003 
1.9 20,85412 0,15400 14,2 163924 6072 
1.6 22.068158 0,1687 14,4 197942 7366 
1.7 24HTTS 0,1849 14,6 239926 8939 
1. 6,8580 0.2027 11.8 IIUN4I 10849 
1.0 29.2427 0,2222 15.0 392621 13165 
U) 31.8487 0.2435 15.2 427598 15981 
5.1 34.609815 0,2669 15.4 518476 19399 

15.6 628795 23548 

15.8 762658 98590 
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nr x 
U ä nm 22 
1 io Ho" (iä)ds Pr x i— Bd 
5 : 0 N S 0 . 
0.00 (0 >.8 1.569836 
0,02 0.100563 6.0 1.569695 
0,04 0.173417 6.2 1.6986 
6.4 1.57004 
01 0,3412165 6.6 1,570184 
0.2 0.547116 0.050103 6.8 sul 
0.3 0.701964 0.031219 7.0 396 
(0.4 0.825993 0.022119 2 472 
0,D 0.927103 0.016707 7.4 534 
0,6 1.011906 0.013111 1.6 D84 
0.4 1.083508 0.010556 7.8 24 
0.8 1.144734 0.008653 8.0 657 
0.0 1,197217 0,007190 8.2 683 
1.0 1,242510 0.006039 8.4 705 
1.1 1.281764 0.005111 sh 122 
2 1.315907 0.004357 8.8 736 
1.3 1.345693 0.003734 90 747 
1.4 1.371745 0.003215 9» 757 
1,5 1.394582 0.O0O2T8O 9.4 64 
1.6 1.414639 0.002411 96 770 
1. 1.432285 0.002096 8 775 
1.5 1.147835 V.001831 10.0 1.570779 
1.9 1.161554 V.0UO1598 10.2 1.5707825 
2,0 1,173675 0.001406 10.4 85 
2,1 1.48439 0.00122 106 87 
22 1.149389 0.001410 10.8 Sg 
2.3 1.350229 000094 11.0 g) 
2,4 1.350975 000084 11. 915 
25 1.516537 - 0,00074 11.4 924 
11.6 1.970793 1 
>26 1.52225 — 0,00262 - 
11.8 dA 
Ar 1.553213 0, 206 
12.0 12 
3.0 1.539095 0). 162 . 
12,2 16 
3.2 1.546159 0, 126 
12.4 19 
3.4 1.55109 -0. 102 Be 
e - . 12,6 2 
3.6 1.359015 000078 k 
a 12,8 4 
3.8 1.355815 0) 64 BE 
rn n 13.0 3b 
1.0 1.56065 0, ) Ze 
ERER 13,2 zZ 
1,2 1.562695 0, 40 ” 
“nimm 13.4 ID 
1.4 1.564295 0, 24 8 
FREE 2 13.6 1.070959 
4.6 1.536553 0, 2 
S 1.356656 0) 2) ER 
h = ® n 14.0 1.>7079606 
).() 1.36739 —(, 17 a zo 
en n . R )l+) 
52 1.56805 () 13 
2 an 14,8 621 
).4 1.356858 0), 10 R ri 
Mr 15.2 625 
AR 1.536901 V.VBODOS ö rd 
15.6 627 
16.0 1.57079629 
X 
im Tu,» ade = -T = 1,570796827 
ım \’ m. 0 \ts)as > IIVNIIDDS A 





xD % 


I. 


0 
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dl q 





. d Y 
ya / er FA LE Ir 
1.5708 00 1.0000 0.147 70 0.0221 
1.4746 01 0.9247 0.142 1.2 0.0208 
1.3856 0.2 0.8955 0.138 1.4 0,0195 
1.230932 03 0.7919 0.1246 7.6 0.0185 
1,2271 0.4 0,7334 0.151 1,8 0.0174 
11565 0.5 0.6796 0.128 S0 0.0165 
1.0910 0.6 0.6301 0,124, 8.2 0.015, 
1.0303 0.7 0.5846 0.12, 8.4 1.014, 
(1.9740 DS 0.5427 0,11, 8.6 0,014, 
0.0217 0.9 0.5040 0,11, 8,0 0,015, 
0.11, 9.0 0.012, 
0.8731 1.0 0,4685 0.11, 2 0,012, 
0.7859 1.2 0.4054 0.10 9.4 0,011- 
0.7108 1.4 3918 0.10, 6 "011, 
0.6446 1.6 0.3060 0.10, O8 0.010, 
Unze 1,9 0,2670 0.10 10.0 0,010, 
0,5375 >00 0,2336 0.099, 102 0,0099. 
0.4037 >. 0.2050 0.097, 10.4 0,0OO9n- 
0,1552 2,4 0,1805 0.095, 10.6 0.0092, 
0.41 I; 0.1904 0.098, 10.8 V,OOSS, 
0.391? IS 0.1412 0.091. 11.0 0,VOSD,, 
5646 0 0,1204 0.090, 11,2 0.0082, 
0.3409 > 0,1118 0,088, 11.4 0.0079, 
0.3197 3.1 0,1000 0.086, 11.6 0.0076. 
V.ZOOS 3.6 0.0897 0.085, 11.8 0.0073, 
VL2SBS 3.0 VV8US LOS, 12.0 0,0071. 
0.2684 1.0) 0.0729 0.082. 1? > 0.0069, 
01.2545 1,2 0661 0.OS1. 1?.4 0.0067, 
0,2419 1. 0,0000 0.079, 12.6 V.OOGA, 
0.2305 1.6 0.0947 D.OTS, 12.8 0.006, 
0.2200 1.8 V,USOO 0.077, 13.0 0.006, 
0.2104 an 0.0459 0.076, 13.2 VOODS 
V>OI6 2 0,0421 0.075, 13.4 0,005, 
0.1936 7.4 0.0380 0.074, 13.6 V.00D, 
0.186] 1. 0.0350 
0.170> 5.8 0.0333 0.072, 14.0 0,005, 


0,1727 
1668 
0,161 
0.196 


0,151 







0.0309 





V,U2SN 
V.O>PHN 


AUMZETE 








0.070, 
0,06>, 
0,06- 
0,06. 


0.06, 
















0.004. 





0,004, 


0,00, 






0.00, 





0,00, 
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KLEINE MITTEILUNGEN 


Allgemeines über Leitkurven in Rich- 
tungsfeldern. \ls Hilfsmittel für die graphische 
Integration einer gewöhnlichen Differentialgleichung 
erster Ordnung y'’=f(r,y) kann man bekanntlich 
in geeieneten Fällen sogenannte Strahlkurven 
verwenden. Unter bestimmten Umständen arten 
diese in Strahlpunkte aus und führen dann zu 
dem Begriff der Leitkurve. Wir knüpfen im 
foleenden an die Darstellunz an. die Willers 
in seinem Buche über .„Methoden der praktischen 
Analysis” (Berlin 1928, S. 299 ff.) diesem Gegen- 
stand geeeben hat. Für eine beliebige Ausgangs- 
kurve g(r,y)=c erhält er die Gleichung der 


zugehörigen Strahlkurve durch Elimination von 
r und 4 aus 
(ut rVx 
:(a.J)=ı 
Ix vu IuV 
!pyT 9x (4 


n(.r.%) y--fJ- 
Ix Fu 


| 
| 
fy x | 
| 


Man sieht. daß bei Gültigkeit von (1) die Identität 


SE | ; 


besteht. die im übrigen die charakteristische Eigen- 
schaft der Strahlkurve unmittelbar zum Ausdruck 
brinet. Willers weist dann darauf hin, dab 
sich die Strahlkurve auf einen Strahlpunkt 
zusammenzieht. falls die Differentialeleichung die 
Form 


. Y 7\g .# Y)) Pr 
1 i i ru 
I hig (rs. y)) 
hat. d. h. falls 
E=hlu), 9 giu) mit u P(t,Y4) . .(8,1) 


ist. Wir bemerken zunächst. daß die Eigenschaft. 
daß die zugehörige Strahlkurve sich auf einen 
Strahlpunkt zusammenzieht, nicht nur der einen 
ausgewählten Kurve zukommt, sondern auch jeder 
anderen Kurve der Schar g(r, y)=const. Denn 
= und n hängen in (1) nieht von 9 unmittelbar. 


ab: das ist aber für alle Kur- 


Auf dliese 


sondern nur von 
14 


ven der Sehar g (.r, 4) = eonst dasselbe. 


Weise erhalten wir für jeden Wert des Vara: 
meters a einerseits eine Kurve 4, der Sehar 
y (r. y) = eonst, andererseits einen Punkt II (£, 7): 


IE (u). n(W)=Il,. Die Punkte II, erfüllen eine 
Kurve X. nämlich die durch die Parameter- 
darstellune (3.1) gegebene. Jeder Punkt Il, ist 
Ausstrahlungspunkt für die Linienelemente in allen 
den Punkten P,, des Ricehtungsfeldes, die auf der 
betreffenden Kurve N, liegen: er ist damit auch 
Ausstrahlungspunkt für die Tangenten an die 
Interralkurven in diesen Punkten P,„. Im Anschluß 
an eine frühere eigene Arbeit (..Ein abbildungs- 
eeometrisches Verfahren zur Darstellung von 
Richtungsfeldern ... .” in „Deutsche Mathematik" 


1938. S. 353ft.) und eine dort angegebene Ver- 
öffentliehunge von Alferov nennen wir Y in 


sinneemäßer Verallgemeinerung des bei linearen 
Differentialgleiehungen erster Ordnung bekannten 
eleichnamigen Begriffes eine Leitkurve des 
Riehtungsfeldes. die Kurven N, die zugehörigen 
Isopunktalen und die Geraden P, II, die 
Leitstrahlen. 


Die in (3) und (3,1) angeführten Bedinzungeen 
für die Ausartunge der Strahlkurven in Strahl- 
punkte sind völlig gleichbedeutend mit der einen 
von Alferov anzerebenen. dab für die Existenz 
einer Leitkurve das identische Verschwinden der 
N l(E 


( m 


Funktionaldeterminante 


hinreichend ist. Aus dieser Formulierune aber 
scheint ebenso wie aus der Bemerkung von Willers. 
dab man „Typen von Differentialgleichungen auf 
stellen kann, für die die Zusammenfassung der 
Richtungen auf bestimmten Kurven Strahlpunkte 
ereibt”, hervorzuzehen. daß bisher die Gültiekeit 
des folgenden Satzes nicht bemerkt worden ist. 

Im Riehtungsfeld jeder Differential 
gleichung y’=f(r,y) läßt sich jede be 
liebige Kurve als Leitkurve wählen und 
die zugehörige Isopunktalenschar un 
mittelbar anzeben. 

Der Beweis dieser Tatsache ist denkbar einfach. 
Sei 


= notwendig und 


O(r,Y) 


& hin). 7 gu) - .. 0,14 


eine Parameterdarstellune der willkürlich als 
lseitkurve gewählten Kurve. dann muß nach (2) 
die Identität 
HN yıu . 
/ /f) ) 
£ h (wu) J 
bestehen. Das ist aber dann und nur dann der 


Fall. wenn wir 
EZ, Y 2 hie)f(lr, y Er h) 


setzen. Damit haben wir 
Schareleiehung der 
als Parameter zewonnen. Man erkennt «das 
entweder unmittelbar oder foleendermaßen: Sei 
.4) die Auflösung von (6) nach vv, so lautet 
die Differentialzeleichunz der Schar (6 


jedoch bereits «die 
Isopunktalen mit % 


Hl 4 Mn 


einerseits 
Se u >. 
4 „ 


aus (6) dureh totale 
x (bei festem «). Nas liefert 


andererseits erhalten wir sie 
Differentiation nach 
zunächst 

() y' hu) /« 7 In *® Yy') a 


und nach Auflösung nach A()-und Kinsetzen in 


(6) die beiden Gleichungen 


I—y 
I«tI1y° Y | 


gW)=y+l- er 


Jı (MH be 77 


- 
2 


die zusammen die Differentialeleiehung «der Schar 
(6) darstellen. (7.1) und (7.3) ergeben vereint 


you 
hin) z + ey ER 
lx 14 lu 1a ”. 
ar . . 
4 vTrrox | 


gu)=y—l- 


IxIy Iy Px 


Die rechten Seiten dieser Gleichungen stimmen mit 
denen von (1) überein. der Zusammenhang zwischen 


den linken Seiten wird dureh (4) vermittelt. Das 
u oO (5, N) 
heißt aber: da weren (4) ohne weiteres 
(1.4) 


erfüllt ist. arten die Strahlkurven der Kurven- 








Mitt 


'Ine 







also In 


der Tat die Isopunktalense ar zu der dureh (4 

Ikürlieh vorzerebenen Leitkurve dar. 

Die damit bewiesene Tatsache, daß es im Rich 
tunesfeld jeder Differentialgleichung / ir, Yy 
l,eitkurven gibt, wird für die P’raxis der graphischen 
Interration immer dann Bedeutunz haben, wenn 
die Isopunktalen einfach genug zu zeichnen sind. 
\lan wird sı var sacTen können. dab das Lie itkurven- 
verfahren infolze seiner besonderen Eigenschaften 
wuch dann d« Vorzue eerenüber dem Isoklinen 
verfahren vertedien: wird, wenn die Schar deı 
Is punk en nieht wesentlich schwerer herzustellen 
st als die der Isoklinen 

Zur weiteren Orientierune im Richtunesfeld kann 
dann noch der folgende Satz dienen: 


l,aeitstrahl 7, 
liert der 


ie 


Il,, «die 
Berührungspunkt / 
Inteeralkurven. d.| 


Berührt zugehörige 


sopunkt le 


| 
tufıdder Wi 


\t l -() z 


uınktskurve « \.h. 


che ' 


er ist selbst Wendepunkt der durch ihn hindureh 
eehenden Interralkurve. 
Im Berührunespunkt muß nämlieh die Riehtung 




















der Isopunktale mit der Riehtung des l,eitstrahles 
übereinstimmen. Wir müssen also in (7.2 y / 
setze Das führt zu der Beziehung 
f. , () 
Schließen wir die Stellen aus. in denen 2 — hin 0) 
und ıinfoleedessen weeen (DD) / 7 nieht mehr 
beschränkt ist. so erhalten wıı 
jj () Sn 
as ist aber erade «die Gleichune der Wende- 
pri kt kuı t 
a 3 
| / 
k / 
| { 
y 
8 v 
Age 
Bild 1 
Wählt man als Leitkurve ireendeine Gerade. 
lann il Shi die Auflösung der Isopunktalen 
eiehune ( nach dem Parameter # besonders 
nfach durehführen. Als Beispiel betrachten wiı 
ie Rieeatische Differentialeleichune 
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\ls Leitkurve nehmen wir die „-Achse und als 
Parameter zunächst n selbst. Dann ist 
hiu (), gu u 
und «die Isopunktalenschar lautet nach (6 
| ö 
// 27 1,.t 7 g“ 7 u. 
Pr 
iR 2: / C 
Sie besteht demnach aus Kreisen um IV, a als 


Mittelpunkt. Es erscheint zweckmäßig, an Stelle 


von» « den Halbmesser dieser Kreise als Para 

meter ? zu wählen. Dann erhalten wir auf der 

Achse als Leitkurve eine nach ? zu beziffernde 
e l- 

Funktionsskala für 7 Das fertige 
a 

System Leitkurve - Isopunktalen ist in Bild 1 


Wert ae 4 Kine dureh 

Inteeration Inteerralkurve 

ist mit einigen Tangenten und der Wendepunkts 
> 2y 

he- 


-) 


dlarzestellt. 


(5) 


für den 


eraphische 


( 


eWONNENE 


( 
SAL 
ui 


kurve die die Gleiehun®e .r? y° 
c+?2y 


Die Et schlossene L,ösun®e ler 
Differentialeleichune 


sitzt. einzezeiehnet. 
R jiceeat 


eerebenen ıschen 


ereibt sieh übrieens in der Form 
Zi | | 
7 J 
f j l 
Ba 
- 


Il die 


r Z, und Z, der Anfangsbedingung angepabte 
/ylinderfunktionen sind 
Breslau. Helmut 67 


Heinrich. 


Resonänzschwingungen von quadra- 
tisch gedämpften Systemen. \chwinzungen 
t quadratischer Dämpfung sind deshalb der 
mathematischen Behandlung schwer zueänglich. 


3 
ill 


weil man die zesamte Schwinezeune nieht dureh eine 
einziee Differentialeleicehune beschreiben kann. Da 
die Reibuneskraft stets der Eewerunze enteeeren 


verichtet ist. Jautet nämlich die Gleiehunge entweder 


Tr N .r: or Pi für () la 
oder 

Y sr? o?ı3 Pi für # () Ib). 
wenn !) die der Schwineung ein prägte Kraft 
ist. Die Lösung des Systems (1) ist im allgemeinen 
überaus sehwierie. Es gibt aber ähnlich wie bei 


«h 


Schwineuneen n 


von 


it linearer Dämpfung eine perio- 
disehe Kraft eanz bestimmtem Verlauf. die das 
System zu rein sinusförmigen Resonanz 
schwingune anreet und diese Schwineunge stationär 
aufreeht erhält, weil in jedem Augenblick 
rade die Reibungsverluste deckt. 

Den Verlaul Kraft kann man leicht 
stimmen, wenn man die Gl. (1) in zwei Teile spaltet: 


seiner 
S]e 


ge- 


dieser he 


T wo“ Tr 0) 124 
Sr Pit 3 
und eine Lösung von 
IB 1sın wf | 
In (3) einsetzt: 
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jild 1. 
l,äßt man diese Kraft, die die Kreisfrequenz 2% 
i ön . 
hat, jeweils nach einer vollen Periode ihr 
() (1) 


Vorzeichen wechseln, so kann das quadratisch ge 
dämpfte System mit seiner Resonanzfrequenz rein 
sinusförmie schwingen, da sich dann seine Be 
werungsrichtung nach einer halben Periode, d. h. 
ändert. 


ebenfalls nach der Zeit Die Amplitude 


(1) 
des schwingenden Systems ereibt sich aus (5) zu 
l Pi 
’ (ı) N) 


Bild 1. 


Kraft- und Bewerungsverlauf zeıgt 
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Zur näherungsweisen Lösung zweier 


Gleichungen mit zwei Unbekannten. <o|| 
man zwei Gleichungen mit zwei Unbekannten 
f(r.y 0 und g(r.y Bez Wi, 1) 
lösen. so muß man sich zunächst durch rohes 
Skizzieren der Kurven 
f(r,y 0(6,) und gir.Y )(E,) 


die ungefähre Lage ihrer Schnittpunkte überlegen. 
Seien .7,.4, die so gefundenen Koordinaten für 
die erste Näherunz eines gewissen Kurvensehnitt- 
punktes. Zur Verbesserung der Näherung wird 
meist das Newtonsche Verfahren oder das Ite- 
rationsverfahren vorgeschlagen. Beide Verfahren 
sind aber in vielen Fällen entweder überhaupt nicht 
konvergent oder ihre Konvergenz zegen den ge- 
suchten Schnittpunkt ist nur ziemlich mühselig 
sicherzustellen. Es ist daher von Interesse, ein 
Verfahren zu entwickeln, das in den meisten prak 
tisch vorkommenden Fällen sicher zum Ziele führt. 
und zwar immer dann, wenn der Sehnittwinkel der 
beiden Kurven &, und &, (wozu also / und g als 
mit stetigen Ableitunzen erster Ordnung begabt 
vorauszusetzen sind) von 0 bzw. zz verschieden 
und die Auszangsnäherunz schon (in einem zleich 
zu präzisierenden Sinn) genügend genau ist. Das 
vorzutragende Verfahren hat den Vorteil. daß sämt- 
liche Näherungspunkte im gleichen von &, und GC, 
bestimmten Qua«dranten liegen und die Näherung 
monoton verläuft in dem Sinne, daß die Entfernung 
der einzelnen Näherungspunkte vom wahren Kurven- 
schnittpunkt monoton nach 0 abnimmt'). Die er- 


I) So ist also dieses Verfahren als eine Übertragung des 
Verfahrens zur Lösung einer Gleichung mit einer Un- 
bekannten anzusehen, wie es an erster Stelle in dem zu 
sammenfassenden Berieht geschildert wird von R v.Mises 
und H. Pollaezek-Geiringer, Praktische Verfahren 
der Gleiehungsauflösung, ZAMM 9 (1929), S. 58 bis 77, ins 
besondere S. 59 bis 61. 


Kleine Mitteilungen 4 


.- 


forderliche Rechenarbeit ist dabei kaum größer als 
die für das Newtonsche Verfahren aufzuwendende. 
wenn man das Verfahren in der hier gewählten 
Kombination zwischen Rechnung und Zeiehnung 
durehtührt. 

Ich berinne mit der Schilderung des Verfahrens: 
Durch Skizzieren der Kurven &, und &, sei ge- 
funden. daß der wahre Kurvensehnittpunkt in 
einem die erste Näherung z,. y, enthaltenden Recht- 
eck R liezen muß. Wir denken uns zunächst an 


die Fläche /(r,.y) in jedem Punkt über AR die 
Tangentialebene gelegt: 
Pi) N 
7 /(r.y) / N Tr) 0/ ( Y). 
O4 O4 


Unter der Fallinie einer Ebene ar +by-+ cz-+d- 
verstehen wir diejenige (nur bis auf Parallelver 
schiebunz bestimmte) Gerade der Ebene, die mit 
der r. y-Ebene den größten Winkel g einschließt. 
Iirsiehtlich ist 


CO> Y 


Suchen wir jetzt unter unseren Tangentialebenen 
diejenige steilster Fallinie, so bekommen wir ihren 
Winkel 9, gzeren die r,4-KEbene aus der Formel 


l 
Bi 2 O/ 2 


+ 


COS 14 


Max | 


Jetzt legen wir durch den Raumpunkt 7. 41.7 (r,, Yı) 
den .Fallkegel”. d. h. den Kreiskegel, dessen Er 
zeurende mit der x, Y-Ebene diesen Winkel y, bilden. 
Sein Sehnitt mit der .r, y-Ebene sei der .„Fallkreis” 
\t,. Analog zeichnen wir (beide Male in entsprechend 
vergröbertem Maßstab) für die Fläche g(r,.y) den 
Fallkreis 8, (Bild 1) 

Die Tangenten aller Niveaukurven f(r. %) 
erhält man mit Hilfe der Formel 


eonst 


O/ 


%) n 


O/ 

O4 

Mit Hilfe dieser und der entsprechenden Formel 
für g(z.u) (oder dureh Skizzieren) sei noch ge 
funden, dab in AR die Tangenten an die Kurven 
/=0 bzw. 970 um zwei mittlere Richtungen 
f, bzw. ft, höchstens um den Winkel a, bzw. «, 
schwanken. #, und f, möwzen (gemessen in den 
7.4, enthaltenden Quadranten) den Winkel «a ein 
schließen und es sei 


Wir zeichnen an die beiden Fallkreise je die beiden 
unzünstigsten Tangenten f,'. ft,” bzw. t,’,t,”. Durch 
Feststellung der Voizeichen der partiellen Ab 
leitungzen läßt sich leicht entscheiden, nach welchen 


Seiten hin die Flächenhöhen /f und 9 abnehmen. 
Bild 1 ist gezeichnet unter der Annahme, dab 
6, unter St, verläuft, &, rechts von N, Durch 


t,” und #, wird dann ein Winkelraum abgegrenzt, 
in dem der wahre Kurvenschnittpunkt liegen mub. 
Denn die wahre Kurve f==0 liert unter dem von 
!,,t,” und dem Zwischenbogen von st, gebildeten 
Kurvenzug, analog 9=0 rechts von dem dureh 
!,.t," und dem Zwischenboge ı von St, gebildeten 
l,inienzug. Nehmen wir den Schnittpunkt von 
t,” und #, als neue Näherung 7. Y,. so liegt diese 
Näherung sicher diehter an dem wahren Kurven- 
schnittpunkt P, da das Dreieck mit den Ecken 
P, #,.4, und 7,.%, bei der letzten Ecke stumpf- 
winklig ist. Dazu muß man noch bemerken, dab 








In Kleine Mitteilunzeen 


Schluß nieht mehr notwendige richtige ist. 
und /,° sieh nieht schneiden (Bild 
Nimmt man jetzt etwa die neue Näherung auf 
Zwischenbogen 
weit links dieht an #, 


dieser 
falls t, 
«dem 
von Mt, an. so könnte P z. B. sehr 
» liegen. so dab die Ausgangs 
besser gewesen wäre als die neue. Im 
Bild hätte man den Berührungspunkt 
von £, als neue wählen können 


näherunge 


Fall vol 


Näherung 














x 
 - 
2 t; 
2, Ya N\ 

A881 1; 

Pz 

“2 
Bild I u ) 

\lle dureh Wiederholun®e (lieses Verfahrens vo 


wonnenen Näherungspunkte liegen im gleichen dureh 
6, und @, gebildeten Quadranten und rücken in 
diesem monoton auf den wahren Kurvenschnittpunkt 
u. Daß die Kontizuration von Bild 1 oder 2 vor 
liert. folgt aus der Voraussetzung (2). da #,” und #,’ 
nach der den Fallkreisen abgzewandten Seite den 
Winkel Q "TG, As T einsehließen. 

\ls Beispiel suchen wir eine Lösung der beiden 
Gleichungen 


j / ) / / Ale r | () 
und 
ia. y\ 2 2? + y° g” ’ 0, 
Um zunächst die Kurven /=0 und 9g=0 zu skiz 
jeren?). schreiben wir: 
() | | 532°’ ı->2r7r ® 
(/ () Pen 0 l 9 173. 22 . 
Wir reehnen : 4, (0 I. und mit Hilfe der unteren 
und oberen festen Skala des Rechenschiebers lesen 
wir aus 
Due 0)) 3 sofort 9, (0 1,86 


Weise rechnen wir mit dem 
Rechenschieber folgende Werte und tragen die 
‚uzehörieen Punkte in Bild 3 auf. indem wir die 
Punkte von /==0 durch Kreuze. die von g=V 
dureh Nreise unterscheiden: 


ab. In der eleiechen 


kleine 


Man könnte übrigens leicht y aus der ersten Gleichung 

ısrechnen und in die zweite einsetzer oder man könnte 

sieh wenigstens die zweite (Gleichung vert infachen. indem 

man den Wert von „3 aus der ersten einsetzt. Das ist 
| } , 


edoceh alles nieht nötige. 


Es 
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y,(l - 18 Yılzı 


( ( um 
(0.97 y.{(0,9) = 


> 


Die richtige Schnittpunktsabszisse des aus Bild 3 





0 x 





-) —— 
v alu 


zu entnehmenden Schnittpunktes liegt also zwischen 


0,9 und 1, die Ordinate zwischen 1.3 und 1.1. 
Wir bilden 

O/ Of \ 

I 2. of) 

0 Q 7 

N N 

( (] z () A “ 

m.2* A ARM » Yu 2y 
N \ } 
( I ( 4 


Als Ausgangspunkt der jetzt (auf Millimeterpapier) 
anzulegenden Zeichnung wählen wir etwa 











dr 0,05 Yı 1.2 
Az, 
12 
7205}, 
r 
ABB14 wi 
0.95 0,955 0,96 
zild 4. 
Dafür ist 
Fi2.%, 2.9712 — 1.728 19 —+1 0.056. 
glaıYı 1.728 — 1,715 1,44 — 0,903 + 3== 0,094. 


Da die Fallkegel außerordentlich steil sind. werden 
die Radien der Fallkreise jedenfalls nicht viel 
erößer sein als 0,01. Für die Tangentenrichtungen 


der Kurven /=tÜonst und 9 = Uonst findet man 
O/ 
i Ol 9.72 2 ‚,‚ °’-+-2r 
// m = R )ZW. Yo z . 
Of 322 ai 3? 2 y 
ÖY 
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4 ; 
Infolze A) und z— <U haben wiır die un 
Mr 


VO X ( 
eünstiesten Tangenten von Bild 1 nach rechts hin 
an die Fallkreise zu legen, und zwar überlegt man 
sich sofort, daß wir nur je die steilstmögliche zu 
zeichnen brauchen. Da auch diese Tanzenten sehr 
steil sind, ist anzunehmen, daß der wahre Kurven 
schnittpunkt eine Abszisse etwa zwischen 0,95 und 
0,96 besitzt, seine Ordinate wird zwischen 1.2 


und 1.21 lieren. Damit ist der Maßstab der 
Zeichnung festgelegt (Bild 4). Wir finden 
1 1 
CO> Y m ma) 
V1-+ (9-0,96? — 2? +9-121° 648 
1 
UO> 4 9 
y1+16-0,96?+2-0,96)?+(3-1,21?— 21,21)? 
1 
Are 


die Radien der Fallkreise werden 


ER M 0.056 = 
fr, Aue? 0.0073. 
1291 V 7,13? — 1 
0.094 
1), 0,0122. 
y 777 — 1 


Buchbesprechunzen 


59 


Weiter wird 


B 41.0,06? ) 146 
4 ‚4b 
. 3.1,2° 

6-0.067 2.0.06 
Y ) ) .)2 ) | ) 1,11 
) =. Z Z 


Aus der Zeichnung lesen wir ab: 


2. = 09814. (» 1.2038 . 
Kinige wesentliche ’unkte der 
überschreiten zwar das der Rechnung 
zugrunde gelerte Gebiet 0,95 a 0,96 und 
1.21 Y 1.2 um ein Geringes, das wird aber 
nichts Wesentliches ausmachen. Tatsächlich rechnet 
sich, erstmals unter Verwendung einer fünfstelligen 
Lorarithmentafel, 


Bemerkung: 
Zeichnung 


/ (Dos Ya 0.0014. 4 (5, Y») VOO32. 


Wir sind also trotz der eben erwähnten Ungenauig- 
keiten doch im gleichen Kurvenquadranten ge 
blieben. Man wird jetzt gegebenenfalls das Ver 
fahren fortsetzen, indem man von 2,.4, ausgeht. 


Braunschweig. Rudolf Ieliseh. 881 
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DINGLER, 0. Prof. a. d. 
München. Die Methode der 


Techn. 
Phy 


sik. 4228. München 1938, Verlae von Ernst Rein 
hardt. Preis geb. 13 M. 
Das Buch bildet einen zusammenfassenden Ab 


schluß der in einer Reihe früherer Veröffentliehun 
een durchgeführten Untersuchungen über die Frage 
nach #iner die Sicherheit und Alleemeineültiekeit 
physikalischer Aussagen gewährleistenden Methode. 
Dineler ist der Ansicht, dab es falsch sei wie es 
die heute herrschende „Schulmeinung" tue . die 
physikalischen Gesetze als allgemeine Eigenschaften 
der Natur zu betrachten und durch messende Ver 
suche auf Grund des Prinzips der Induktion er 
schließen zu wollen: damit komme man nieht zu 
einem eindeutireen. die Sieherheit der Aussagen ver 
bürrenden System, sondern zu einer Vielheit nicht 
miteinander im Einklang befindlicher Theorien. 
Ir verlangt wie für die Mathematik so auch für die 
P’hvsik einen axiomatischen Aufbau, der zum „ein 
deutie-methodischen System” der Physik führe. 
Die Methode. die «das Jeiste, müsse auf «dem 
Prinzip der einwandfreien und lückenlosen Begrün 
dungx beruhen (also unwiderleglieh sein). das Prinzip 
der Ordnung befolgen (d. h. jeder Schritt dürfe nur 


Schritte voraussetzen, die ihm vorangzehen) und da- 


mit die Forderung der Eindeutigekeit (oder, was da 
mit oft äquivalent sei, der Einfachstheit) er- 
füllen. Die Grundlagen dieses Systems werden von 
Aussagen „vorsvstematischer" Natur gebildet. die 
nach Dinelers Meinung unmittelbare Sicherheit in sieh 
tragen: dazu gehören vor allem die sog. Relations 
erkenntnisse (Gleichheits-. Ähnliehkeits- und Ver 
schiedenheitserkenntnisse). Zwischen (diese vor 
systematischen Aussagen (die im wesentlichen De 
tinitions- und Handlunges- oder Willenssätze sind 
und das auf ihnen dureh die Anwendung der Ge 
setze der Logik aufgebaute System schiebt sich 
aber nach Dingler noch der „Realisierungsprozeb“, 
der Bereich manueller Betätigung, der für die 
(Grundlerung der Wissenschaft ebenso erundlegen(d 
sei wie Experiment und Theorie was die heutige 
’hysik leider nicht berücksiehtige: sonst könnte 
sie z. B. nieht behaupten, die Frage der Krümmung 
des Raumes sei mit Hilfe von Messungen zu ent 
scheiden. da ja aller Instrumentenbau auf der An 
wendung der Euklidischen Geometrie beruhe, man 


Hilfe solcher Instrumente nie zu «einer 
andern Form der Geometrie gelangen könne als zur 
Kuklidischen. Bei Grundlerunge der Geometrie (die 
der der Physik vorauszurehen habe) kommt Ding 
ler vor allem durch Verwendung des Prinzips 
der Eindeutiekeit, das ihm die Definitionen für die 
eeometrischen Grundbegriffe liefert zur Kukli 
(wometrie als der im «#indeutige-metho 
lischen System allein möglichen. Beim axioma 
tischen Aufbau der Dynamik erhält er dureh wie 
derholte Anwendung «des Prinzips der Einfachstheit 
ohne Zuhilfenahme oder Ausnutzune irgendwelcher 
messenden Versuche als „Elementargestalt“ der Dy 
namik das Newtonsche Attraktionsgesetz. Mit «dem 
im Kuklidischen Raume geltenden Newtonschen 
Gesetz ist nach Dineler die Grundlage für 
systematischen Aufbau der Phvsik eine 
Grundage, die durch kein Experiment erschüttert, 
lurch keine neue Entdeekune in Frage 
werden könne. Die Verdienste, die sieh neuere 
Theorien wie Quanten- und Relativitätstheorie um 
den Fortschritt physikalischer Erkenntnis erworben 
haben. werden von Dingeler zwar zugregeben,. doch 
seien jene Theorien nur vorläufige Notbehelfe der 
vorwärtsstürmenden eehören der 
„Frontphysik” an). Die gegenwärtige Situation sei 
eben nur eine Zwischenperiode der Forschung, die 
beendet sein würde, sobald der Aufbau des „eindeu 
ir methodischen Systems” einmal ernsthaft in An 
eriff genommen würde: in diesem System würden 
dann auch die Ergebnisse jener Theorien ihren Platz 
finden, heute eeiibte Ver 
fahren, auf Grund e 


also mit 


dischen 


einen 


eeorehben. 


gestellt 


lorschun® sie 


Völlle abwerie sei das 
inzelner Versuchsereebnisse 
wie etwa des Michelsonversuchs) «die alten klassı 
sehen (d. h. Newtonschen) Grundlazen dureh 
neue ersetzen zu wollen. Bei der Darleruner und 
Verteidieune dieser Anschauungeen. «(die von ver 
schiedenen Seiten immer von neuem beleuchtet wer 
den. eeht es nieht ohne scharfe Aneriffe auf die 
herrschende „Schulmeinung“ ab: aber von den „po 
pulären Phantasien der Einstein, Planek, Heisen 
bere. Kddinzton“ zu sprechen geht doch wohl 


etwas weit. Auch hätte man an Stelle von allge 
mein «ehaltenen Andeutungen („da werden wir 
belehrt“, „es wird darauf hinzewiesen" usw.) ee 


nauere Angaben gewünscht, etwa über die Werke, 
in denen „man behauptet, Galilei habe die Galilei 
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transformationen FOSt halfen” N, 295 E Kin kleiner 
Irrtum findet sieh auf S. 168: nieht der Historiker 
les Materialismus, F. A. Lanee, sondern Ludwie 


lange war der Schöpfer des Inertialsystems. Noch 
eine Bemerkung zu dem Dinglerschen System selbst. 
Sobald! man den axiomatischen Auseangsforderun 
een Dinrlers zustimmt, wird man auch die Folge 
runzen und damit den ganzen Aufbau der Physik. 


so wie er ihn sich denkt, bejahen müssen. Ob frei 
lich «das Fundament des „eindeutig-methodischen 
Systems" wirklich jenen Grad von unmittelbarer 


Sicherheit und Ausschließlicehkeit hat. von dem er 


überzeugt ist, und ob man sein Prinzip der Ein 


leutiekeit bzw. Einfachstheit so ausleren und an 
wenden muß, wie er es vorschreibt, oder ob man 


iicht vielmehr in dem Aufbau eines eindeutigen 
Systems der Physik ein Ideal sieht, dessen Er 
reichunge sieh nur auf Grund der Fortschritte unse 
rer Erfahrungen in einem asymptotischen Prozesse 
anstreben. nicht aber lurch diktatorische Maßnah- 
men in seinen Grundlagen wie in seiner Richtung 
erzwineen läßt. das ist eine andere Frage, die hier 
nieht zur Erörterung steht. 


Berlin. K. Mosch. 840 


Dr.-Ine. ERWIN MAYER, VDI. Einfluß deı 
Querscehnittsverformung auf die Ent 
wieklune der Geschwindiekeits- und 
Druckverteilung bei turbulenten 
Strömunzeen in Rohren. (VDI-Forschungs 
heft 389, Beilare zu „Forschung auf dem Gebiete 
les Inzenieurwesens“, Ausgabe B. Bd. 9 März 
\pril 1938.) 20 S. m. 64 Abb. Berlin 1958, VDI 
Verlae G.m.b.H. Preis br. 5M. 

Ks handelt sieh hier um sehr sorgefältie durch 
eeführte praktische Untersuchungen über den Ge 
schwindiekeits- und Drucekverlauf in einem Rohr 
von konstanter Quersehnittszeröße, aber veränder 
licher Quersehnittsform (Überzanz vom Kreis auf 
das Rechteck und umezekehrt). Die Erzebnisse 
entsprechen in jeder Hinsicht dem, was auf Grund 
alleemeiner strömungesteehnischer Überlerungen 
sich einstellen muß. Eine zenaue mathematische 
Bereehnunge der Vorzängze ist bis heute noch nicht 
mörlieh. Von Wiehtirkeit ist, daß nirzends in 
ien Strömuneen Ablösuneserscheinunzen beoh 
achtet wurden, selbst dort nicht, wo die Strömung 
nicht unerheblieh verzörert wird. 

Dresden. KH. Sörensen VDI. 858 


Dr.-Ine. ALBRECHT HUSSMANN, Rechne 
rische Verfahren ur h 
\nalvse und Synthese mit Schablonen für 
eine Reehnunge mit 12. 24. 36 oder 72 Ordinaten. 
2S S, mit 24 Abb. im Text und auf 2 Tafeln sowie 
I0 Bereehnunestafeln. Berlin 1938. Verlar von 
Julius Springer. Preis steif geh. 9,60 M. 


4 l 
armonlsenen 


Für die in vielen Gebieten der Physik und 
lechnik gebrauchte harmonische Analyse enthält 
das vorlieeende Heft eine Reihe sehr brauchbarer 
Rechenblätter, die dem vor 25 Jahren von Runge 
und Emd«ı herauseereebenen Rechenformular 
ähnlich sind. Sie beschränken sich aber nicht wie 


dieses Formular und wie die bekannten Scha 
blonen von Zipperer und wie die von Tere 
besiı auf die Benutzune von 24 Ordinaten. son 


dern sehen die Verwendune von 12. 24. 36 oder 
72 Ordinaten vor. Im letzten Fall ermöglicht die 
in die Rechnune eineehende erobe Zahl von 
Ordinaten die Bestimmune von Oberschwinruneen 
höherer Ordnung, die oft von Wichtiekeit sind. 

Im einführenden Text werden die Grundlaren 
der harmonischen Analyse und Synthese ausein 
andereesetzt, ler \ufbau ler Rechenschemata 
erörtert und die Durehführunz sowohl der Ana 
Ivse wie der Synthese unter Benutzung der 
techenblätter an einem Beispiel vorgeführt; ins- 








hesondere werden die in die Formulare auf 
eenommenen Kontrollen erläutert. Weiter werden 
der Fehler des Näherungsverfahrens und Metho 
den zu seiner Abschätzung und Verbesserung be- 
sprochen. Zur Berechnung einzelner Teilwellen 
höherer Ordnung wird die Methode von Fischer- 
Hinnen empfohlen. Der Anhang enthält die 
I0 Rechenformulare, die so aufeebaut sind. daß 
tie Rechnung auch von Hilfskräften durchgeführt 
werden kann. Den Gebrauch denkt sich Verfasser 
so, dab diese Formulare mit durchsichtigem 
’apier überspannt und auf diesem an den ent- 
sprechenden Stellen die Zahlen notiert werden. 
Praktischer und übersichtlicher wäre es nach 
meiner Ansicht, wenn man die Rechnungen auf 
den Blättern selbst ausführen könnte. Um das zu 
ermöglichen, müßte sich allerdings der Verlag 
entsehließen, die Formulare einzeln zu billigen 
Preisen in den Handel zu bringen. 

Dresden. Wiıllers. 877 


B. RIEMANN, Partielle Differential 
eleiehungen und ihre Anwendungen 
aufphysikalische Fragen, herausgegeben 
von Karl Hattendorff. m. einem Vorwort von 
l’rof. Dr.-Ing. e. h. Fritz Emde (Stutteart), XI 
DD >. Braunschweie 1938. Verlae Friedr. 
Vjewee & >Sohn. Preis eeh, IH0 M. 

Nach Riemanns Tode eab Hattendorff 
unter Benutzune «einer zusammenfassenden Aus 
arbeitune Riemanns auf Grund seiner eieenen 
Nachsehrift die Vorlesune von B. Riemann über 
partielle Differentialgleicehungen, die dieser zuerst 
1854 und weiterhin noch zweimal zehalten hatte. 
heraus, Sehon in dieser ersten Bearbeitung ist wie 
in den späteren etwa die Hälfte, nämlich die drei 
ersten Abschnitte: Bestimmte Integrale, Unenidlliche 
Reihen (Fouriersche Reihen) und Differential 
eleiehunzen den mathematischen Grundlagen. die 
andere Hälfte: Bewerune der Wärme in festen Kör 
pern, Schwingungen elastischer fester Körper und 
bBeweeune der Flüssiekeiten den Anwenduneen ee 
widmet. Das Bueh erlebte in dieser Form drei Auf 
lagen, die letzte 1882. Von dieser gibt jetzt Emde 
einen unveränderten Abdruck heraus. Nach Hat 
tendorffs Tode übernahm Heinrich Weber 
die Bearbeitung. Die vierte von ihm  heraus- 
oeeebene Auflare erschien 1900/1901. Nachdem 
nachseinem Tode noeh ein unveränderter Abdruck des 
Riemann-Weber herauseekommen war, wurde 
vor etwa zehn Jahren das Werk unter Mitwirkun«e 
einer eroben Zahl von Bearbeitern vollkommen neu 
vestaltet. Dadureh wurde aber der Charakter des 
Buches «zänzlich zeändert. aus einem Lehrbuch 
wurde ein Nachschlaeewerk. Der Umfane,. der in 
er Hattendorffschen Auseabe etwas über 
300 Seiten betraren hatte und der schon in der Be 
arbeitune von Weber auf rund 1000 Seiten an 
gewachsen war, schwoll jetzt auf über 2000 Seiten 
an. Es ist daher zu beerüßen. daß das Werk in 
ler ursprünglichen Fassung wieder allgemein zu 
eänelieh gemacht ist. Wenn man natürlich heute 
auch manches anders darstellen würde, wird viel 
leicht doch mancher zern zu dieser klassischen Dar 
stellune von Riemann greifen. 
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FRANUESCO TRICOML Funzioni Anali 
tiche ‘Consielio Nazionale Delle Ricerche, Mo- 
nogratie di Matematiea Applieata, 1936/XIV). VI 
II) S. Bolorna 1936. Nicola Zanichelli Editore. 
Preis 35 Lire. 

FRANUCESCO TRICOMIL Funzioni Ellit 
tiche (Consielio Nazionale Delle Rieerche. Mono 
rrafie di Matematiea Applicata. Bologna 1937/XV, 
IN + 274 S. Boloena 1937. Nicola Zanichelli Edi 
tore. Preis 45 Lire. 
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Das Erscheinen (dieser Bücher in einer Reihe, die 
ausdrücklich der angewandten Mathematik «dienen 
soll. kann ein leiehtes Staunen hervorrufen, wenn 
man beim Durehblättern feststellt. «dab 
I) die übliehen Gerenstände der (deutschen) Univer- 
sitätsvorlesuneen über alleemeine Funktionentheorie 
umfaßt. während in 2) über 200 Seiten der Theorie 
und nur 40 den Anwendungen der elliptischen Funk 
tionen eewidmet sind. Indes zeigt eine genauere 
Prüfune, daß die Einstellunz des zweiten Werkes 
das erste soll nur eine recht gute Kinleitung 
dazu sein) zur Anwendune hin zwar nicht so sehr 
äußerlich. dafür aber um so mehr in der Ausführun« 
les klassischen Stoffes zum Ausdruck kommt: Von 
\nfane an ist mit Zeichnungen und besonders den 
Emdeschen  Reliefdarstellungen die Anschauung 
unterbaut. dann aber in aller Einzelauswahl des 
Stoffes und in der Durchführung immer zielbewubt 
der Grundsatz verfolgt, alles bis zum Beginn der 
reinen numerischen Rechnung Notwendige über 
siehtliceh bereitzustellen: Die Formeln, mit denen die 
Zurückführung auf tabulierte Funktionen gelingt. 
sind bequem greifbar gemacht: und für die Fälle, 
wo tabulierte Funktionen noch nieht zur Verfürung 
stehen, ist der Anschluß an bestkonvergente Reihen 
entwieklungen weboten. die dem Rechner mit dem 
Bleistift oder der Maschine keinen erheblichen Zeit- 
aufwand verursachen. Soweit Stichproben es zeigen 
können, ist das Buch in der Wiedergabe der 
"ormeln zuverlässie. 

In der Darstellung fällt noch auf, daß «der Verf. 
auf die praktisch wichtigen Sonderfälle der für 
reelles 3 reellwertirgen Funktionen sründlich ein 
reht und sie durch Kurvenbild und -bespreehung 
dem Leser nahebringt. Sehr verdienstlich ist das 
eute Schlagwortverzeichnis. eine breitere Formel- 
sammlung und last not least eine Übersicht über (die 
verschiedenen abweichenden Bezeichnungen der 
Verfasser von Werken über den Geeenstanld. 

So elauben wir, daß auch andere, die mit Anwen 
duneen der Mathematik zu tun haben. «das vor- 
lierenle Werk als das z. Zt. der angewandten 
Mathematik nächste über den Gerenstanıd dankbar 
anerkennen werden. 


ersten 


(sieben. K. Ullrieh. 885 

Dr. GUIDO HOHEISEL, a. 0. Prof. a. d. Univer- 
sität Greifswald. Gewöhnliche Differen- 
tialeleichungen. 3. neu bearb. Aufl. (Samm- 
lune Göschen Bd. 920 126 S. Berlin 1938. Ver- 
lae Walter de Gruyter & Co. Preis geb, 1.62 M. 

Die neue Auflage des Bändchens von Hoheisel, 
die im einzelnen an verschiedenen Stellen umge- 
arbeitet und ergänzt ist, brinet im wesentlichen den 
eleichen Stoff wie die voriee Auflage und zwar ge- 
eliedert in drei Kapitel, deren erstes den Differen- 
tialeleichungen erster Ordnung gewidmet ist. 
während das zweite kurz Gleichungen höherer Oril- 
nung und lineare Gleichungen, das dritte Randwert- 
aufeaben behandelt. In einem Anhang wird weiter 
die Abhängiekeit der Lösungen von Anfangs- 
werten und Parametern untersucht. Die Dar- 
stellunze beschränkt sieh dabei durchgehend auf 
reelle Variable. Als erste Einführunz in die Theo 
rie der gewöhnlichen Differentialeleichungen kann 
das kleine Bändehen durchaus empfohlen werden. 
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Dr. OTTO LANKES und Dr. LUDWIG BAUM- 
GARTNER, Wirklichkeit und Formel, 
einmathematisches Lesebuch für die 
Jueendunddas Volk. Tell. VIN 112 8. 
m. 3 Bildern von Josef Bergmann u. 45 Fig. Bam- 
bere 1935. €. C. Buchners Verlag. Preis broseh. 


1.40 M. Teil II. IX +146 S. m. einem Bild von 
Josef Bergmann und 46 Fig. Bamberg 1938. Preis 
brosch. 2,40 M. 


Die Verfasser sind der Meinung, daß, „wie im 
Mittelalter das Latein. so heute die mathematische 
Formelsprache die Gemeinschaft der nach Erkennt 
nis Strebenden zerklüftet“. Um Kluft zu 
überbrücken, haben sie in den beiden vorliegenden 
Bändehen „es gewagt, dies Latein der Formeln in 
unser eeliebtes Deutsch zu übertra&en“, indem sie 
„von der Wirkliehkeit her die Welt der mathema 
tischen Formeln erschließen“. Gegen dieses nicht 
verade neue Verfahren ist erundsätzlieh dureh 
aus niehts einzuwenden: im Geeenteil: die Wahl 
der Beispiele, die Reiehhaltiekeit des Materials, das 
vorsiehtiee Vorwärtsschreiten bei Einführunze schwie 
rieer Begriffe verdienen «durchaus Anerkennung, 
l,eider bevorzugen aber die Verfasser bei der Dar 
bietun® (der Beispiele Breite der Darstellune. 
die das Wesentliche zelerentlich allzusehr in den 
Hintergrund dränet, so dab mancher Leser schon 
etwas ermüdet sein könnte, ehe er an den Kern der 
Sache kommt. „Jedenfalls helfen die wortreichen 
Schilderungen in keiner Weise zur Erzielung bes 
seren Verständnisses für die mathematische 
Zeichensprache: wie zur Zeit des Euklid gibt es 
auch heute noch keinen Könieswee zur Mathematik. 
Den Verfassern war es, wie sie sagen, hauptsäch 
lieh um die Klärung der Beeriffe des Grenzwerts, 
(les Differentialguotienten und des Integrals zu tun. 
Hier fallen einige Entgleisungzen unangenehm auf, So 
wird im I. Teil S, 86 .„.beliebie eroß* und „unend 
lich eroß" als dem Sinne nach identisch hingestellt. 
Auf S. 108 wird das Leibnitzsche d.r als „der belie 
bir kleine, oder, was nichts anderes bedeutet, der 
unendlich kleine Koordinatenunterschied" definiert. 
Im 11. Teil. S. 112. heißt es sogar, dab „die noch 
ereifbaren Größen Ax, Ay erst im Dahinschwin 
den zu den Differentialen d.r, dy werden“. Das 
Ist eine Mystik, von der man in den grundlegenden 
\bhandlungen Leibniz’ etwa der „Nova Metho 
dus pro Maximis et Minimis“, wo der Differential 
beeriff klar entwickelt wird nichts findet. 
Erfreulich ist die Fülle von Anwendungen auf phy 
sikalische und technische Probleme. denen wir be 
sonders im Il. Bändchen bereenen und die didak 
tisch meist recht geschickt behandelt werden. 
Der junge Lehrer, der sich mit der Didaktik und 
Methodik des schwierizen Gebietes beschäftirt, wird 
manche Anrerune für seinen Unterricht finden. 


Mosch. 84 


diese 


eine 


Berlin. K, 


GERHARD KOWALEWSKI, Magische Qua 
drate und magische Parkette. Mit 
I0 Fir. u. einem Titelbild. Seientia Deleetans Heft 2. 
7S S. Leipzig 1937. Verlag K. F. Koehlers Antiqua 
rium. Preis brosch. 2 M. 


Den seit Jahrhunderten viel behandelten ma 
eischen Quadraten gewinnt Verfasser neue Auf 
vaben ab, indem er aus den S80 vierreihiren ma 
eischen Quadraten, die Frenicele de BressY 
als mörlich aufgezählt hat, durch besondere Beilin 
ungen Klassen aussondert. Sind die Zahlen 1-16 
in einer quadratischen Matrix a,, id —1,...4: k 
l...4) anzeordnet,. so ist die Bedingung für ein 
marisches Quadrat, daß die Summen jeder Horizon 
tal-Vertikal- und Diagonalreihe zleieh sind (— 34). 
K. schreibt noch eine Kreuzbedingzung vor: 
Aa; s. Indem die Zahlen als Ge 
wiehte gedeutet werden, ist so ein links geriehtetes 
Hauptkreuz ausgezeichnet, dessen Arme gleich 
belastet sind. Daraus ergibt sich ein rechts g» 
richtetes Geeenkreuz: Asg + Ass Agı T Ag u. 
Bei dem von Dürers „Melancholie“ bekannten Qua 
drat ist die Kreuzbedineune erfüllt und außerdem 
s Ss” 17. Marische Quadrate dieser Eigenschaft 
nennt Verf Dürerquadrate zum Unterschier 
von den Sonnen dquadraten, bei denen zwischen 
Belastung des einen und anderen Kreuzes der gröbt 
mögliche Unterschied 48 besteht. 


Ayo T (ya (Ion 








/. rer Untersu na sowie ler 48 Düre 
nad e We lie Zahlen 1—16 dur: 0—15 ı 

setzt un 4 IArweste lie Matrix deı 
l ‚ | It \ | rd IM le1 « I Oben 
) us \ lanı Pa ette eebilck irre 
ll sind nieht parkettbildend: von den Sonne 

‘4 | le VIer. \erf. kommt ı einem neuartieen 
’arkett aus dunklen Quadraten und hellen Rhon 

S e] E \ itekten empfie Er Das ebene 

I) N ra I | Like Its t eil 


Dr. HANS KELLERER, Mathematik und 
id Kkı4 r, eine Jlebensnahe Kinführung 
lie Methoden der Statistik. (Math. 
: hrse. v. W. Lietzmann u. 
\. Wittine. Bd. 094 IS S. m. 20 Abb. und 20 Aufe. 


B. (1, Teubner. Preis 


\n Hand von Beispielen, die den verschiedensten 
(ebjeten des \Verkehrslebens entnommen sind, führt 
Verfassen | einfachsten Rechenverfahren 

| Darstellunesformen der Statistik ein. Er glie 
len Stoff in drei Absehnitte. Im ersten be 


nlelt er Mittelwerte. Streuunesmabe und Ver 


niszahlen. im zweiten Häufiekeits- und Summen 

\ Häufiekeits-, Summen- und Konzentrations 

ven, Belastunesbänder und Streuungesbilder uni 

ı letzten statistische Bearbeitune zeitlich ze 

ordneter Zahlenreihen. Da die Darstellune eanz 

enta ılten ist. eienet sieh das Bändchen 

für eine erste Einführune in «las Gebiet der 
Statistik. 
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U, PESONEN, Messune der Basis Bala 
schov in Rußland, im Jahre 1935 (Baltische 
Sonderveröffentlichung 


’. 16 8. Helsinski 1938. Verla@ Osakevhtiö 


reed liitiIseile Kommission 


Geodätischen 


Kommission beivetreten \\ LT. wi le au! \nreeunge 


I 
von Il. Bons«dorff in dem baltischen Dreiecks 
olvLron len bereits gemessenen sieben 
(rundlinien nor eine weitere, auf russischem Ge 
let, mM ler Näln von Balaschov oplerene (‚rund 
]t FOTOS Die \essune le nelien (Grundlinie 
t ) te Ilırı lest } Beohba: ie] V. R. (0) an leı 
ınd U, Pesonen sowie mit (denselben Geräten 
wie die der älteren Grundlinien. Zur Verwendung 
kamen acht, je 24 m Jaı Invardrähte des Finni 
N n Geodätisehen Instituts, von denen je vier in 
verschiedenen ‚Ja erzestt wurden, so «ab 
o renänderungeen zu erwarteı 
vareı I) (‚rund ENMESSUNg"“ mit Invar 
' Su u, tire Bestimmun ler Dral 
ins olet vor und r Messune auf Grund 
li & strecke ei Nummela leren 
N] 65 Pra O eirae \uberden wurden 
) t Lil vi] nd na | VMessu or ie] 
ru | / lt I SSISt 1] \o ılmab n Mos 
I verg en. Diese loppelte | ınz ergal 
It ax X | ersehleid visehe le nnischen 
nd (den SSıS l.äng nheit Die Messung 
Tr (rund eriolete ler be Messu | 11 
In | \nordnung Der tier: 
Fi t l | re « PIWiAs 4 \ 0) N no 
(‚rund 0 1.9 
Stu l |’ Wi IK ( Ss Tte7I N72 


tia Kauations AVI+416 38. 
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Ferner sind bei der Schriftleitung foleende Bücher 


eingegangen (ausführliche Besprechung bleibt vor- 


\meriean Mathematical Soeletv. 
Semiecentennial Publieations Bd. 1 
RAYMOND CLARE ARCHIBALD, A Semicen 
tennıal History ofthe American Ma 
thematieal Society 1888-1938, with Bio 
eraphiesofthe Past Presidents ?2#N. 


Bd. II, Semicentennial Addresses of 


ıe American Mathematical Societv. 
»15 S. New York 1938. Verlae American Mathema 


tie] Soeletv, 


\W,M. SMART, M. AX.. D. Se.. Reeius Professor of 
\stronom\ in the Universit\ of (laseow. 
ar Dvnamies, VII+ 434 8. Cambridee 1938. 


"he University Preis 30 Schilling. 

REGLER. Höh. 
| lehranstalt Augsburg, Physikalische 
\ufeabensammlune für höhere Schule 
| leipzig 1939. Verla® Jo 


‘ 1 
Stel 


Pı emn 


JOSEF Studienprofessor a. d. 


eehn. 
ınd Studium. 1209 8. 
hann Ambrosius Barth. Preis 6.30 M. 


NA\RI,R MARBE. Das \useleichsprinzip 
n der Statistik und verwandte Pro 
bleme. \ 164 S. München 1938. G. H. Beck 


handlune, Preis broseh. SM. 


GEORG SCHEFFERS, Lehrbuch 
ler Mathematik zum Selbstunterrieht 
ınd für Studierende \ 
und der Technik. 
führunein di« 


seI)ıe \erlaesbuceh 


se] a ftten 


0 Differential- und Inte 
eralreehnung und in die analvtische 
Geometrie. VIIT+743 S. m. 438 Fig. Ber 


| 
Iın 1938. Verlae Walter de Gruvter & (Üo, 


Dr. phil. OTTO HAUPT, Prof. a. d. Univ, Erlan 
een, unter Mitarbeit von Dr. phil. GEORG AUMANN, 
Prof. a. d. Univ. Frankfurt a. M.. Differential 


und Inteeralreehnune unter beson 


lerer Berücksichtieune neuerer Er 
vebniss« (öschens  Lehrbücherei. Bd. 24). 
I. Bl. Einführung in die reelle Analx 
Sıs 196 S. m. 2 Fie. Berlin 1938. Verlae Walter 


ru Preis geb. 11.20 M. 

Del. Bd. 25. IH. Bd. Differentialrech 
nune 186. m. 1 Fie. Berlin 1938. Verlae Walter 
le Gruyter & Co, Preis zeb. 9,80 M. 

Del. Bl. 26. II. Bd. 
in 1938, Verlag 
10.60 M. 


Dr. teelın. ADOLF PUCHER, Die Momenten 
ınflußfelder rechteekieer Platten 
„Deutscher Ausschuß für Eisenbeton“, Heft 90). VI 
Berlin 1938. Verlae Wilhelm Ernst & Sohn. 
Preis 5.50.M. 

R. A. FRAZER, PB.A. D.Se. F.R.Ae.S, of the 
\erodlvnamies Department. the National Physical 
‚aboratorv., 

W. J. DUNCAN, D. Se.. A.M. I. Mech. E.. F. 
Ss, Wakefield Professor of Aeronauties at the Uni 
versitv of Hull. 

\. R. COLLAR, B. \.. B. Se. of the Aerodvnamies 
Department, the National Physical Laboratorv. 
KBlementary Matricesand Some Appli 

ations to Dvnamies and Differen 
| ) Cambridge 


Inteeralreehnune., 
Walter de Gruyter & 


AS N 


| 


i 


1938. The University Press, Preis zeb. 30 Schilling. 


T. G. ROOM, M. \.. 
‘ University of Sidney, The 
Determinantal l,oci 
| Universit\ 


Professor of Mathematies in 
‚eometr\ of 
XXVII +43 8. 


Press. Preis eb. 
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NACHRICHTEN 


Arnold Sommerfeld zum 70. Geburtstag. 
Wochen 


beeineen die deutschen 


Geburtstages 


Vor wenieen 
Physiker die Feier des siebzigsten 
ihrer erößten Theoretiker, A. Sommer 
feld. Auch die Vertreter der Anzewandten Me 
ehanik haben allen Grund, bei Anlaß mit 
Dankbarkeit des Mannes zu gedenken. der in 
jungen Jahren einmal einer der ihrigen war. >Som- 
merfell war nach seiner Studentenzeit Assistent 
von Felix Klein in Göttingen und wurde von 
liesem mit der Ausarbeitune und «dem weiteren 
\usbau einer Vorlesung über die Theorie des 
Kreisels betraut. Daraus ist das bekannte Kreisel- 
buch entstanden !),. Es war wohl wesentlich (dieses 
schöne Buch, das ihm. nachdem er ein paar Jahre 
die Studenten der Bereakademie Clausthal Mathe 
matik zelehrt hatte, 1900 den Ruf auf den Me- 
ehaniklehrstuhl der Teehnischen Hochschule Aachen 
einbrachte. Wie sehr er sieh mit der Schaffenslust. 
die seine eanze Lebensarbeit kennzeichnet, in diese 
ihm neue Aufzrabe vertieft hat. zeigen seine Ver 
öffentliehungen aus dieser Zeit. Am klarsten geht 
seine Auffassune von den Aufgaben der techni- 
schen Mechanik aus dem zusammenfassenden Be 
richt hervor, den er auf der Naturforscherversamm- 
lune in Kassel 1903 erstattet hat?).  Liebevolle 
Versenkune in eine Einzelaufgeabe zeigt eine AÄr- 
beit ..Zur Theorie der Eisenbahnbremsen“ 3), in der 
der Bremsvoreane bei Bremsdruck für 
ein experimentell eeschwindiekeits 
abhäneiees Reibungesgesetz in allen Einzelheiten 
durehdiskutiert ist und die Bremswege in Ah 
häneiekeit vom Bremsdruck für normale Bremsung 
und auch für Übereane zu fest gebremsten, auf den 
Schienen eleitenden Rädern numerisch verfolgt 
ist. Von seinen Ideenzeänzeen und seinen anschau- 
lichen Demonstrationsversuchen aus dieser Zeit 
iindet man vieles in: den Berichten des Aachener 
Bezirksvereins Deutscher Ingenieure. Wegweisenil 
zeigt sich hier sein Vortrag „Zum dynamischen 
Ausbau der Festiekeitslehre“ (3. Juli 1901). in dem 
er sich vor allem mit Schwingungserscheinungen 
befaßt und u. a. seinen „Wackeltisch”“ demon- 
striert. einen Tisch. dessen Platte dureh einen mit 
einem exzentrischen Gewicht versehenen Elektro 
motor in Resonanzschwingungen versetzt wird und 
durch eine Rückwirkung auf den Schwingungs 
erreeer die Drehzahl des Motors in sehr auffälliger 
Weise beeinflußt. Ein anderer Vortrag befaßt sieh 
nit der Bestimmung der Knicklast Stabes 
durch Beobachtung Eigenschwingungszahl 
bei steieender Last. Eine sehr schöne, zum Teil 
allerdings auf das physikalische Gebiet  zielende 
Arbeit befaßt sich mit Lissajous-Figuren und Re 
sonanzschwinzungen bei Schraubenfedern und ihre 
Verwertune zur Bestimmung des 
Verhältnisses %#). Eine Arbeit aus der Elastizitäts 
theorie, die zwar das erstellte Problem nicht voll 
erleliete. damals als erste Arbeit auf diesem Gebiet 
jedoch sehr viel Beachtung fand, befaßt sich mit der 
Knieksicherheit der Stege von Walzwerkprofilen >). 


elNeSs 


diesem 


eeerebenem 


festeestelltes. 


eines 
seiner 


Poissonschen 


I, F. Klein und A. Sommerfeld: Über die Theorie 
des Kreisels, Heft I (Einführung) 1897, Heft II (schwerer 
Kreisel) 1808, Heft III (störende Einflüsse, astronomische 
und geophysikalische Anwendungen) 1903, Heft IV (teeh- 
nische Anwendungen, mit F. Noether) 1910, 


2), Die naturwissenschaftlicehen Ergebnisse und die Ziele 
der modernen teehnischen Mechanik, abgedruckt in Phys. 
2L., Bd. 4 (1903), S. 773 bis 782. 


3), Denksehrift der Technischen Hochschule zu Aachen, 
verfaßt aus Anlaß der Industrie-A usstellung in Düsseldorf 
1902, S. 58 bis 71. 


4), Festschrift Adolf Wüllner gewidmet zum siebzig 
sten Geburtstag, Leipzig 1905, S. 162 bis 193. 


5, Z. f. Mathematik und Physik, Bd. 54 


(1906), 8. 
153 und 318 bis 324. 


113 bis 


Arbeit zur technischen 


Zeit, 


\uech ein 
stammt aus 


Hyidrodynamik 
vielbeachtete und 


dieser seine 


auch heute noch viel erwähnte Rechnung zur hyıdro 


ıyvnamischen Theorie der Schmiermittelreibung ®). 


unreeht. darüber zu klaren. 


nach 


Ks wäre natürlich 
dab Sommerfeld schon sechsjähriger Tätig 
keit in Aachen dureh Berufune auf die Professun 
für theoretische Physik an die Universität München. 
wo ihn Fülle von Aufeaben 
technischen Mechanik fast eanz verlorengerangen 
ist. Denn was er in der Physik geleistet hat. ist 
so eewaltie, daß wir von der Mechanik uns fürliech 
mit dem Verlust abfinden müssen, den unsere Gilde 
durch seinen Forteane erfahren hat. Das Inter 
esse für (die Fraren der Mechanik une speziell der 
Hydrodynamik hat aber auch weiter bei 
feildl noch nacheewirkt. Seine auf dem Mathe 
matikerkonereß in Rom 1008 voreetraoene Xrbeit 
„Ein Beitrag zur hvdrodynamischen Erklärung (der 
turbulenten Flüssiekeitsbeweruneen“ "). an die sieh 
noch einige hydrodynamisehe 
sehlossen,. ist eine Frucht 
seiner Technikzeit. 


eine erwartete, «der 


Sommer 


Dissertationen an 


dieser Nachwirkung 


beereiflieh. daß er bei der 
Entwicklung, die die Physik 
hat. eanz von ihr absorbiert 


is ist eewaltigeen 
seitlem genommen 


worden ist und wir 


nicht hoffen dürfen, daß er noch einmal zu den 
Problemen der Mechanik zurückkehrt. Wir wün 
schen ihm. daß er noch auf viele Jahre seine 
Schaffenskraft zum Besten der deutschen Wissen 


schaft bewahren möchte und füren für uns den be 
Wunsch zu. daß er die Entwieklunge., 
ach seit seinem Ausscheiden aus ihm 
eenommen hat, weniestens von ferne mit dem alten 
Interesse verfoleen 


scheildlenen 


die unser 


more, 


l.. Prandtl 9% 


Stephan Timoschenko zum 60. Geburtstag. 


Ks eibt in der Mechanikliteratur so manche Ar 
beit. lie rein mathematische Ziele verfolgt und die 
mechanische Fraeestellune nur als schmückenden 
Rahmen benutzt. und auf der andern Seite auch 
manche, die mit unzureichenden mathematischen 
Mitteln an ein Problem herangeht und kaum mehr 
als eine für den Augenbliek zeeschaffene Notlösung 
liefert. Zwischen diesen beiden Extremen liert «das 
fruchtbare Gebiet technisch-wissenschaftlicher For 
schung. «+ 


In Timosehenko. (dem wir heute unsere Glück 
wünsche zu seinem 60. Geburtstae brineen, ver 
ehren wir einen wissenschaftliehen Führer, der ein 
Menschenleben Jane Klastizitätstheorie in 
fruchtbaren Sinne zetrieben hat und einer 
ration junger Ingenieur-Forscher den Wer gewiesen 
hat. der zwisehen den beiden Extremen hindureh 
zum Erfole führt. 


Am Anfange seiner Forschertätickeit 
nur wenie bekanntzewordenen russischen Arbeiten. 
in denen schon die Fraren der elastischen Stabi 
lität die Hauptrolle spielen, und die schon der Ent 
wieklune der energetischen Methode der Stabilitäts 
theorie gewidmet sind. In den letzten Jahren vor 
lem Krieee wandte er sich dureh Aufsätze in 

französischen Zeitschriften an eine 


diesem 


(ene 


stehen (die 


deutschen und | 
erößere Öffentlichkeit und gab 1013 in den Annales 


des Ponts et Chausseces einen ersten zusammen 


6), Z. f. Mathematik und Physik, Bd. 50 (1904), S. 97 bis 
155, wieder abgedruekt in Ostwalds Klassikern, Heft 218, 
S. 108 bis 201. 

7”, Attidel IV. Congresso Internazionale dei Matermatiei, 
Rom 1909, S. 116 bis 124. 
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lassen len Bericht übeı die eesuamte damaliee Sta 
hilitätstheorie. 

\ls er nach dem Kriege, bereits ein Mann von 
Weltruf, nach Amerika übersiedelte, fand er eine 
nNelle \ufeabe vor sieh: lie (rün lune eines eroben 
Schülerkreises, Wie befruchtend er damit auf die Ent 
vieklune der Elastizitätstheorie in Amerika einzewirkt 
hat, erkennt man daran. daß der überwierende Teil 
ier junzen amerikanischen Elastiker seine Schüler 
sind. Neben vielen Einzeluntersuchungen und den 


Beiträren zum Handbuch der Physikalischen und 
lsehnisehen Mechanik hat uns Timoschenko in 
lieser Zeit vor allem seine l,ehrbücheı eeschenkt. 
le ihren Wee lureh lie vanze Welt eemacht 
haben und vielfach übersetzt worden sind. Mehrere 
befassen sieh mit der Festiekeitslehre im Ganzen, 
Kins ist «der technischen Scehwingeungslehre ze 
widmet une elot die oleiche Tendenz wie SEINE 
\rbeiten zur >Stabilitätstheorie: die Entwieklung 


der auf energetischen Betrachtungen aufgebauten 
Näherungsverfahren, die dem Ingenieur ein Mittel 
in lie Hand vzeben. auch in schwierigen Fällen 
schnell, anschaulich und dabei auf zuverlässigem 
Boden zu einem Zahlenerzeebnis zu kommen. Sein 
letztes eroßes Werk ist wieder eine zusammen 
tassende Darstellun®e der Stabilitätstheorie. Ver 
eleieht man sie mit der französischen Arbeit von 
i013. so ermißt man, welche Entwieklunge auf 
liesem Gebiet inzwischen statteefunden hat, und 
wer (die Literatur im einzelnen kennt. weib,. einen 
vie eroßen Anteil Timoschenko selbst. seine Schüleı 
und die (dureh seine Arbeiten aneereeten Forscher 
aller Völker «daran gehabt haben. 
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Internationaler Mathematiker-Kongreß 1940. 


Der nächste internationale Mathematiker-Kongreb 
iindet in der Zeit vom 4. bis zum 12. September 
1940 in Cambridre Massachusetts USA, statt. 
\ußer eingehenden Besprechungen über ‚Algebra, 
Maßtheorie, Wahrscheinliehkeitsreehnung, Mathe 
natische Logik und Topologie werden kürzere 
Vorträre z„ehalten werden. Für diese sind zur 
Zeit sechs Sektionen voreesehen: 1. \leebra und 
>. Analvsis, 3. Geometrie und Topo 


lorie, {. Wahrscheinliehkeitsreehnung, Statistik. 


Zahlentheorie. 


Versicherunesmathematik und  Volkswirtschafts 
lehre, 5. Mathematische Physik und angewandte 
Mathematik. 6. Logik. Philosophie, Geschichte unıd 
Pädaeoeik. Nähere Auskunft über Programm. 
Unterkunft. Geselliere Veranstaltungen usw, erteilt: 
Fhe \n eriean Mathematical Soclet\ >31 West 
116th Street. New York City. UN.A. 220 


Neu berechnete Tafileln ! 


FE. L. INCE, Tables of the Elliptie-eylinder Fune 


tions: and Zeros anı Turninz Points of the Ellip 


tie-evlinder Funetions. Proc. R. Soe. of Edinburgh 
1931- 32), vol. LIIL, part IV (Nos. 22, 23), p. 359... 
133. Als Sonderdruck käuflich. Preis: 7 Shilling. 


die Tafeln betreffen (die 12 periodischen lL,ösuneen 


nn > A \ (a we » di) > » =) 
Ce&,n (7.9), SontılT.F9) Ckan+ı!lT, ). S@on I 


d? 7) 


für n=0.1.2 der Differentialrleichung 2 
d .r: 


1) Kiinftie wird hier über nen berechnete Tafeln höherer 


Funktionen berichtet werden. 
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(a —29cos2r) y=0 „Mathieusche Funk- 
tionen”), die gemäß |y?dr sc, also auf den 
() 


quadratischen Mittelwert Y 0,5 normiert sind. 


Für =0(0,1,2... 910,12... 18,20,24... 36,40 
werden gegeben: 


l. Die „charakteristischen Zahlen” a, für ce, 
und d,, für se,„. bei denen eine periodische Lösung 


vorhanden ist. mit 7 Dezimalen. (Für © 0) ist 
dm Din m“, Für 9 I) ist U Fr b, > 67.6: 
As 2 b, = — 1.73: lb, —_ —- 41.43.) 


) 


2. Die Fourierschen Koeffizienten 4A, für ce, 
und 7, für se,, mit 7 Dezimalen. soweit sie in der 
‘. Dezimale von 0 verschieden sind (k < 24). 


3. Die Nullstellen z und Extremstellen z in 
(‚raden mit 4 Dezimalen. (Für se, bei © 10: 
Nullstellen = 545° und 734°: Extremstellen = 
10,5°: 64,3: 81,8 °,) 


Ferner werden für © 1.2, 2... 9,10 und z==0°., 
1" ...89°, 90° die Werte der 12 ersten Mathieu 
schen Funktionen mit 5 Dezimalen und ihre ?. 
Differenzen gereben. 


m) 


ce; | 09-10 
aan EEE Im Aw | . 
ce 


| } 1 156; ! 156, | 
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Bild 1. 


Um eine Vorstellune von dem Ganze der 
Mathieuschen Funktionen zu zeben. sind in 
Bild lin dem Bereich £r=0...9%° die Kurven für 
Ce,. Se,. Ce,. Se, gezeichnet. und zwar bei © 10. 


Die Interpolation nach © bei festem .r stößt noch 
auf Schwieriekeiten. 


Für die Bereehnunz wurden folrende 4 Maschinen 
lenutzt: eine Brunsvira nova III. eine Trinks-Tri 
plex-Brunsviga. eine Facit und «eine Sundstram 
Buchungesmaschine. Die Bereehnunesarbeit hat sich 
anscheinend etwa über fünf Jahre erstreckt. 


Fritz Emde. 928 


Persönliches. 


Der 0. Prof. für Bauelemente der Luftfahrzeure 
an der Techn. Hochsehnule Berlin Dr. Herbert Wae- 
ner trat vom Lehramt zurück. 


l’rof. Dr. Hamel. Berlin. wurde zum ordentlichen 
Mitelield der Preußischen Akademie der Wissen- 
schaften ernannt. 

Der n. b. a. o. Prof. an der Techn. Hochschule 
harlsruhe Dr. Wolfeane Haack wurde zum o, Prof. 
für Mathematik an dieser Hochschule ernannt. 





Für den Textteil verantwortlich: Professor Dr. Fr. A.Willers, Dresden-A.20, Dorotheenstr. 12. — Printed in Germany. 


Druck von A. W. Ziekfeldt, Osterwieck am Harz. 
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Sonderheft Flugabwehr 


Um einen Überblick über den Stand und die Entwicklungsrichtung 
der Flugabwehrtechnik zu geben, wurden die in den letzten Jahren 
in der „Zeitschrift des Vereines deutscher Ingenieure” und der 
„Rundschau Deutscher Technik” von A. Kuhlenkamp VDI, Flieger- 
Stabsingenieur im Reichsluftfahrtministerium, veröffentlichten Auf- 
sätze zusammengefaßt in dem 


Sonderheft Flugabwehr. Herausgegeben vom Verein deurcher 
Ingenieure. 1938. DIN A4. 52 Seiten mit 108 Abbildungen. 
Broschiert RM. 2,75 (VDI-Mitglieder RM. 2,50). 





Inhalt: Geräte und Waffen der neuzeitlichen Flugabwehr 7 Physikalische und 

technische Grundlagen der Abhörgeräte für den Flugzeugsähall / Scheinwerfer 

der deutschen Luftwaffe 7 Die Feuerleitung von Flakbatterien 7 Die Getriebe in 

Feuerleitgeräten 7 Das Flakkommandogerät „Vickers“ / Visiere für die Flug- 

abwehr 7 Flugabwehr-Kanonen / leichte Flugabwehr-Kanonen / Sachwort- 
verzeichnis. | | 
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